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ABSTRACT 


The  microwave  properties  of  partially  magnetized  materials  have 
been  investigated.  The  real  parts  (p',  k'  ,  p' )  of  the  components  of  the 
permeability  tensor  have  been  measured  as  a  function  of  the  ratios  of  satu¬ 
ration  magnetization  (4irM  )  to  angular  frequency  (u)  and  average  magneti- 
zation  (4rrM)  to  angular  frequency  and  location  on  the  magnetic  hysteresis 
loop  on  an  yttrium  iron  garnet  (YIG)  and  magnesium  manganese  ferrite 
(MgMnF).  The  off-diagonal  component,  gives  good  agreement  with 
Rado's  prediction  that  k'  ^yA-irM/u.  The  diagonal  components  p',  p'  do 
not  fit  Rado’s  prediction,  and  the  cause  of  the  failure  of  Rado's  analysis 
forp',  p'  is  suggested.  The  dependence  of  p\  p'  upon  •y4ttM/u),  y4ttM  /u> 
has  been  empirically  fitted  by  certain  mathematical  functions.  The  im¬ 
aginary  components  (p",  p:')  have  been  measured  for  0.  y4irM  /uS  1 

z  s 

using  circularly  polarized  cavities.  For  -y4itM  /<j  <  0.  8,  it  has  been  found 

s 

that  k"  «  p,;  and  that  the  loss  of  the  partially  magnetized  state  can  be 

characterized  by  a  single  value  of  p'\  that  for  the  completely  demagnetized 

state.  For  values  of  y4irM  /u*  0.7  large  rectangular  samples  and  exact 

s 

the  solution  to  the  waveguide  cavity  mode  have  been  utilized  to  obtain  p" 
for  the  demagnetized  state. 

The  high  power  threshold  has  been  measured  on  YIG  as  a  function 
of  sample  shape,  the  polarization  of  the  rf  magnetic  field,  and  the  orientation 
of  the  rf  field  to  the  static  magnetization.  The  Suhl  and  Schlomann  theories 
for  high  power  nonlinearity  have  been  extended  to  include  arbitrary  ellip¬ 
soids,  linearly  polarized  rf  magnetic  field,  anti-Larmor  circularly  polarized 
rf  magnetic  field,  and  arbitrary  orientation  between  the  rf  magnetic  field 
and  the  static  magnetization.  The  experimental  data  and  the  extended  theory 
give  good  agreement  for  saturated  samples.  For  thin  disk  geometries,  it 
has  been  found  that  the  threshold  is  very  critical  to  alignment  when  the  dc 
magnetic  field  is  perpendicular  to  the  plane  of  the  disk.  If  a  needle-like 
anti- parallel  domain  structure  is  assumed  for  the  demagnetized  state,  then 
the  experimental  results  for  the  demagnetized  state  can  be  explained  in 
terms  of  the  results  of  the  saturated  state. 
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I.  GENERAL  INTRODUCTION 


Great  impetus  has  been  given  to  ferrite  device  development  in 
recent  years  by  the  phased  array  radar  development.  The  average  phased 
array  antenna  requires  several  thousand  phase- shifting  elements.  These 
elements  must  be  cheap,  reliable,  all  identical,  and  of  course  have  ac¬ 
ceptable  operating  characteristics  such  as  low  insertion  loss ,  a  certain 
peak  power  capability,  low  drive  power  for  phase  setting,  rapid  phase 
switching,  bandwidth,  and  a  capability  to  operate  over  a  wide  temperature 
environment.  The  drive  power  requirement  has  essentially  limited  con¬ 
sideration  to  those  devices  where  the  magnetic  material  is  partially  mag¬ 
netized.  Devices  in  which  the  ferrite  is  latched  into  the  remanent  state  have 
been  particularly  popular.  These  devices  have  been  waveguide,  coax,  strip¬ 
line,  and,  most  recently,  microstrip  geometries.  This  latter  geometry  is 
receiving  a  great  deal  of  attention  today  because  of  its  application  in  micro- 
wave  integration  and  miniaturization. 

The  requirements  on  these  devices  offer  a  challenge  to  both  the 
microwave  designer  and  materials  physicist.  Results  to  date  indicate  that 
there  is  a  real  need  for  clever  design  and  for  better  materials.  The  re¬ 
sults  also  indicate  there  is  a  need  for  a  better  understanding  of  how  material 
characteristics  affect  device  performance.  But  this  is  not  possible  if  the 
material  characteristics  of  partially  magnetized  materials  are  not  known. 
Therefore,  the  goal  of  our  program  is  to  experimentally  determine,  and 
when  possible  to  create  a  theoretical  model  of  the  significant  microwave 
characteristics  of  partially  magnetized  materials. 

A  question  frequently  posed  by  the  microwave  designer  to  the 

materials  physicist  is:  what  is  the  optimum  ratio  of  saturation  magnetization 

(4rrM  )  to  operating  frequency  (<j)  to  obtain  the  minimum  insertion  loss?  Or, 

in  general,  how  does  insertion  loss  vary  with  y4ttM  /u>  in  a  360°  phase  shifter? 

s 

We  have  a  vague  idea  how  this  may  vary  in  some  nonreciprocal  geometries 
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but  no  inkling  of  how  it  may  vary  in  a  reciprocal  geometry.  This  is  a  par¬ 
ticularly  difficult  question  since  changing  y4TrMg  means  changing  composition 
and  ceramic  microstructure.  These  are  two  variables  to  which  magnetic 
loss  is  very  sensitive.  Consequently,  in  this  program  we  have  put  a  strong 
emphasis  on  temperature  and  frequency  measurements  in  which  composition 
and  ceramic  microstructure  are  held  constant. 

Designing  for  low  insertion  loss  is  further  complicated  by  peak 
power  requirements.  Since  peak  power  threshold  and  insertion  loss  are 
dependent  upon  similar  relaxation  mechanisms,  low  insertion  loss  usually 
means  low  peak  power  performance.  Techniques  for  improving  the  peak 
power  capability  of  a  material  frequently  result  in  an  increase  of  insertion 
loss.  Therefore,  it  is  important  to  understand  the  source  of  the  peak- power 
nonlinearity  in  each  type  of  device  structure.  To  do  this  we  must  first  have 
some  idea  of  the  field  distribution  in  the  device  and  second  must  understand 
how  the  nonlinearity  threshold  for  partially  magnetized  materials  depends 
upon  the  type  of  rf  polarization  (linear,  elliptical,  circular),  the  angle  be¬ 
tween  the  rf  field  and  the  dc  magnetization  and  the  sample  geometry.  Clearly 
device  configurations  that  locate  material  in  regions  of  low  peak  power  thresh¬ 
old  and  low  phase  shift  activity  are  undesirable.  Such  a  situation  can  occur 
in  latching  phase  shifters  where  it  is  necessary  to  include  magnetic  material 
within  the  microwave  structure  for  the  purposes  of  providing  a  closed  dc 
magnetic  path.  Such  material  can  contrioute  little  to  the  phase- shifting 
properties  but  can  be  the  source  of  the  peak  power  instability 

An  eventual  goal  of  this  progr  am  is  a  handbook  of  the  microwave 
properties  of  partially  magnetized  materials  which  could  be  used  in  computer- 
aided  device  design  Theory  now  exists  that  permits  some  such  design  on 
nonreciprocal  waveguide  phase  shifters  and  junction  circulators.  Computer 
design  of  these  geometries  and  other  geometries  for  which  a  suitable  theo¬ 
retical  approach  is  developed  can  optimize  performance  more  quickly  and 
cheaply  than  the  present  laboratory  approach.  In  such  a  handbook,  it  would 
be  desirable  to  classify  the  phase  shifting,  insertion  loss,  and  peak  power  per¬ 
formance  as  a  function  of  all  the  compositional  and  microstructure  parameters. 
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Unfortunately,  our  present  understanding  of  the  effect  of  microstructure 
and  composition  upon  loss  and  critical  power  level  is  not  sufficiently  well 
developed  to  allow  such  a  classification  within  the  resources  of  the  present 
program.  A  more  realistic  approach  is  to  measure  the  phase  shifting,  loss, 
and  critical  power  characteristics  of  the  most  popular  commerical  phase 
shifting  materials.  Then  the  handbook  will  permit  computer-aided  design 
with  today's  materials.  Such  design  will  be  a  valuable  feedback  for  future 
material  development. 

This  report  describes  the  progress  during  the  first  year  of  this 
program  and  is  divided  into  two  parts.  The  first  part  deals  with  the  low  power 
determination  of  the  real  and  imaginary  parts  of  the  components  of  the  per¬ 
meability  tensor,  the  second  with  high  power  measurements  as  a  function  of 
geometry,  polarization,  and  pumping  angle.  In  both  phases  of  our  investiga¬ 
tion  established  experimental  techniques  were  employed.  However,  certain 
modifications  and  innovations  were  required  either  because  the  material  was 
partially  magnetized  or  because  of  the  novel  geometry.  Since  these  modi¬ 
fications  and  innovations  comprised  a  significant  portion  of  our  effort  during 
the  first  year  of  this  program,  they  are  described  in  detail.  The  remainder 
of  our  effort  involved  measurements  on  three  materials:  G-113  (YIG), 

TT1-390  (MgMnF),  and  MMF-10  (MgMnF).  In  a  certain  sense  these  three 
materials  have  been  trial  horses.  The  usual  problems  of  drift,  temperature 
control,  sensitivity,  reproducibility,  etc.  ,  have  been  worked  out  in  the  process 
of  measuring  these  materials. 


II.  PERMEABILITY  TENSOR  OF  PARTIALLY  MAGNETIZED  MAGNETIC 
MATERIALS 


A.  Introduction 

Much  work  has  been  done  in  recent  years  to  develop  microwave 
devices  such  as  phase  shifters  and  circulators  that  operate  in  partially 
magnetized  states.  These  states  include  all  those  from  the  completely 
demagnetized  state  to  the  remanent  state  and  even  slightly  higher  using 
small  additional  bias  fields.  In  addition  to  the  experimental  investigations 
of  these  devices,  various  geometries  have  been  theoretically*  examined 
anJ  computer  programs  exist  which  describe  the  performance  of  the  device 
in  terms  of  the  componepts  of  the  permeability  tensor  of  the  ferrimagnetic 
materials  in  them.  Unfortunately,  information^  available  in  the  literature 
describes  only  in  broad  generalities  how  these  tensor  components  depend 
“  on  the  state  of  magnetization,  and  the  saturation  magnetization  of  the  ma¬ 
terial  and  the  operating  frequency.  However,  these  computer  programs 
need  numerical  values,  preferably  in  a  functional  form,  so  that  numerous 
cases  can  be  evaluated  and  device  geometry  and  magnetic  materials  can 
be  optimized.  Thus  one  of  the  objectives  of  this  program  is  to  describe 
quantitatively  how  the  permeability  tensor,  including  both  its  real  and 
imaginary  parts,  varies  as  one  moves  about  the  magnetic  hysteresis  loop 
and  as  one  varies  frequency  and  saturation  magnetization. 

When  a  magnetic  material  is  placed  in  a  small  rf  magnetic  field 
an  rf  magnetization,  4irrnr,  is  produced,  having  components  related  to  those 
of  the  driving  rf  magnetic  field,  "R,  by  a  tensor  permeabilty  fT 

4Ttm  =  (pi  -  1)  •  H 
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Since  a  polycrystalline  material  exhibits  rotational  symmetry  about  the 
direction  of  magnetization,  which  we  ma;>  take  as  the  z  direction,  this 
tensor  takes  the  simpler  form: 


having  only  three  independent  quantities  p,  k,  and  r  .  Each  of  the  com- 

z 

ponents  of  these  tensors  will,  in  general,  be  complex: 


R  =  R*  -  JR" 

K  =  K1  -  j K" 


R 


z 


JR 


If 

z 


The  real  parts  of  these  components  produce  dispersion,  while  the  imaginary 
parts  produce  absorption. 

If  the  material  is  saturated  by  a  variable  internal  magnetic  field  H 

in  the  z  direction,  the  elements  p  and  «  undergo  a  resonance  at  a  field  u/v 

where  y  is  approximately  2.  8  Oe/MHz.  The  line  shape  of  this  resonance 
3 

has  been  analyzed  and  is  outside  the  scope  of  this  work.  We  are  interested 
only  in  the  region  where  H  approaches  zero  and  the  sample  is  not  completely 
magnetized.  If  the  internal  magnetic  field  could  be  reduced  to  exactly  zero 
without  any  demagnetization  taking  place,  the  theories  applicable  to  saturated 
materials  would  still  be  applicable  and  we  would  have 


r'  =  1 

p'  =  1 

z 

*'  =  where  =  ^4"MS 
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In  practice  vever,  this  state  can  never  be  attained  experimentally, 
since  at  z(  internal  field  the  magnetization  has  dropped  from  its  satura¬ 
tion  value  to  its  remanence  value,  which  is  usually  substantially  lower. 
Values  of  magnetization  somewhat  higher  than  remanence,  but  still  below 
saturation,  can  be  obtained  with  magnetic  fields  small  enough  not  to  con¬ 
tribute  significantly  to  the  permeability  themselves. 

In  the  partially  magnetized  state  the  imaginary  parts  of  perme¬ 
ability  components  become  very  large  if  u>  2  1.  producing  what  is  called 
"low  field  loss.  "  This  effect  has  been  explained  by  Polder  and  Smit.4  Be¬ 
cause  of  it,  ferrite  devices  that  operate  in  the  partially  magnetized  state 
always  operate  with  <  1  • 

A  theory  of  Rado5  states  that  if  u^/u  «  1 ,  p'  and  p'  are  both 
equal  to  1.  0  for  all  unsaturated  states  and  k'  =  y4nM/u  where  4irM  is  the 
microscopic  magnetization  averaged  over  the  domains.  Although  it  has 
been  generally  accepted  that  Rado's  theory  for  k'  is  a  good  approximation 
even  as  u>M/u>  gets  fairly  close  to  1.0,  the  theory  for  p'  and  p^,  is  valid 
only  in  the  limit  -*  0.  Rado's  results  are  derived  by  averaging  the 
torque  equation 

4ttM  =  ■y4rrM  X  H  (1) 

over  the  sample.  By  making  several  ptions,  including  /u  «  1. 0, 

on  the  relative  size  of  various  components  of  the  magnetization  4itM  and 
the  magnetic  field  H,  he  reduced  the  torque  equation  at  a  given  point  in  the 
material  to 


\4nrM 

4irm  =  i - -  HXu 

w 


(2) 


Here  u  is  a  unit  vector  in  the  direction  of  the  static  magnetization  in  the 
domain  in  question,  4Trm  is  the  rf  magnetization,  and  H  is  the  rf  magnetic 


.  jwm*  *■ 


field.  The  field  H  includes  a  spatially  constant  part  consisting  of  the  applied 
rf  drive  field  less  the  demagnetizing  fields  from  the  sample  boundary  and  a 
spatially  fluctuating  component  due  to  demagnetizing  fields  from  poles  on  the 
domain  and  grain  boundaries.  Although  the  spatial  average  over  the  sample 
of  this  fluctuating  field  is  zero,  it  is  the  source  of  the  deviations  of  p.'  and 
from  1. 0  as  well  as  the  small  deviations  in  k1  from  y4irM/u  which  we  have 
observed  in  some  materials.  Unfortunately  Rado  did  not  take  into  account 
the  effect  of  this  spatially  fluctuating  field. 

When  “jyj/w  «  1  *  these  spatially  fluctuating  fields  become  small 
and  hence  Rado's  results  become  more  accurate.  The  condition  «  1 

was  not  required,  however,  at  the  point  where  he  invoked  it,  namely  in  sim¬ 
plifying  the  torque  equation. 

In  order  to  derive  a  bulk  susceptibility  from  the  point  relationship. 
Eq.  (2),  one  must  take  a  spatial  average  of  Eq.  (2)  over  the  whole  sample 

y4trM 

<  4"m  >av  =  i - —  <  K  x  u  >ay  (3) 


By  writing 

K  =  <E>  +  6K 

av 

u  3  <u>av  +  <5u 

where  <5E  and  <5u  are  the  spatially  fluctuating  components  of  K  and  u,  one 
obtains  with  a  small  amount  of  manipulating 

<E  x  u  >  =  <E>  x  <u  >  +  <  <5E  x  <5u  >  (4) 

av  av  av  av 

The  error  in  Rado's  theory  is  his  assertion  that  <5E  and  <5u  are  uncorrelated 

and  hence  <  6E  x  6u  >av  =  0.  Using  this  assumption,  and  the  fact  that 

<u  >  =  4ttM/4ttM  ,  Rado's  results  that  n'  =  n'  -  1  ■  0  and  *'  -  4irM/u  are 

av  s  z 

easily  obtained. 
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It  will  be  shown  below  that  the  assumption  that  6u  and  <5E  are 
uncorrelated  is  incorrect.  The  demagnetizing  fields  from  poles  on  domain 
and  grain  boundaries  far  from  the  point  in  question  average  to  zero  and 
hence  do  not  contribute  to  <5K.  The  fields  from  poles  close  to  the  point  in 
question  but  not  on  its  own  domain  or  grain  boundary  contribute  to  <5E  and 
may  be  reasonably  uncorrelated  with  6u.  It  is  the  field  due  poles  on  the 
boundary  of  the  domain  or  grain  in  question,  which  is  a  substantial  part  of 
6E,  that  is  strongly  correlated  to  <5u.  These  poles  are  produced  by  the  rf 
magnetization  of  the  domain  which  is  precessing  about  a  direction  deter¬ 
mined  in  large  part  by  u  (as  well  as  by  the  domain  shape).  If  we  change 
the  direction  of  the  u,  the  precessional  path  of  the  rf  magnetization  will 
change  with  it  and  hence  demagnetizing  field  of  that  domain  or  grain  which 
is  contributing  to  6E  will  char.ge.  Thus  6E  is  correlated  with  6u. 

If  we  look  at  an  oversimplified  picture  of  a  single  domain  we  can 

see  why  this  neglect  of  correlation  produces  errors  primarily  in  p'  and  p' 

z 

rather  than  in  k'  We  will  consider  the  domain  to  be  an  ellipsoid  of  revolu¬ 
tion  about  the  direction  of  u.  We  will  also  consider  only  the  contribution  to 
<5E  is  from  the  pole?.:  on  its  own  boundary.  All  other  demagnetizing  fields, 
including  those  from  adjacent  domains  and  grains,  are  assumed  to  sum  to 
zero.  The  precessing  magne+ization  driven  by  the  applied  field  <E>iv  will 
have  two  orthogonal  components  in  the  plane  perpendicular  to  <5u.  These 

are  in  the  directions  of  6u  x  <E>  and  5u  x  (6u  x  <E>  ).  The  rf 

av  av 

magnetization  will  produce  a  6E  in  the  same  direction  as  the  rf  magnetiza¬ 
tion  and  related  to  it  by  the  transverse  demagnetizing  factor  of  the  ellipsoid. 

In  computing  the  average  <<5u  x  5h>  we  will  thus  have  terms  proportional 

sv 

to 


<5u  x  (<5u  x  <E>av)  and  x  (0U  x  (<5u  x  <E>av))  (5) 

We  now  assume  further  that,  for  each  domain  like  the  one  under  considera¬ 
tion,  there  is  one  of  the  same  shape  with  opposite  6u.  This  assumption  is 
probably  not  exactly  true  since  only  the  average  value  of  6u  is  zero.  It  is 
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not  necessary  there  be  an  equal  number  of  domains  with  tdu  of  each  size 
and  shape.  However,  within  the  framework  of  this  assumption  we  see  the 
following.  If  we  change  the  sign  of  6u  in  the  term  above  containing  the  series 
of  cross  products  with  three  du's,  the  sign  of  the  term  changes  and  hence 
its  average  over  the  sample  is  zero.  However,  the  average  <6u  x  dK> 

8V 

still  contains  the  remaining  term  <  <5u  x  (fiu  x  <K>av)>av  which  is  not  zero. 
Applying  a  vector  identity  to  this  term  gives 

<5u  x  (<5u  x  <H>  )  =  (6u  •  <E>  )  <5u  -  ($u  •  3u)  <E>  •  (6) 

81 V  8V 

We  can  define  the  x  direction  so  that  <K>  is  in  the  x-z  plane.  Taking  the 

3V 

average  of  the  first  term  on  the  right  side  of  Eq.  (6)  above  over  all  domains 
gives 


< fiu  •  <E>  )  <5u>  =  <du2>  x  +  <6u2>  <h>  z 

av  av  x  av  z  av  av 


(7) 


where  x  and  “z  are  unit  vectors  in  the  x  and  z  directions.  The  average 
of  of  the  second  term  in  Eq.  (6)  is  just 

<(aE  •  fiE)  <K>  >  -  <  l«ul2>  (<h  >  x  +  <h  >  z)  .  (8) 

'  av  av  1  1  av  '  x  av  z  av 

Combining  these  averages  with  Eqs.  (3)  and  (4)  gives  contributions  to  the 
rf  magnetization  of 


<4nm 


> 

x  av 


<h 


> 

x  av 


<4irm 


> 

z  av 


(n;-i)  <h 


> 

z  av 


(9) 


where 

(n1-  1 )  oc  <  5  U2  >av  -  <  j  <5u  1 2  >av 

(K^-l)  oc  <  6u“>av-  <  [au|2>av  (10) 
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The  constant  of  proportionality  depends  on  the  appropriate  sus¬ 
ceptibility  term  of  an  isolated  magnetic  moment  <5u  times  its  average 
transverse  demagnetizing  factor.  No  additional  contribution  to  k'  is  made 
in  this  approximation.  It  must  be  emphasized  again  that  very  rough  ap¬ 
proximations  are  made  in  this  paragraph  and  that  *'  may  be  expected  to 
deviate  from  Rado's  theory,  but  the  deviation  is  less  than  that  for  p'  and  p^. 

It  should  also  be  emphasized  that  there  is  no  a  priori  reason  for 
the  components  of  the  permeability  to  be  unique  functions  of  the  dc  magnet¬ 
ization.  States  having  the  same  magnetization  but  on  opposite  sides  of  B-H 
hysteresis  curve  have  different  domain  structures,  and  hence  the  spatial 
average  <  6u  x  6E>  should  not  be  expected  to  be  the  same  for  them.  The 
differences,  however,  are  second-order  effects  as  are  the  deviations  of 
from  Rado's  theory. 

Other  authors^*’ ^  have  theoretically  considered  the  permeability 
of  partially  magnetized  ferrites  from  a  superficially  different  point  of  view. 
They  start  with  the  tensor  permeability  of  a  randomly  oriented  domain  and 
average  over-all  orientations  arriving  at  the  same  result  as  Rado.  This 
is  to  be  expected  since  their  tensor  permeability  is  derived  from  the  same 
torque  equation  Rado  used,  and  they  too  neglected  the  spatially  fluctuating 
rf  magnetic  field  within  the  material. 

It  has  not  yet  been  possible  to  evaluate  the  average  iu  x  6h  > 
in  terms  of  a  physical  model  of  the  domain  structure  and  hence  give  explicit 
formulae  for  p'  and  p^  In  this  report  we  therefore  present  graphs  of  our 
measurements,  along  with  phenomenological  formulae  for  use  in  device 
design. 


A  separate  theoretical  analysis  of  p'  and  p'  in  the  totally  demag- 

2  g 

netized  state  (where  they  are  equal)  was  made  by  Sandy  and  Green  in  which 

a  somewhat  different  approach  was  taken.  The  spatially  fluctuating  rf  fields 

4 

were  treated  in  the  manner  used  by  Polder  and  Smit  to  explain  low-field  loss 


when  u^j/u  >  I  -  0. .  They  showed  that  these  demagnetizing  fields  cause  the 
domains  to  undergo  ferromagnetic  resonance  at  frequencies  from  zero  to 
While  we  are  concerned  with  the  region  u^/u  <1.0,  it  is  the  high  frequency 
tails  of  these  resonances  that  make  p’  and  p'  differ  from  1. 0.  Using  this 
approach,  the  domains  were  assumed  to  have  an  unknown  distribution  function 
D(Up)  of  resonant  frequencies  between  zero  and  and  that  each  of  these 
resonances  was  Lorentzian  in  shape.  By  integrating  various  guesses  of 
D(Upj)  times  the  Lorentzian  lineshape,  it  was  found  that  the  frequency  de¬ 
pendence  of  the  permeability  above  was  not  strongly  dependent  on  the 
choice  of  distribution  function,  and  reasonable  functions  were  found  that 
fit  experimental  data  on  the  frequency  dependence  of  the  initial  permeability 
of  several  ferrites.  However,  no  single  distribution  function  fit  all  ferrites, 
indicating  the  presence  of  significantly  different  types  of  domain  and  grain 
structures. 

In  considering  interactions  between  magnetic  moments  and  mag¬ 
netic  fields,  it  is  generally  valid  to  scale  all  magnetizations  and  Helds  to  u/y. 
Thus  the  question  immediately  arises  whether  '  and  p’  depend  on  4ttM, 

4irM  ,  and  u  independently  or  only  on  the  two  ratios  y4ttM/w  and  y4trM  /u. 

3  3 

This  suggested  dependence  was  confirmed  experimentally  for  \4irM  /u  in 

8  ® 
the  work  cited  earlier.  Further  confirmation  for  y4nM/u  will  be  given 

later  in  this  report. 

We  have  also  studied,  experimentally,  the  imaginary  parts  of  the 
permeability  tensor,  p",  p"  ,  and  k".  In  addition  to  all  factors  important 
for  the  real  parts,  the  imaginary  parts  depend  on  the  intrinsic  loss  mech- 
ansim  of  the  material  and  hence  on  chemical  composition,  impurities,  and 
temperature.  They  also  depend  strongly  on  the  degree  of  internal  strain 
in  the  crystal  structure.  This  shows  up  as  large  surface  losses  in  machined 
samples  that  have  not  been  annealed  to  relieve  strains  caused  by  machining. 

Our  measurements  to  date  have  been  primarily  on  two  commonly 
used  commercial  ferrites  made  by  Trans  Tech:  G-113,  an  yttrium  iron 
garnet  with  4irMg  =  17  50  gauss  at  room  temperature;  TT1-390,  a  magnesium 

1 1 


manganese  spinel  with  4ttM  =  2150  gauss  at  room  temperature;  and  on 
MMF-10,  another  commercial  magnesium  manganese  spinel,  made  by 
Airtron,  having  a  4irMg  =  2400  gauss  at  room  temperature.  This  latter 
material  was  included  because  it  contains  no  second  phase  while  a  sub¬ 
stantial  amount  is  present  in  TT1-390,  and  it  was  desired  to  determine 
if  this  was  a  major  factor  in  any  of  the  observed  properties.  Measure¬ 
ments  were  also  made  on  YIG  with  1  per  cent  Dy  in  order  to  improve  the 
sensitivity  of  the  loss  measurement  Most  of  our  measurements  were 
made  at  5.  5  GHz  where  was  close  to  1. 0  for  all  three  materials  at 

room  temperature.  By  heating  the  samples,  u^/u)  could  be  continuously 
reduced  as  the  Curie  temperature  was  approached.  It  is  assumed  that  the 
prime  effect  of  heating  on  p' ,  p^,,  and  k'  arises  from  the  changes  in  w^/u. 
Other  less  important  changes  may  occur  due  to  changes  in  factors  that 
affect  the  domain  structure.  Some  evidence,  including  measurements  on 
other  materials  at  9.20  GHz  and  room  temperature,  will  be  presented 
showing  that  these  effects  are  small.  The  direct  effect  of  temperature  on 
intrinsic  loss  parameters  on  the  other  hand,  is  not  small.  However,  the 
dependence  of  p",  p^),  and  k"  onu^/u  is  so  drastic  when  *s  only  slightly 

<  0.1  that  the  temperature  dependence  of  the  intrinsic  loss  parameters  causes 
a  relatively  small  change  in  them.  To  eliminate  this  effect,  measurements  are 
also  being  made  on  p"  in  the  demagnetized  state  at  room  temperature  with 
frequency,  <j,  varied  and  held  constant.  In  that  case  the  losses  will 
contain  the  frequency  dependence  rather  than  the  temperature  dependence 
of  the  relaxation  rate 

B.  Measurements  of  p'  and  * ' 

To  measure  p'  and  k' ,  we  have  used  a  cylindrical  TM^jq  cavity 
similar  to  that  used  by  LeCraw  and  Spencer  2  The  sample,  a  0.  075  in. 
diam  rod,  is  mounted  along  the  axis  of  the  cavity.  The  cavity  supports 
two  degenerate  modes  which  are  split  into  left  and  right  circularly  polarized 
modes  whenever  k'  *  0.  The  frequencies  of  these  two  modes  are  related  to 
p.'  +  k'  and  p’  -  k  1 ,  respectively.  The  cavity  has  four  symmetrically  placed 
irises.  By  feeding  two  orthogonal  ports  ±  90°  out  of  phase,  either  mode 
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could  be  excited  separately.  To  improve  the  degeneracy  of  the  cavity, 
a  ring  of  eight  tuning  screws  was  placed  in  one  end  of  the  empty  cavity, 
as  seen  in  Fig.  1 . 

The  cavity  was  mounted  with  the  sample  rod  between  the  pole 
faces  of  a  variable  gap  electromagnet  as  shown  in  Fig.  2.  By  closing  the 
gap  until  the  pole  faces  nearly  contacted  the  rod,  the  demagnetizing  poles 
on  the  ends  of  the  rod  were  nearly  eliminated  and  uniformity  of  the  mag¬ 
netization  is  improved.  The  rods  used  for  the  measurements  at  5.  5  GHz 
were  3  in.  long;  those  used  at  9.  2  GHz  were  2  in.  long.  The  5.  5  GHz 
measurements  were  made  as  a  function  of  temperature  in  order  to  vary 

4irM  and  hence  u,  ,/u.  - - -  - 

3  M 

In  order  to  know  the  magnetization  of  the  rod,  as  the  field  of  the 
electromagnet  was  swept,  a  2000  turn  pickup  coil  was  placed  around  the 
rod  just  outside  the  cavity.  This  coil  was  wired  in  series  with  an  adjacent 
2000  turn  coil  not  containing  the  sample,  wound  in  the  opposite  direction 
in  order  to  buck  the  flux  of  the  magnet's  H  field.  The  two  coils  were  con¬ 
nected  to  an  integrator  whose  output  was  proportional  to  the  magnetization 
of  the  rod.  By  mounting  the  rod  and  coil  in  the  magnet  gap  without  the 
cavity  and  moving  the  coil  along  the  rod,  we  found  the  magnetization  over 
the  parts  of  the  rod  normally  in  the  cavity  or  the  coil  was  constant  within 
a  few  percent  for  the  3  in.  rods  and  somewhat  worse  for  the  2  in.  rods. 

The  cavity  frequency  was  obtained  as  a  function  of  magnetization 
as  the  magnetization  was  cycled  about  a  major  hysteresis  loop,  by  making 
the  cavity  part  of  a  pound  discriminator  locking  the  klystron  to  the  cavity. 

The  correction  voltage,  fed  to  the  klystron  reflector,  was  calibrated  to  the 
klystron  frequency  and  plotted  vs  the  output  of  the  integrator  on  an  x-y  re¬ 
corder.  This  was  done  for  both  circular  polarization  modes  of  the  cavity.  From 
these  two  frequencies,  the  values  of  and  k'  could  be  computed  as  a  function 
of  magnetization.  This  could  be  done  either  by  perturbation  theory  as  de¬ 
scribed  by  LeCraw  and  Spencer^  or  by  an  exact  solution  given  by  Bussey 
o 

and  Steinert.  The  two  methods  are  nearly  equivalent  at  5.  5  GHz  where 
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Fig.  1  Cylindrical  TM 


most  of  our  measurements  were  made.  However,  at  9.2  GHz  the  0.  075  in. 
rods  are  too  large  for  perturbation  theory  to  be  accurate.  We  have  thus 
written  our  data  reduction  computer  programs  to  use  the  exact  solution. 

Measurements  of  p.'  and  *'  on  MMF-10  were  obtained  as  a  by¬ 
product  of  measurements  on  n;'  and  k"  using  a  different  technique  described 
in  Sec.  D. 


The  results  of  the  measurement  of  k1  vs  magnetization  on  G-113 
were  reasonably  consistent  with  Rado's  theory.  The  measured  values  of 
k'  were  generally  a  bit  higher  than  theoretically  expected  with  the  greatest 
deviations  (less  than  0.  05)  occurring  at  y4ttM/u  «  0.  4.  This  can  be  seen 
in  Fig.  3  for  k  at  room  temperature  at  5.  5  GHz  and  9-  2  GHz.  The  devia¬ 
tions  were  similar  at  all  temperatures  at  5.  5  GHz,  and  no  real  dependence 
of  u^/u  could  be  found.  Of  significance,  however,  was  the  fact  that  for 
G-113  there  was  no  dependence  of  on  hysteresis.  Although  not  equal 
exactly  to  y4irM/u;,  it  was  a  unique  function  of  it. 


The  measured  dependence  of  and  y4nM/u  and  y4rrMg/u  for 
G-113  at  5.  5  GHz  is  shown  in  Fig.  4.  4rrM  is  varied  by  changing  tempera- 

S 

ture.  Again  there  is  no  dependence  on  hysteresis  history. 


The  solid  curve  in  this  figure  is  the  function 


|i>  =  ^  +  (1  - 


tanh  (1. 25  (M/Mgr) 
tanh  (1.25) 


(ID 


where  is  the  value  of  when  4ttM  =  0.  It  can  be  seen  that  this  function 
fits  the  data  reasonably  well  over  the  entire  range  of  values  of  y4irM/u  and 
Y4irMg/ u.  While  this  functional  relationship  has  no  physical  significance, 
it  is  intended  to  be  useful  in  optimizing  microwave  ferrite  devices  with  the 
aid  of  computer  programs.  The  deviation  of  the  data  from  the  solid  curve 
at  the  highest  temperature  (lowest  w^/u)  is  at  least  in  part  due  to  tempera¬ 
ture  gradients  in  the  cavity.  The  point  on  the  outside  of  the  cavity  at  which 
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our  thermocouple  was  attached  may  have  been  a  few  degrees  cooler  than 
the  sample  inside  the  cavity.  As  a  result  the  end  point  of  the  solid  curves 
is  at  too  high  a  value  of  y4ttM/w  for  the  actual  u^/u  of  the  sample. 

In  making  use  of  the  values  of  u^/u  quoted  for  each  curve  in  this 
and  subsequent  figures,  it  should  be  remembered  that  due  to  uncertainties 
in  the  materials'  g  factor,  and  hence  in  y,  and  due  to  errors  in  the  measure 
ment  of  4irMg,  there  will  be  errors  of  a  few  percent  in  u^/u. 

In  an  earlier  report  we  fit  this  data  with  a  simpler  function: 

I*'  *  '  d-^)(M/Ma)3/2  .  (12) 

The  function  in  Eq.  (11)  does  not  give  a  better  fit  than  the  one  in  Eq.  (12). 
However,  it  does  give  a  reasonable  fit  to  the  data  for  TT1-390  and  MMF-10 
as  well  as  for  G-113  which  the  simpler  function  failed  to  do. 

At  the  highest  values  of  4trM  shown  for  each  curve  in  this  figure, 
the  internal  field  is  about  50 Oe.  As  a  result,  there  is  a  small  contribution 
to  p'  from  the  low-field  tail  of  the  ferromagnetic  resonance  curve.  This 
reduces  p'  below  that  produced  by  its  domain  structure  alone  and  prevents  p' 
from  approaching  1 . 0  as  4ttM  approaches  4irMg. 

The  measurement  made  on  TT1-390,  a  magnesium -manganese 
spinel,  showed  quite  different  results.  At  values  of  close  to  1. 0 

both  p'  and  k'  showed  significant  dependence  on  hysteresis  history  rather 
than  depending  on  magnetization  aione.  This  can  be  seen  in  Figs.  5  and  6. 
This  dependence  decreases  as  decreases  but  a  residual  amount  was 

detectable  at  all  temperatures .  It  should  be  noted  that  on  major  hysteresis 
loops  the  two  crossings  at  4ttM  -  0  yield  the  same,  value  of  p’,  while  the 
values  of  k’  are  equal  in  magnitude  but  opposite  in  sign.  This  results  from 
a  symmetry  of  the  permeabijity  tensor  which  requires  that  p  be  an  even 
function  of  magnetization  and  its  hysteresis  while  k  is  an  odd  function. 
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(a) 


(b) 

Fig.  5  k'  vs  \4ttM/ u  on  TT1-390 
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In  Fig.  7  we  show  both  the  experimental  p'  and  that  of  Eq.  (11)  as 
a  function  of  -y4TrM/u  for  several  values  of  obtained  by  heating  the 

sample.  Only  one  branch  of  each  hysteresis  curve  is  shown  in  this  figure 
for  clarity.  The  other  branch  is  the  mirror  image  of  it  as  can  be  seen 
from  Fig.  6.  In  this  and  all  other  figures  where  the  full  hysteresis  curve 
is  not  shown,  the  direction  of  the  branch  shown  is  from  positive  magnetiza¬ 
tion  to  negative  magnetizations. 

It  had  been  suggested  that  the  reason  hysteresis  occurred  in  p' 
and  k'  vs  y4wM/u  for  TT1-390  and  not  for  G-113  was  that  TT1-390  had 
a  substantial  amount  of  second  phase.  Measurements  have  thus  also  been 
made  on  a  single-phase  magnesium  manganese  spinel,  MMF-10,  manu¬ 
factured  by  Airtron.  The  results  for  p'  are  shown  in  Fig.  8  and  are  not 
qualitatively  different  from  TT1-390. 

As  a  confirmation  that  our  results  for  varying  by  heating  the 

sample  at  constant  frequency  are  valid  and  also  as  a  confirmation  that  fre¬ 
quency  is  not  an  independent  variable  in  addition  to  y4irM/u  and  y4irM  /u, 

s 

we  show  in  Fig.  9  p'  vs  y4irM/u  at  9.2  GHz  for  G-113,  R151  (a  Raytheon 
ferrite  very  similar  to  TT1-390)  and  MZMF,  a  Raytheon  magnesium,  zinc, 
manganese  ferrite  with  a  4irMg  -  3000  gauss  at  room  temperature.  The 
shapes  of  these  curves,  including  their  deviations  from  the  solid  curves 
showing  Eq.  (11),  are  similar  to  the  comparable  curves  at  5.  5  GHz.  The 
only  significant  difference  occurred  for  the  MZMF;  in  this  material  p'  is 
noticeably  larger  than  the  functional  form  which  fits  the  other  materials 
quite  well.  In  addition  its  hysteresis  in  p'  vs  y4itM/u>  has  the  opposite 
sign  to  that  of  TT1-390  and  MMF-10. 

C.  Measurements  of  p' 

- cz 

Measurements  of  pz  were  made  on  samples  placed  in  the  center 
of  a  rectangular  TE1Q2  reflection  cavity.  Originally  we  used  rod-shaped 
samples  with  the  rf  and  dc  magnetic  fields  parallel  to  each  other  and  per¬ 
pendicular  to  the  rod.  This  technique  proved  unsatisfactory  since  a  rod 
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vs  v4uM/u>  on  MMF-10  at  5.5  GHz 


-113 


tends  to  magnetize  along  its  axis  if  there  is  the  slightest  misalignment  of 
the  dc  field.  The  measurements  reported  here  were  made  on  spherical 
samples.  Although  this  greatly  reduced  the  cavity- filling  factor,  and  hence 
the  sensitivity,  it  eliminated  the  serious  errors  due  to  misalignment. 

Since  moderate  internal  dc  fields  parallel  to  the  rf  field  do  not 

cause  the  sample  to  approach  ferromagnetic  resonance,  it  was  possible 

to  measure  (jl  at  magnetization  substantially  closer  to  4wM  than  we  could 
z  s 

when  measuring  |i  and  k. 

Measurements  were  made  point  by  point,  varying  the  dc  magnetic 
field  from  values  large  enough  to  saturate  the  sample  in  one  direction 
through  zero  to  values  large  enough  to  saturate  it  in  the  negative  direction. 
For  magnetizations  below  remanence  4"M  =  3  H^.  For  high  magnetizations, 
where  the  internal  magnetic  field  is  significant,  values  of  4ttM  vs  ob¬ 
tained  from  a  separate  measurement  of  magnetization  were  used.  Apparent 
values  of  n!  and  p "  were  obtained  from  the  shift  in  cavity  resonant  fre¬ 
quency  and  reflection  coefficient,  using  the  value  obtained  with  a  9  5  kOe 

field  as  a  reference.  Internal  values  of  u  were  then  obtained  from  the 

z 

apparent  p.  using  the  formula 
z 

=  apparent  _  ;  -  -  <i3a) 

*  3  ^z  apparent 

The  results  for  MMF-10  and  YIG  with  1  percent  Dy  are  shown  in 
Figs.  10  and  11.  The  solid  curves  are  the  function 

(1-lM/M  >5/2) 

<13b> 

£>  4  O 

Again  it  should  be  emphasized  that  this  functional  relationship  is  given 
because  it  fits  the  data  reasonably  well  and  hence  would  be  useful  in  computer 
aided  design  programs  for  ferrite  devices.  No  real  physical  basis  for  this 
relationship  is  implied. 
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2 


vs  yUM/u,  on  YIG  1  percent  Dy  at  5.  5  GHz 


on  MMF-10  at  5.  5  GHz 


In  this  report  we  have  discussed  fitting  our  data  to  three  different 
functional  forms 

Nl 

Fj  =  p'Q  +  (l-^HM/Mg)  l) 


(14) 

Each  of  these  functions  is  identically  p^  at  4ttM  =  0  and  1. 0  at  4rrM  =  4nMg 
It  is  worthwhile  to  mention  how  the  functions  F2  and  differ  in  shape  for 
0 < 4rrM <  4rrMg  from  the  more  simple  function  F^.  Basically  F^  and  F^ 
with  constant  N2  and  N3  are  the  same  as  Fj  with  a  variable  Nj.  F2  is  the 
same  as  Fj  with  Nj  =  N2  at  4irM  =  0  and  an  Nj  that  decreases  as  4ttM  ap¬ 
proaches  4rrMg.  F^  on  the  other  hand  is  the  same  as  Fj  with  the  Nj  =  N3 
when  4irM  =  4TrMg  and  an  Nj  that  decreases  as  4ttM  approaches  zero.  As 
each  N  is  increased,  all  three  functions  become  flatter  at  4irM  =  0  and 

increase  more  sharply  as  4ttM  approaches  4irM  .  It  is  found,  however,  in 

s 

fitting  our  data  to  F,  that  if  we  choose  Nj  so  that  the  function  fits  our  data 
for  small  4ttM,  then  the  function  is  too  large  for  p'  and  too  small  for  p'  at 

4 

large  4nM.  Thu^  we  use  the  functional  for  F2  ofr  p'  and  F^  for  p^,. 

The  question  of  whether  there  is  hysteresis  in  p^,  i.  e  ,  whether  p 
is  a  unique  function  of  magnetization  or  depends  on  which  side  of  the  B-H 
hysteresis  loop  the  sample  is  on,  is  not  as  clear  cut  as  it  was  for  p'  and  k’ 
because  our  measurement  technique  is  less  sensitive  than  that  for  p'  and  «' 
In  addition,  being  a  point -by- point  measurement  rather  than  a  continuous 
measurement,  it  is  subject  to  small  drift  errors.  No  hysteresis  was  de¬ 
tectable  in  \±'z  for  the  YIG.  A  small  amount  of  hysteresis  is  apparent  is 
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apparent  for  MMF-10  at  low  values  of  4irM  and  values  of  u^/u  >  0*  9.  This 

values  of  Uj^/u  is  greater  than  the  maximum  value  for  YIGT  however.  The 

hysteresis  in  n'  above  this  value  shows  up  as  smaller  values  of  ji'  at  nega- 
z  z 

tive  magnetizations  than  at  the  corresponding  positive  magnetization. 

D.  Measurements  of  (*"  and 

The  imaginary  parts  of  p.  and  k  were  measured  with  the  same  cavity 
as  the  real  parts.  However,  since  it  is  not  possible  to  obtain  a  continuous 
plot  of  Q  vs  magnetization  as  the  rod  is  cycled  about  the  B-H  hysteresis 
loop,  point- by- point  measurements  must  be  made.  The  cavity  Q  and  fre¬ 
quency  were  recorded  as  a  function  of  dc  magnetic  field  since  magnetization 
could  not  be  measured  directly  point  by  point.  A  separate  plot  of  4trM  vs  H 
was  then  made.  Using  these  results,  one  could  calculate  =  p"-«"  as 
well  as  p^  =  p'-x'  as  a  function  of  either  field  or  magnetization  using  per¬ 
turbation  theory.  This  was  done  only  for  our  preliminary  measurements 
on  G-113  used  to  test  the  system.  For  subsequent  measurements  the  exact 
“solutions  for  a  rod  in  a  cavity  were  used.  This  solution  required  a  knowledge 
of  cavity  Q  and  resonant  frequency  for  both  senses  of  circular  polarization 
at  each  magnetization.  This  did  not  mean  that  measurements  had  to  be 
made  with  both  senses  of  circular  polarization  since  the  values  for  one 
sense  are  the  same  as  those  of  the  opposite  sense  when  the  rod  is  oppositely 
magnetized  both  in  magnitude  and  hysteresis  history.  Since  measurements 
were  made  at  points  with  uniformly  spaced  magnetic  fields,  data  was  not 
directly  available  at  points  with  equal  and  opposite  magnetization.  There¬ 
fore,  we  drew  smooth  curves  through  the  raw  data  and  used  the  interpolated 
values  from  the  smoothed  curves  rather  than  the  raw  data  to  do  the  sub¬ 
sequent  calculations  The  data  reduction  using  the  exact  solution  yields 
(i1,  h",k',  k"  directly  rather  than  and  p“  which  are  obtained  from  per¬ 
turbation  theory. 

One  of  the  most  striking  effects  which  we  observed  in  our  early 
loss  measurements  was  a  large  surface  loss  caused  by  strains  induced  by 
grinding  the  sample  to  shape  Although  we  were  aware  of  the  existence 
of  surface  loss,  our  first  preliminary  loss  measurements  were  made  on 
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a  rod  of  G- 113  which  had  not  yet  been  annealed  to  remove  the  surface  strain 
and  loss.  We  had  expected  the  effect  would  be  to  increase  n”  =  p"  -  <c"  by  a 
constant  scale  factor  related  to  ratio  of  surface  loss  to  bulk  loss  and  the 
fractional  volume  of  the  strained  surface  layer.  We  were  therefore  sur¬ 
prised  to  find  the  results  for  p+'  vs  shown  in  Fig.  12.  In  this  figure 
we  also  show  the  magnetic  hysteresis  loop  5  vs  H.  It  can  be  seen  that 

vs  H  does  not  close  until  much  larger  magnetic  fields  than  those  neces¬ 
sary  to  close  the  B  vs  K  hysteresis  loop  This  phenomenon  results  from 
the  fact  that  the  strained  surface  layer,  whose  loss  is  actually  the  dominant 
loss  observed  in  spite  of  its  small  volume,  also  has  a  very  high  coercive 
force.  The  observed  thus  has  the  very  different  hysteresis  properties 
of  the  surface  layer  rather  than  the  average  measured  hysteresis  properties 
which  correspond  primarily  to  the  bulk  material.  The  results  obtained  after 
this  sample  was  annealed  are  shown  in  Fig.  13.  The  magnitude  of  is 
greatly  reduced,  the  two  peaks  occur  at  the  coercive  fields,  and  the  loops 
close  at  the  same  field  as  the  B-H  loop  shown  in  Fig.  12  closed. 


The  scatter  shown  at  negative  fields  is  due  to  the  rf  demagnetizing 
fields  on  the  sample  surface.  These  demagnetizing  fields  necessitate  a 
correction  of  the  following  form  to  convert  the  apparent  circular  suscepti¬ 
bility  (p  -1)  which  is  proportional  to  the  change  in  complex  frequency 
(f  +  if/2Q)  of  the  cavity  to  true  susceptibility  (p  -1) 


u  -  1  * 


2(u 

+  app 


-1 


'**+  app"”^ 


(15) 


Writing  p+  and  p+  p  in  complex  form ,  separating  real  and  imaginary 
parts  and  noting  that  p  +app  <<:  ^4- app  gives 


“i-  1 


2Kapp-1} 

2*  '^app-1' 

^  +  app 
app- 1)1 


p"  (1+p'  )2 

r  +  app  v 


il6) 


(16) 
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Since  the  system  sensitivity  is  inversely  proportional  to  the  minimal  de¬ 
tectable  change  in  l/Q  and  this  is  proportional  to  ,  the  scatter  in  p" 
is  proportional  to  (l+p'+)^.  Since  1  +  p'+  varies  almost  from  1. 0  to  3.  0, 
the  change  in  sensitivity  is  almost  an  order  of  magnitude  as  the  direction 
of  magnetization  changes.  Although  the  sensitivity  was  improved  for  sub¬ 
sequent  measurements,  it  was  felt  that  differences  between  p"  and  x"  and 
their  dependence  on  y4ttM/cj  could  be  more  accurately  observed  in  an  in¬ 
trinsically  lossier  material.  We  therefore  used  a  Y1G  doped  with  1  per¬ 
cent  Dy  rather  than  G-113  in  subsequent. 

The  results  of  p"  and  x!l  for  the  YIG  with  1  percent  Dy  are  shown 
in  Figs.  14  and  15.  This  and  all  other  materials  referred  to  in  the  rest 
of  this  and  the  two  following  sections  have  been  annealed.  The  variation 
in  was  obtained  by  heating  the  sample.  Although  this  also  changes 

the  relaxation  rate,  the  shapes  of  p"  and  x"  with  (j^/u  at  all  values  of 
y4ttM/u  is  primarily  due  to  the  drop  in  u^/u.but  does  depend  to  a  notice¬ 
able  degree  on  the  change  in  the  relaxation  rate  as  UM/W  becomes  sub¬ 
stantially  less  than  1.0.  It  is  interesting  to  note  that  as  gj^/u;  decreases, 
the  peak  in  x"  at  4irM  =  0  disappears  and  then  becomes  a  minimum.  It  is 
believed  that  this  is  caused  by  the  following.  The  observed  losses  are  due 
to  the  high-frequency  tails  of  ferromagnetic  resonance  curves.  Near 
remanence  these  curves  are  centered  about  a  resonant  frequency  of  zero 
since  the  net  internal  dc  magnetic  field  including  anisotropy  and  demag¬ 
netizing  fields  is  close  to  zero.  In  the  demagnetized  state.  Polder  and 

4 

Smit  have  shown  that  these  resonances  occur  at  all  frequencies  from 
zero  to  When  u  is  close  to  u>^,  the  resonances  close  to  predom¬ 

inate  and  make  p"  larger  in  the  demagnetized  state  than  near  remanence. 
When  u  »  u^,  all  of  the  resonance  tails  contribute  equally  and  there  is 
no  tendency  for  p"  to  peak  at  4ttM  =  0.  However,  in  the  demagnetized  state, 
one  third  of  the  domains  are  parallel  to  the  rf  drive  field  and  do  not  con¬ 
tribute  to  the  ferromagnetic  resonance  tail.  Thus  p"  can  have  a  minimum 
at  4irM  =  0  for  sufficiently  low  values  of  u^j/gj.  This  of  course  assumes 
that  p"  «  p"  for  4ttM  =  4ttM  and  an  internal  field  equal  zero.  There  is 

Z  S 

however  a  very  small  but  finite  p'^  due  to  subthreshold  parallel  pumping. 
This  is  discussed  in  the  next  section. 
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percent  Dy  at  5.  5  GHz 


Fig.  15  k"  vs  y4ttM/u>  ior  YIG  1  percent  Dy  at  5.5  GHz 


It  should  also  be  noticed  in  Figs.  14  and  15  that  for  all  values  of 
obtainable  with  YIG,  k"  «  p".  Since  k  is  an  odd  function  of  4rrM, 
k "  =  0  when  4ttM  =  0.  It  is  also  zero  theoretically  in  the  fully  saturated 
state  if  H^  =  0.  However,  a  drop  of  k"  toward  zero  as  the  material  ap¬ 
proached  maximum  magnetization  was  not  observed  for  this  material. 

The  results  of  n"  and  k"  vs  y4nM/u  for  various  values  of  u^/ui 
are  shown  for  TT1-390  and  MMF-10  in  Figs.  16  through  20.  These  include 
higher  values  of  cj^/u  than  could  be  obtained  for  YIG  at  5.  5  GHz  and  at 
room  temperature  and  above.  However,  since  they  are  intrinsically  low- 
loss  materials,  it  was  not  possible  to  go  to  as  low  values  of  u^/u  as  with 
the  Dy-doped  YIG.  In  Fig.  16  the  full  hysteresis  loop  is  shown  for  p"  and 
k".  In  the  remaining  figures  only  one  branch  is  shown,  the  other  branch 
being  its  mirror  image.  At  large  values  of  p"  is  found  to  have  a 

double  peak.  This  results  from  a  phenomenon  first  observed  by  Le  Craw 
and  Spencer!^  in  which  xil  peaks  at  negative  magnetizations 

rather  than  at  zero.  They  referred  to  this  as  "anomalous  ferrimagnetic 
resonance"  and  gave  a  qualitative  explanation  of  it. 

For  both  TT1-390  and  MMF-10,  we  found  that  |x"|  peaked  at 
some  value  of  4irM  between  zero  and  saturation  as  expected  theoretically. 

As  Uj yj/u  dropped  below  1. 0,  k"  becomes  much  less  than  p". 

Both  these  materials  also  showed  substantial  hysteresis  in  p." 
and  k"  vs  4ttM.  The  character  of  their  hysteresis  is  different,  however. 
For  magnetization  decreasing  from  positive  to  negative,  the  maximum  n" 
occurs  for  negative  magnetization  in  TT1-390  and  positive  magnetization 
in  MMF-10.  The  maximum  of  |  k"  [  occurs  for  positive  magnetizations  in 
TT1-390  and  negative  magnetizations  in  MMF-10.  This  is  merely  a  mani¬ 
festation  of  the  fact  that  their  domain  structures  behave  differently  as  the 
magnetic  hysteresis  loop  is  traversed. 

In  Figs.  21  through  25,  pi"  is  plotted  as  a  function  of  for 

4ttM  =  0  and  4ttM  -  2/3(4rrMs)  for  each  of  the  three  materials  measured. 

The  magnetization  2/3(4rrMs)  is  assumed  to  be  roughly  the  remanence  value- 
The  values  of  k"  vs  at  2/3(4rrMg)  and  at  their  peaks  are  also  plotted. 
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Fig.  16  hi"  and  k"  vs  v4itM/u  for  TTI-390  at  5.  5  GHz 
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Fig.  19  (1"  vs  \4irM/u  for  MMF-IO  at  5. 


0.03 
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Fig.  20  k"  vs  y4ttM/u  for  MMF-10  at  5.5  GHz 
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Fig.  22 


H-"  vs  u  /u  for  TT1-390  at  5.  5  GHz 


E.  Measurements  of  p" 

1  1  Z 


Measurements  of  fi"  were  made  on  the  same  system  as  y.1  using 
the  change  in  reflection  coefficient  of  the  cavity  between  9400 Oe  and  the 
measurement  field.  The  results  for  Y1G  with  1  percent  Dy  and  MMF-10 
are  shown  in  Figs.  26  and  27.  The  ranges  of  for  these  two  materials 

do  not  overlap  significantly. 

For  all  measured  values  of  for  both  materials,  has  a 

single  peak  at  4trM  =  0.  However,  the  behavior  of  p"  compared  to  p"  as 
4irM  approaches  remanence  and  higher  magnetizations  is  quite  different 
for  the  two  materials.  For  the  Dy- doped  YIG,  p"  levels  off  near  or  some¬ 
what  below  remanence,  while  p^  continues  to  drop  rapidly  as  saturation  is 
approached.  In  this  region  p"  is  less  than  p".  For  MMF-10  both  p"  and 
p"  drop  rapidly  as  magnetization  increases  but  the  drop  in  p"  is  not  as 

2  Z 

sharp  and  at  high  magnetization  p"  is  greater  than  p".  It  is  not  clear  if 

z 

these  differences  are  due  to  the  different  u^/u  values  for  the  two  materials 
or  due  to  an  inherent  difference  in  the  materials. 

There  are  two  effects  which  contribute  to  p^  in  saturated  samples 

which  should  be  mentioned  in  any  discussion  of  p".  One  of  these  is  Larmor 

1112  z 
excitation  ’  by  inhomogeneities,  the  other  is  subthreshold  parallel 

pumping.*^  The  Larmor  excitation  results  because  the  local  magnetization 
and  rf  field  are  not  exactly  parallel  even  though  the  average  magnetization 
and  rf  field  may  be  parallel.  This  lack  of  local  parallelism  is  caused  by 
inhomogeneities  such  as  porosity  and  crystalline  anisotropy.  In  dense  poly¬ 
crystalline  materials  the  effect  is  very  small  and  only  noticeable  for  applied 
* 

dc  fields  which  place  the  sample  in  the  vicinity  of  ferromagnetic  resonance. 

The  subthreshold  parallel  pump  absorption  results  from  rf  magnetic 
field  pumping  of  standing  spinwaves  which  are  at  half  the  pump  frequency. 
These  spinwaves  have  a  net  moment  along  the  axis  of  their  precessional 
cone  and  can  couple  to  a  uniform  rf  field.  The  moment  is  a  result  of  the 
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Fig.  27  n"  vs  y-*"M/u  for  MMF-10  at  5.5  GHz 


elliptical  precessional  motion  and  the  coupling  mechanism  is  the  same  as 
that  which  causes  the  parallel- pump  instability  at  large  rf  fields. 

The  curves  in  Fig.  28  have  been  computed  using  the  theoretical 
equations  of  Ref.  13.  The  exchange  constant  of  YIG  has  been  assumed  for 
these  curves;  this  is  valid  assumption  only  for  YIG.  Nevertheless,  the 
results  should  be  a  reasonable  approximation  for  other  garnets  and  spinels 
whose  magnetization  differs  from  that  of  YIG  but  whose  Curie  tempera¬ 
tures  are  close  to  the  Curie  temperature  of  YIG.  The  curves  in  Fig.  28 
give  the  subthreshold  p"  as  a  function  of  frequency  for  several  values  of 
4rrMg.  It  has  been  assumed  that  the  temperature  is  room  temperature, 
that  the  sample  is  saturated,  and  that  the  internal  dc  field  is  zero.  It  is 
necessary  to  plot  p"  in  this  fashion  since  it  is  not  a  function  of  only  the 
ratio  For  values  of  w^/u  <  0.  8,  the  p^'  is  in  the  low  10  or  in 

the  10"  5  region.  When  u^/w  becomes  small  this  subthreshold  absorption 

may  be  one  of  the  dominant  contributions  to  p"  and  p"  in  the  partially  mag- 

z 

netized  state. 

In  comparing  values  p"  and  p^  in  the  demagnetized  state  for  the 
same  values  (where  they  are  presumably  equal),  it  should  be  remem¬ 

bered  that  will  have  errors  of  a  few  percent  due  to  errors  in  magnetiza¬ 
tion  measurement  and  deviation  of  y  from  2.  8  MGHz/Oe. 

F.  Low  Magnetic  Losses 

Certain  ferrimagnetic  materials  have  magnetic  losses,  p",  that 
are  much  less  than  0.  001  in  the  demagnetized  state,  providing  u^/u  is  not 
close  to  one.  However,  values  of  p"  as  low  as  a  few  times  10"4  are  difficult 
to  measure  accurately  by  cavity  perturbation  techniques,  because  a  large 
sample  is  required  in  order  to  produce  a  measurable  change  in  cavity  Q. 

If  the  sample  is  made  too  large,  the  results  of  perturbation  theory  are  not 
valid.  For  large  samples  an  exact  solution  of  the  new  cavity  mode  is  neces¬ 
sary.  This  can  be  done  only  for  certain  specific  geometries.  We  have 
chosen  to  use  a  rectangular  cavity  containing  a  thick  slab  that  fills  the  entire 
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Fig.  28  Subthreshold  n”  vs  Frequency 
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cross  section  at  one  end  wall.  A  diagram  of  the  cavity  is  shown  in  Fig.  29. 
The  method  of  calculation  of  resonant  frequency  and  cavity  A  as  a  function 
of  the  dielectric  constant,*,  and  complex  permeability,  n,  is  outlined  below. 


V 


Consider  a  slab  of  thickness  <5  in  a  rectangular  waveguide  cavity  of  length  L- 
The  cutoff  frequencies  of  the  waveguide  are: 

f  =  in  the  empty  part 

=  - - -  in  the  filled  part 

2a  1 « fi 

where  a  is  the  width  of  the  waveguide.  The  propagation  constants  in  the 
two  parts  of  the  waveguide  are  given  by: 


Y 


o 


in  the  empty  part. 


-  <T>2 


in  the  filled  part. 
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Fig.  29  Rectangular  TE1Q2  Cavity  for  Measuring  Very  Low 


The  characteristic  impedance  of  each  part  of  the  waveguide  is: 


ZQ  =  \J  i  -  (-j£)2  in  the  empty  part, 

Zf  =  \fl-  (-y  )2  in  the  filled  part. 

The  cavity  is  made  by  placing  a  short  at  each  end  of  the  waveguide.  Standing 
at  the  interface  between  the  ferrite-filled  and  empty  waveguide,  and  iooking 
down  the  empty  waveguide  toward  one  short,  one  sees  an  impedance 

Z  =  ZQ  tanh  (v(L-6)) 

Looking  back  into  the  ferrite-filled  waveguide,  one  sees  an  impedance 
Z  =  Zf  tanh  (-  y6)) 

Equating  these  two  values  for  the  impedance  gives 
ZQ  tanh  (y(L-<5))  +  Zj.  tanh  (y<5)  =  0 

Using  this  equation  and  the  preceding  equation  for  ZQ,  Z^,  yq>  Yp  fQ> 
and  fj.  one  can,  with  the  aid  of  a  digital  computer,  solve  for  any  one  of 
< ,  p,  L,  and  the  resonant  frequency  f  in  terms  of  the  other  three.  Natu¬ 
rally  if  p  is  complex  the  resonant  frequency  must  also  be  complex,  its 
imaginary  part  being  f"  =  f'/2Q 

In  practice  the  procedure  used  is  the  following.  The  frequency  of 
the  empty  cavity  is  first  measured  and  used  to  compute  the  effective  length 
of  the  cavity  to  the  accuracy  of  the  frequency  measurement.  Then  the  sample 
is  inserted  and  a  large  dc  magnetic  field  is  applied  parallel  to  the  rf  mag¬ 
netic  field  in  the  sample.  This  makes  p  =  1  and  allows  the  effective  di¬ 
electric  constant,  «,  to  be  computed.  This  effective  dielectric  constant 
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may  be  slightly  less  than  the  real  dielectric  constant  since  there  must  be 
some  very  small  air  gap  between  the  sample  and  broad  wall  of  the  wave¬ 
guide.  This  gap  acts  as  small  capacitance  in  series  with  the  large  capaci¬ 
tance  of  the  high  dielectric  constant  ferrite  and  the  net  effective  dielectric 
constant  is  given  by 


eff 


=  «  (1 


d« 


where  d  is  the  gap  thickness  and  b  is  the  height  waveguide,  0.  62  in.  at 
5.  5  GHz.  If  d  =  0.  002  in.  and  «  =  16  the  correction  factor  is  0.  95.  In  ad¬ 
dition,  the  measurement  of  the  effective  dielectric  constant  on  the  sample 
in  question  automatically  takes  into  account  any  variation  in  dielectric  con¬ 
stant  from  sample  to  sample. 


The  dc  field  is  then  reduced  to  zero  and  the  sample  demagnetizes 

in  its  own  demagnetizing  field.  The  resonant  frequency  and  cavity  Q  are 

again  measured  Using  the  change  in  f'/Q  for  f"  in  the  preceding  equations, 

we  can  calculate  m- '  and  p".  While  this  technique  has  been  used  successfully 

-4 

to  measure  cases  of  p"  <  1 0  ,  there  is  one  precaution  that  must  be  taken. 

We  have  neglected  the  skin  losses  in  the  waveguide  by  using  the  change  in 
(f'/Q)  from  high  to  zero  field  for  f.  In  doing  this,  we  have  assumed  that 
the  skin  losses  are  unchanged  from  high  to  zero  field  and  hence  could  be 
ignored.  In  fact,  this  is  not  true,  since  p'  and  therefore  resonant  frequency 
change,  the  skin  depth  and  losses  also  change.  Furthermore,  the  rate  of 
energy  loss  into  the  waveguide  walls  adjacent  to  the  material  depend  directly 
on  the  dielectric  constant  and  permeability  of  the  adjacent  material.  This 
also  causes  a  change  in  Q  from  high  to  zero  field.  Finally,  the  relative 
energy  stored  in  the  two  parts  of  the  cavity  changes  as  p1  changes  and  causes 
a  further  change  in  Q.  All  these  changes  can  be  calculated  in  straightforward, 
although  tedious,  manner  since  the  field  patterns  can  be  obtained  exactly 
throughout  the  cavity. 

The  error  in  pl!  is  given  by  (9  log  Q/9  log  u-  1)  Ap'/2Q  where  Ap'  is 
the  change  in  p’  from  high  to  zero  field  and  9  log  Q/9  log  w  is  the  logarithmic 
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change  in  Q  with  frequency  as  **'  is  varied.  For  a  TEjQ2  cavity  made  of 
1. 37  X  0.  62  in.  copper  waveguide  2.  43  5  in.  long  whose  resonant  frequency 
is  5.  5  GHz,  the  calculated  values  of  the  factor  (3  log  Q/3  log  u-  1J/2Q 
are  given  below  for  various  thicknesses  of  ferrite  slab  with  an  t  =  15. 


Thickness 

{3  log  Q/3  log  u  - 

0.  20 

-0.  3  X  10'4 

0.  040 

0.  2  X  10*4 

0.  060 

0.  7  X  10‘4 

0.  080 

0.  6  X  10'4 

0.  100 

2.  7  X  10‘4 

0.  120 

3.  8  X  10‘4 

Since  Ap.'  is  typically  0.  01  to  0.  1  for  very  low  loss  materials  with  u^/w  «  1 
the  error  in  will  generally  be  insignificant.  However,  as  the  slab  thick¬ 
ness  increases  and  more  of  it  is  in  the  region  where  the  electric  field  is  not 
zero,  its  large  dielectric  constant  causes  substantial  changes  in  the  stored 
energy  and  wall  loss.  Thus  this  factor  can  be  seen  to  increase  rapidly. 

Some  care  must  be  taken  not  to  use  an  excessively  large  sample  unless  n" 
is  small  enough  to  require  it  and  An'  is  small  enough  to  allow  it. 

We  have  made  preliminary  measurement  on  several  materials  at 

frequencies  from  5.5  to  16.5  GHz.  These  include  G-113,  TT1-390,  MMF-IO, 

G-800,  G- 1 004,  and  a  Raytheon  yttrium  iron  aluminum  garnet  (YFA)  with  a 

magnetization  of  4rrM  -  300  gauss.  The  results  are  summarized  in  Fig.  30. 

s 

Some  of  the  data  can  be  refined  further.  Many  more  points  are  needed  on 
several  of  these  samples  to  completely  evaluate  their  frequency  dependence. 
The  most  extensive  data  shown  is  for  G-113  covering  the  range  from  5.  5  to 
13  GHz.  The  results  obtained  at  16.  5  GHz  were  too  small  to  be  considered 
meaningful  for  this  material.  In  this  case  the  correction  described  abov'e 
for  the  change  in  cavity  wall  loss  with  frequency  overwhelms  the  material 
loss. 


57 


We  have  shown  earlier  in  this  report  that  for  u^/u  not  close  to 
1. 0  k"  «  k"  and  n"  is  not  a  strong  function  of  4irM.  We  can,  therefore, 
characterize  the  transverse  loss  as  a  function  of  <j^/u  by  the  single  value 
given  in  Fig.  30.  We  see  that  for  any  given  the  loss  of  G-113  is 

less  than  that  of  the  magnesium  manganese  ferrites  which  in  turn  is  less 
than  that  of  the  substituted  garnets. 
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Fig.  30  vs  uM/u  at  4»rM  =  0  for  Low  Loss  Materials 
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III.  THE  HIGH  POWER  PROPERTIES  OF  PARTIALLY  MAGNETIZED 
MATERIALS 


A.  Introduction 


Considerable  interest  in  microwave  devices  incorporating  partially 

magnetized  magnetic  materials  has  developed  in  recent  years.  To  obtain 

acceptable  high  power  performance,  it  is  necessary  that  concepts  exist 

which  can  be  used  during  the  design  to  optimize  peak  power  performance. 

However,  the  polycrystalline  materials  of  importance  in  device  design 

usually  possess  a  complex  domain  structure  when  partially  magnetized 

or  demagnetized,  and  the  characterization  of  the  high  power  properties 

is  extremely  complicated.  Virtually  no  systematic  study  of  the  high  power 

properties  of  such  systems  has  been  reported.  Experimental  work  to  date 

has  used  small  spherical  samples  magnetized  to  saturation  either  parallel 

or  perpendicular  to  the  direction  of  the  linearly  polarized  microwave  field. 

Furthermore,  the  theoretical  work  concerning  the  spinwave  instability 

processes  which  determine  high  power  performance  has  limited  applica- 

14 

bility.  The  Suhl  theory  of  instability  processes  applies  only  to  samples 
of  ellipsoidal  shape,  magnetized  to  saturation  along  an  axis  of  rotational 
symmetry,  and,  for  a  circularly  polarized  microwave  field,  perpendicular 
to  this  axis  rotating  in  the  proper  sense  to  excite  ferromagnetic  resonance. 
Schlomann's  theory  of  parallel  pumping*5  is  valid  for  a  linearly  polarized 
microwave  field  in  the  direction  of  saturation  for  ellipsoidal- shaped  samples. 
In  short,  the  published  work  concerning  high  power  properties  of  ferro¬ 
magnetic  insulators,  both  experimental  and  theoretical,  does  not  provide 
an  adequate  basis  for  optimizing  peak  power  performance  of  devices  in¬ 
corporating  partially  magnetized  materials. 

The  general  objective  of  this  investigation  was  to  initiate  a  program 

to  understand  the  parameters  which  determine  the  high  power  properties 

of  such  systems,  both  experimentally  and  theoretically.  The  approach  taken 

4 

in  this  initial  study  is  based  on  the  original  Polder-Smit  concept  that,  for 
partially  magnetized  systems,  the  local  internal  fields  are  determined  by 
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domain  distributions  and  domain  shapes,  in  much  the  same  way  that  they 
are  determined  by  the  sample  shape  for  saturated  materials.  Consequently, 
two  specific  related  objectives  have  been  followed:  1)  To  characterize 
the  peak  power  properties  of  saturated  materials  as  a  function  of  sample 
shape  and  microwave  field  configuration;  2)  To  relate  the  conclusions  for 
saturated  systems  to  results  for  partially  magnetized  samples  on  the  basis 
of  domain  shape  considerations.  The  problem  has  been  approached  essen¬ 
tially  from  a  microscopic  point  of  view  where  the  partially  magnetized 
system  is  viewed  as  a  collection  of  domains,  each  of  which  behaves  as  a 
saturated  sample  with  a  shape  defined  by  the  domain  geometry.  By  re¬ 
ducing  the  partially  magnetized  sample  to  a  collection  of  saturated  sample¬ 
like  regions,  the  problem  has  been  reduced  to  tractable  form  and  can  be 
treated  on  the  basis  of  theoretical  considerations  for  saturated  systems. 

This  approach  has  proved  quite  successful. 

It  is  important  to  indicate  what  types  of  high  power  measurements 
are  likely  to  yield  useful  and  pertinent  information  concerning  peak  power 
limitations  for  partially  magnetized  and  saturated  materials  of  interest. 

The  peak  power  handling  capability  of  microwave  magnetic  materials  is 
determined  by  the  magnitude  of  the  microwave  susceptibility  at  high  power 
levels.  Prior  to  this  work,  three  phenomena  had  undergone  extensive 
experimental  investigation:  saturation  of  the  main  resonance,  subsidiary 
absorption,  and  parallel  pump  absorption.  For  each  of  these  phenomena, 
the  nonlinear  coupling  which  exists  between  the  normal  modes  of  the  mag¬ 
netic  system  (spinwaves)  leads  to  a  parametric  excitation  of  certain  modes 
when  the  microwave  field  amplitude  exceeds  a  critical  threshold  value.  The 
possibility  of  such  excitations  must  be  taken  into  account  in  device  design. 
The  saturation  effect  consists  of  a  decrease  in  the  susceptibility  for  ma¬ 
terials  biased  near  ferromagnetic  resonance.  This  saturation  causes  a 
decrease  in  the  isolating  ability  of  a  material  used  in  a  resonance  isolator 
if  the  critical  power  is  exceeded.  Subsidiary  absorption  corresponds  to 
an  increase  in  the  susceptibility  for  materials  biased  well  below  ferro¬ 
magnetic  resonance,  and  leads  to  an  increase  in  insertion  loss  at  high  power 
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for  a  phase- shifting  device  using  the  material.  Parallel  pump  absorption 
is  similar  to  subsidiary  absorption,  but  occurs  for  the  microwave  field 
parallel  to  the  static  biasing  field. 

In  considering  these  and  related  phenomena,  their  relative  order 
of  importance  depends  on  the  specific  application  and  the  particular  high- 
power  requirements  of  the  device.  In  order  to  minimize  the  insertion  loss 
for  phase-shifting  devices  using  partially  magnetized  materials,  the  ma¬ 
terials  are  usually  chosen  with  a  magnetization  sufficiently  low  that  the 
local  internal  fields  are  much  less  than  required  for  ferromagnetic  reso- 
nace  at  the  frequency  of  operation.  Insofar  as  no  part  of  the  material  in 
the  remanent  device  is  biased  near  resonance,  saturation  of  the  resonance 
need  not  be  considered  in  characterizing  the  high- power  properties  of  such 
materials.  The  relevant  phenomena  of  importance  for  partially  magnetized 
magnetic  materials  are  related  to  subsidiary  absorption  and  parallel  pumping 
absorption,  which  do  occur  at  fields  less  than  required  for  resonance.  Both 
are  associated  with  so-called  first-order  spinwave  instability  processes 
In  the  actual  devices,  moreover,  more  complicated  microwave  field  con¬ 
figurations  may  exist,  and  first-order  instability  processes  for  such  con¬ 
figurations  must  also  be  considered.  The  experimental  approach,  therefore, 
has  been  to  examine  the  instability  threshold  for  samples  biased  far  below 
resonance  and  to  interpret  the  results  in  terms  of  first-order  instability 
processes. 

The  investigation  can  be  conveniently  divided  into  three  parts 
concerned  with  1)  the  characterization  of  peak  power  properties  as  a  func¬ 
tion  of  static  field,  sample  shape,  microwave  field  configuration,  and  state 
of  magnetization;  2)  the  theoretical  description  of  spinwave  instability 
processes  for  more  general  sample  shapes  and  pump  configurations  than 
previously  reported,  and  3)  the  application  of  the  theory  to  the  various 
experimental  results. 

Data  were  obtained  for  different  sample  shapes  and  pump  arrange¬ 
ments  which  were  expected  to  be  directly  related  to  the  magnetic  environments 
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for  partially  magnetized  systems:  1)  The  influence  of  gross  changes  in 
sample  shape  on  the  threshold  was  examined  by  considering  perpendicular 
pumping  for  axially  magnetized  rods,  spheres,  and  normally  magnetized 
disks.  The  rod  data  were  found  to  be  particularly  important,  since  many 
partially  magnetized  systems  consist  of  elongated  rod- shaped  domains; 

2)  Parallel  and  perpendicular  pumping  in  thin  disks  as  a  function  of  field 
and  pump  direction  were  also  examined.  Thin  disks  were  examined  in 
detail  for  two  reasons.  First,  many  remanent  devices  incorporate  ma¬ 
terials  with  slab  geometry,  so  that  the  results  are  likely  to  be  directly 
applicable  to  device  design.  Second,  the  domain  structure  for  thin  disks 
can  be  approximated  by  simple  models,  such  as  proposed  by  Kittel,^ 
and  the  results  for  partially  magnetized  disks  should  be  related  to  the 
theory  for  saturated  systems  according  to  the  above  reasoning.  Because 
the  actual  pump  environment  for  partially  magnetized  systems  is  likely 
to  be  more  complicated  than  simple  parallel  or  perpendicular  pumping, 
other  pump  configurations  were  also  investigated;  3)  Pumping  with  a  cir¬ 
cularly  polarized  field  rotating  in  the  proper  sense  to  excite  ferromagnetic 
resonance  and  in  the  opposite  sense  were  examined.  Pumping  with  a  linearly 
polarized  microwave  field  at  an  angle  between  zero  and  90*  to  the  static 
field  was  also  considered.  The  results  of  these  experiments  were  used 
in  two  ways.  Data  for  samples  magnetized  to  saturation  were  used  to  deter¬ 
mine  the  influence  of  sample  shape  and  pump  configuration  on  the  threshold. 
Data  for  partially  magnetized  samples  were  used  to  establish  a  connection 
between  domain  structure  and  the  threshold  fields,  guided  by  the  shape  in¬ 
formation  obtained  for  the  saturated  samples. 

The  general  theory  of  first-order  instability  processes  is  developed 
in  Appendix  1  Specific  equations  are  developed  which  apply  to  the  experi¬ 
mental  situations  described  above.  This  theory  provides  the  basis  for  a 
quantitative  comparison  of  the  various  experimental  results.  The  data  for 
saturated  systems  can  be  explained  in  terms  of  the  above  theory.  The 
theory  is  applicable  to  partially  magnetized  materials  if  domain  structure 
is  taken  into  account.  Additional  results  (for  oblique  pumping  in  spherical 
samples  as  a  function  of  angle  between  the  static  and  microwave  field)  are 
presented  in  Appendix  II,  these  results  are  in  agreement  with  the  above 
conclusions. 
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In  the  next  section,  the  instrumentation  which  was  developed  partly 
in  conjunction  with  this  program  is  described,  and  the  experimental  technique 
is  briefly  summarized.  A  large  amount  of  data  was  required  in  the  course 
of  the  investigation,  and  an  automated  system  to  determine  the  instability 
threshold  was,  therefore,  developed.  This  system  is  described  in  the  fol¬ 
lowing  three  sections;  the  experimental  results  are  presented  for  each  of 
the  above  experiments  and  discussed  in  terms  of  theory.  In  each  of  these 
presentations ,  the  emphasis  is  on  shape  and  pump  configuration  effects  for 
saturated  systems  and,  where  possible,  the  relation  to  domain  shape  effects 
on  hcrit  for  partially  magnetized  systems. 

B.  Experimental  Technique  and  Instrumentation 

Although  experimental  techniques  for  determining  the  instability 

17 

threshold  in  ferromagnetic  insulators  are  well  known,  several  aspects 
of  the  present  technique  and  the  instrumentation  developed  partly  in  con¬ 
junction  with  this  program  depart  significantly  from  previous  procedures. 

The  departures  consist  primarily  in  1)  automating  the  instrumentation 
used  to  determine  the  threshold  amplitude,  and  2)  introducing  sample  align 
ment  procedures  for  nonspherical  samples.  The  first  was  considered  neces¬ 
sary  because  of  the  large  amount  of  data  required  in  the  present  program. 

The  second  was  required  because  of  the  emphasis  of  peak-power  properties 
as  a  function  of  sample  shap  \  and  was  particularly  important  for  samples 
with  slab  geometry.  Both  of  these  will  be  considered  in  detail  below.  In 
discussing  the  automation  procedure,  we  will  briefly  review  the  basic  con¬ 
siderations  concerning  microwave  susceptibility,  cavity  fields,  and  micro- 
wave  power  which  are  necessary  in  applying  microwave  measurement 
techniques  to  the  understanding  of  physical  phenomena. 

The  basic  procedure  for  determining  the  first-order  instability 
threshold,  hcr^,  is  to  measure  the  imaginary  part  of  the  diagonal  com¬ 
ponent  of  the  effective  or  external  susceptibility  tensor,  which  is  conven¬ 
tionally  denoted  as  x  as  a  function  of  the  applied  microwave  field 
amplitude,  h,  at  the  sample  position.  Using  standard  cavity  perturbation 
techniques  reviewed  in  Ref.  17,  xg"  can  be  calculated  to  be 
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where  V  and  V  are  the  cavity  and  sample  volumes,  respectively,  r  and 
To  are  reflection  coefficients  for  the  cavity  with  and  without  the  magnetic 
specimen,  and  g  is  a  geometric  factor  related  to  energy  storage  in  the 
cavity  given  by 

«-i  [*-A]  <•« 

L  n  a  -1 

for  a  TEjgn  rectangular  cavity  where  d  is  the  cavity  length,  a  is  the 
cavity  width,  and  n  is  the  mode  index  corresponding  to  the  number  of 
half-wavelengths  along  the  length  of  the  cavity.  At  low  microwave  power 
levels,  the  susceptibility  xe"  is  independent  of  the  microwave  field  ampli¬ 
tude  h  and  Tg  is  constant,  for  fixed  static  field.  As  the  microwave  power 
level  is  increased  to  the  point  where  h  at  the  sample  position  exceeds  the 
threshold  value  for  the  excitation  of  spmwave  instabilities,  xe”  increases 
sharply.  Experimentally,  this  phenomenon  is  evidenced  by  an  increase 
in  the  reflection  coefficient  Tg  for  the  cavity  containing  the  sample.  The 
standard  experimental  procedure  is  to  measure  T g  as  a  function  of  the 
microwave  field  amplitude  h,  calculate  xe"  from  Eqs.  (17)  and  (18),  display 
Xg"  vs  h  graphically,  and  determine  hcr^  as  the  break  in  the  curve  which 
separates  the  constant  xe"  region  from  the  increasing  xe"  region.  Such  a 
procedure  is  extremely  time-consuming,  because  considerable  data  are 
required  to  determine  a  single  value  of  hcr^.  Values  of  hcrit  as  a  function 
of  a  number  of  parameters  (such  as  sample  shape,  sample  orientation, 
microwave  field  configuration,  and  static  field)  require  extensive  data. 

A  complete  point -by- point  characterization  of  high-power  properties  would 
be  impractical.  A  semiautomatic  instrumentation  system  has,  therefore, 
been  developed  to  simplify  the  above  procedure  considerably.  The  system 
allows  curves  of  log  Tg  to  be  plotted  as  a  function  of  h  .  Values  of  h  ^ 
can  be  easily  and  rapidly  determined  from  the  curves  generated. 
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The  system  is  designed  around  a  standard  reflection- cavity  micro- 
wave  spectrometer;^  the  microwave  instrumentation  is  indicated  sche¬ 
matically  in  Fig.  31.  The  basic  system  consists  of  a  2J51  magnetron  source 
operated  in  a  pulsed  mode,  which  feeds  a  TEjq^  rectangular  reflection 
cavity  through  a  series  of  high- power  attenuators  and  a  high- power  isolator. 
A  peak  power  of  1  kW  incident  upon  the  cavity  could  be  obtained.  Direc¬ 
tional  couplers  sample  the  incident  power  into  reference  arm  A  and  the 
reflected  power  into  cavity  arm  B.  A  commercial  self-balancing  micro- 
wave  power  meter  with  a  temperature- compensated  thermistor  mount  was 
used  to  monitor  the  average  power  level  in  the  reference  arm  and  to  pro¬ 
vide  a  voltage  output  proportional  to  this  power  for  the  x-axis  of  the  x-y 
recorder.  A  second  directional  coupler  in  the  reference  arm  feeds  the 
crystal  detector  through  one  side  of  the  microwave  switch.  The  detector 
monitors  the  reflected  power  in  the  cavity  arm  through  the  other  side  of 
the  switch.  Attenuators  are  provided  in  both  the  reference  and  cavity  arms 
to  adjust  the  signal  levels,  and  to  facilitate  point- by- point  measurements 
when  the  automatic  feature  is  not  desired.  The  electronic  circuits  provide 
appropriate  pulses  to  control  the  magnetron  output  and  to  drive  the  micro- 
wave  switch.  The  system  output  consists  of  a  dc  voltage  proportional  to 
the  logarithm  of  the  cavity  reflection  coefficient  Tg.  This  voltage  is  applied 
to  the  y-axis  of  the  recorder.  In  operation,  the  output  voltage  is  constant 
and  independent  of  the  power  level,  as  long  as  the  microwave  field  amplitude 
at  the  sample  does  not  exceed  the  instability  threshold  amplitude.  In  this 
field  region,  the  recorder  trace  consists  of  a  horizontal  line,  as  indicated 
on  the  left  side  of  the  x-y  recorder  plot  shown  in  Fig.  31.  As  the  field 
amplitude  exceeds  some  critical  value,  log  Tg  increases,  the  recorder 
curve  exhibits  a  break  and  then  increases  with  positive  slope.  The  break 
point  corresponds  to  the  instability  threshold. 

Before  considering  the  electronic  instrumentation  in  more  detail, 
it  is  necessary  to  indicate  the  various  functions  which  such  instrumentation 
must  perform  in  order  to  obtain  the  recorder  output  described  above.  The 
sequence  and  relative  timing  of  the  various  pulses  are  indicated  in  Fig  32. 
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Fig.  32  Schematic  Representation  of  the  Relevant  Pulse  Sequences 
for  Operation  of  the  Semiautomatic  System 
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Each  trigger  pulse  initiates  a  4psec  magnetron  pulse,  which  is  monitored 
by  the  reference  arm.  The  reflected  power  from  the  cavity  is  monitored 
by  the  cavity  arm.  The  pulse  shapes  for  this  reflected  arm,  shown  in 
Fig.  32,  indicate  the  typical  shapes  expected  when  the  microwave  field 
amplitude  is  above  threshold.  The  crystal  detector  is  switched  to  sample 
each  arm  on  alternate  pulses.  As  shown  in  Fig.  32,  the  detector  output 
will  correspond  to  the  incident  power  on  one  pulse  and  then  to  the  reflected 
power  on  the  next  pulse.  Information  concerning  susceptibility  can  be  ob¬ 
tained  from  the  reflection  coefficient  Tg(  which  is  equal  to  the  ratio  of  the 
reflected  pulse  amplitude  VQ  to  the  incident  pulse  amplitude  V^.  In  order 
to  convert  these  amplitudes  to  dc  levels,  two  sample-and-hold  circuits, 
with  appropriate  gates,  as  indicated  by  the  bottom  pulse  train  in  Fig.  32, 
must  be  used  to  sample  the  crystal  output  on  alternate  pulses.  The  quantity 

log  r  can  then  be  conveniently  obtained  by  feeding  the  sample-and-hold  - 

circuit  outputs  through  logarithmic  converters  and  then  to  a  differential 
amplifier.  This  quantity  (log  Tg  =  log  Vg  -  log  V^)  should  be  constant  as 
long  as  the  microwave  field  amplitude  is  below  the  threshoJd  value  and 
should  increase  as  the  threshold  amplitude  is  exceeded.  A  sharper  increase 
could  be  obtained  by  converting  the  voltage  output  to  Tg  using  an  anti- 
logarithmic  converter  after  the  differential  amplifier,  but  in  this  investi¬ 
gation  the  thresholds  were  quite  sharp  without  such  conversion. 

The  electronic  instrumentation  developed  to  derive  the  y-axis 
deflection  voltage  proportional  to  log  r  is  indicated  in  Fig.  33.  A  rate 
generator  synchronized  to  the  line  frequency  is  used  to  provide  trigger 
pulses  for  the  magnetron  pulse  generator,  the  pulse  generator  for  the 
sample-and-hold  circuit  gate  pulse,  and  a  bistable  divider.  The  bistable 
divider  drives  an  electronic  switch  which  applies  the  sample-and-hold 
gate  pulse  to  the  'A"  and  "B"  sample  and  hold  circuits  on  alternate  mag¬ 
netron  pulses,  respectively.  It  also  drives  the  microwave  switch,  so  that 
the  crystal  detector  output  applied  to  the  sample-and-hold  circuits  is  de¬ 
rived  from  reference  arm  A  when  the  sample-and-hold  gate  is  applied  to 
the  "A"  circuit,  and  from  the  cavity  arm  B  when  the  gate  is  applied  to  the 
"B"  circuit.  Both  voltages  are  derived  from  the  same  crystal  detector 
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output  so  that  the  tracking  error  which  would  result  from  using  different 
crystals  to  monitor  power  levels  in  the  reference  and  cavity  arms  is 
eliminated. 


To  determine  the  critical  threshold  field  from  the  break  point  in 

the  experimental  curve,  the  only  calibration  required  is  for  the  horizontal 

or  power  axis  in  terms  of  the  microwave  field  amplitude  in  the  cavity  at 

the  sample  position.  Calibration  of  the  peak  power,  P.  ,  incident  on  the 

me 

cavity,  in  terms  of  the  average  power  monitored  by  the  power  meter  in 
the  reference  arm  can  be  easily  obtained  by  mounting  the  thermistor  power 
probe  in  place  of  the  cavity,  noting  the  conversion  factor  between  this  reading 
and  the  reference  arm  power,  and  determining  the  magnetron  duty  cycle.  As 
reviewed  in  Ref.  17,  the  microwave  field  amplitude  hQ  can  be  calculated  as 
a  function  of  P^nc  by  considering  energy  transfer  and  storage  for  the  cavity 
and  the  microwave  field  distribution  inside  the  cavity.  The  result  is 
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where  r  ***  is  the  reflection  coefficient  of  the  empty  cavity  in  db.  The 

O  JL 

cavity  Q  without  the  sample, Oq,  and  rQ  are  measured  separately 


Data  were  obtained  for  a  particular  sample  configuration  by  gen¬ 
erating  x-y  recorder  curves  of  log  Tg  as  a  function  of  average  incident 
power  with  the  static  magnetic  field  as  a  parameter.  The  static  magnetic 
field  was  measured  using  a  Hall  effect  magnetic- field  probe  calibrated 
with  proton  nuclear  magnetic  resonance.  For  simplicity,  the  electro¬ 
magnetic  and  Hall  probe  are  not  indicated  in  Fig.  31.  Typical  data  are 
shown  in  Fig.  34.  A  horizontal  scale  calibration  for  h  ,  obtained  in  the 
manner  described  above,  is  also  indicated  in  that  figure.  The  instability 
threshold  field  amplitude  corresponds  to  the  point  where  each  curve  breaks 

upward.  Previously,  each  curve  shown  in  Fig.  34  would  have  to  be  obtained 

db 

from  time-consuming  point-by- point  measurements  of  as  a  function 

of  microwave  power  level.  All  the  above  curves  could  now  be  obtained  in 
a  matter  of  minutes. 
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Fig.  34  Sample  Recorder  Outputs  of  the  Cavity  Reflection 

Coefficient  as  a  Function  of  Microwave  Field  Ampli¬ 
tude.  At  successive  static  field  values,  the  curves 
were  shifted  upward  so  that  the  different  curves  could 
be  displayed  conveniently  on  the  same  figure. 
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The  second  aspect  of  the  experimental  technique  to  be  considered 
is  the  sample  alignment  procedure.  The  results  for  thin  disks  are  quite 
sensitive  to  slab  orientation.  In  order  to  ensure  accurate  sample  align¬ 
ment,  an  apparatus  was  designed  to  facilitate  accurate  sample  positioning 
with  respect  to  the  static  magnetic  field  direction.  A  photograph  of  this 
apparatus,  installed  on  the  magnet  pole  face  and  with  the  microwave  cavity 
in  position,  is  shown  in  Fig.  35.  The  cavity  is  mounted  near  the  center  of 
a  circular  lucite  plate  and  connected  to  the  microwave  system  by  a  section 
of  flexible  waveguide.  The  plate  assembly  is  mounted  in  the  gap  of  the 
electromagnet  so  that  the  polarization  direction  for  the  microwave  mag¬ 
netic  field  of  the  cavity  is  either  parallel  or  perpendicular  to  the  static 
field  direction,  depending  on  the  pump  orientation  desired.  The  plate  was 
pivoted  at  its  center  and  provided  with  two  adjustment  screws  spaced  90° 
apart  on  the  edge  of  the  disk  for  orientation  purposes.  The  cavity  plate 
was  held  under  tension  against  the  adjustment  screw  support  points  and 
center  pivot  point  with  rubber  bands  located  near  the  adjustment  screws, 
as  well  as  on  the  opposite  side  of  the  disk.  Sample  alignment  was  usually 
accomplished  at  low  power  with  the  X-13  klystron  microwave  source  shown 
in  Fig.  31.  By  minimizing  or  maximizing  the  reflected  power  from  the 
cavity  with  the  sample  biased  on  the  low-field  side  of  ferromagnetic  reso¬ 
nance,  depending  on  the  sample  shape  and  orientation,  accurate  alignment 
to  ±  0.  r  could  be  accomplished. 

Alignment  was  extremely  critical,  particularly  for  disk- shaped 

samples.  To  illus7rate  this  point,  data  for  perpendicular  pumping  are 

shown  in  Fig.  36  for  a  normally  magnetized  0.  020  in.  thick  and  1  cm  diam 

disk  of  polycrystailine  YIG  at  9.2  GHz.  Data  were  obtained  with  the  static 

magnetic  field  at  angle  of  0,  5,  and  15°  to  the  disk  normal.  The  threshold 

field  h  .  is  shown  as  a  function  of  external  static  field.  These  data  show 
cnt 

that  orientation  has  a  significant  effect  of  hcr-t,  particularly  for  the  par¬ 
tially  magnetized  state  (corresponding  to  static  fields  less  than  about  17  50G 
in  Fig.  36).  At  a  nominal  field  of  500  G,  hcrit  is  reduced  by  almost  a  factor 
of  two  by  the  misalignment.  This  reduction  corresponds  to  a  factor  of  four 
reduction  in  the  peak-power  limitation  for  the  material. 
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Fig.  3  5  Photograph  of  the  Alignment  Apparatus  Utilized  in  Obtaining 
Data  on  Disk-Shaped  Samples.  The  two  knurled  knobs  allow 
the  sample  to  be  adjusted  normal  (or  parallel)  to  the  static 
field  within  ±0.1*. 
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C.  Sample  Shape  and  the  Perpendicular  Pump  Threshold 


Large  changes  in  sample  shape  have  a  profound  influence  on  h 

In  Fig.  37  the  threshold  field  hcr^t  is  shown  as  a  function  of  static  external 

field  Hq  for  perpendicular  pumping  in  spheres,  axially  magnetized  40-mil 

diam  rods,  and  normally  magnetized  10-mil  thick  slabs  of  polycrystalline  YIG. 

The  data  were  obtained  at  9.  2  GHz  and  room  temperature.  Data  are  shown 

both  before  and  after  a  1  hr  anneal  at  1000*C.  The  threshold  for  the  axially 

magnetized  rod,  indicated  for  unannealed  and  annealed  rods  by  the  open  and 

solid  circles  respectively,  decreases  with  increasing  field,  has  a  minimum 

of  about  50e  at  a  static  field  of  about  1300  0e,  and  then  increases  sharply. 

The  results  for  the  spherical  samples,  indicated  for  the  unannealed  and 

annealed  materials  by  the  open  and  solid  triangles  respectively,  show  a 

slow  increase  of  H  less  than  about  600  Oe,  then  a  decrease  to  a  50e  mini- 
o 

mum  at  a  static  field  of  1900  0e  and  then  a  sharp  increase.  Except  for  the 
shift  to  higher  static  field  and  the  decrease  in  h  below  600 Oe,  the  ex¬ 
perimental  curve  for  the  spherical  sample  is  quite  similar  to  that  for  the 
axially  magnetized  rod.  The  data  for  the  10-mil  thick,  normally  magnetized 
slab  are,  however,  quite  different.  In  general,  hcrit  is  significantly  larger 
for  the  slab  than  for  either  the  sphere  or  the  rod  geometry.  The  threshold 
increases  from  12  0e  at  zero  external  field  to  a  maximum  of  15  to  18  0e  at 
an  external  field  of  1700  0e,  decreases  sharply  to  a  local  minimum  of  90e 
at  about  23000e,  exhibits  a  second  local  minimum  at  2900 Oe,  and  then  in¬ 
creases  sharply  at  higher  field.  The  discussion  below  will  treat  the  three 
important  features  of  the  data  evident  from  Fig.  37:  1)  The  influence  of 
geometry  on  h^rit  for  saturated  samples,  2)  Behavior  of  hcrit  for  partially 
magnetized  samples,  3)  Annealing  effects. 

It  is  well  known  that  the  internal  magnetic  fields,  both  microwave 

and  static,  are  influenced  by  geometry  and  that  the  action  of  these  fields  on 

the  local  magnetization  determine,  in  part,  the  high-power  properties  of 

the  material.  The  increase  in  field  value  for  the  minimum  h  ..  from  rod 

crit 

to  sphere  to  slab  data. is  connected  with  an  increase  in  the  demagnetizing 
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factor  N""  in  the  static  field  direction,  in  going  from  a  rod  shape  (Ng  »  0) 
to  a  sphere  (N  *  1/3)  to  a  slab  (N  »  1).  The  corresponding  increase  in 
the  minimum  h  from  4.60e  (rod)  to  4.  90e  (sphere)  to  8.  50e  (slab)  is 
associated  with  a  decrease  in  the  efficiency  with  which  the  microwave  field 
can  excite  the  uniform  precession  mode. 

These  points  can  be  demonstrated  more  quantitatively  by  conridering 

theoretical  curves  of  h  ..  vs  static  field  H  for  the  sample  shapes  of  interest. 

crit  o 

The  theoretical  curves  for  comparison  with  the  data  in  Fig.  37  as  well  as  for 
comparison  with  the  data  to  follow  in  the  next  two  sections,  were  obtained 
according  to  the  following  procedure:  1)  Data  on  spherical  samples  and  the 
theory  given  in  Appendix  I  were  used  to  evaluate  the  spinwave  linewidth  as 
a  function  of  spinwave  number  k  and  the  spinwave  angle  0^,  2)  The  depen¬ 
dence  of  AH^  on  k  and  0^  was  then  incorporated  into  the  theory  for  the  rod 
and  disk  configurations  and  hcr^t  vs  Hq  was  calculated  for  these  geometries. 

In  other  words,  sphere  data  were  used  as  a  standard  to  obtain  AH^  values 
which  were  then  utilized  to  obtain  theoretical  curves  for  all  other  configura¬ 
tions.  A  comparison  of  the  data  for  all  configurations  on  the  basis  of  derived 
AH^(k,  0^)  values  might  seem  to  be  a  more  consistent  approach.  Inasmuch  as 
AHk  is  actually  a  derived  quantity,  a  direct  comparison  of  threshold  field 
values  seemed  to  be  more  appropriate. 

The  AH^ik,  0^)  results  for  polycrystalline  YIG  were  obtained  from 
parallel  and  perpendicular  pumping  data  summarized  in  Fig.  38,  using  the 
theory  given  in  Appendix  I.  The  results  of  this  determination  are  shown  in 
Figs.  39a  and  b,  where  AH^  is  shown  as  a  function  of  k  and  0^  for  both 
parallel  and  perpendicular  pumping. 

It  is  important  to  realize  that  the  portion  of  the  threshold  field  vs 
static  field  curve  of  primary  interest,  the  portion  for  static  field  less  than 
that  at  the  minimum  hcr^t  position  HQmin,  corresponds  to  k  >  0,  the  spin- 
wave  linewidth  is  given  by  the  results  in  Fig.  39b,  which  indicate  that  AH^ 
is  essentially  an  increasing  function  of  k  and  does  not  differ  appreciably 
from  parallel  to  perpendicular  pumping,  even  though  0^  is  quite  different 


78 


(•0)^°M  CH3U  aiOHS3dHl 


Sl3 
*  c 
oJ 

*?%  «-* 

<L>  rtj 

*« 
r2  u 

g£ 

If 

am 

.H  "o 
£  e 
r)  <g 
*•» 

c/3  r* 

v  • 
jC  M 

•_  b. 

®  v.  , 

c  o  w 

.2  «  S' 

7!  d.*— 
c  £  cfl 

3  S  u 

&■«  a 

M  « *5 

m  o  5 
to  •-  o 
b  Q 
■*■•  ^ 
‘Cjs  a 

j&S 

o  ~  u 
’C  4)  J2 

C  3 
■D  C  O 

0<  = 
£0C 
oi  jr  0) 
O  -*-»  c. 

,  $•* 
j=  r  <u 


79 


for  the  two  configurations.  Below  the  minimum  h  ^  position  on  the 
h  ..  vs  H  curve,  corresponding  to  H  <  H  min,  these  results  can  be 
used  to  calculate  theoretical  h  ^  values  for  the  other  field  configurations 
and  sample  shapes  of  interest.  The  situation  for  >  HQmin  is  more 
complicated.  In  this  field  region,  the  unstable  spinwave  generally  cor¬ 
responds  to  k  »  0,  and  the  AH^  determinations  summarized  in  Fig.  39a 
must  be  used  to  calculate  the  theoretical  curves.  These  determinations 
indicate  that  AH^  is  not  a  single-valued  function  of  0^  for  k  -  0,  but  de¬ 
pends  on  the  pump  direction.  For  parallel  pumping,  AH^  has  a  maximum 
at  0^  83  50°.  For  perpendicular  pumping  AH^  increases  with  decreasing  9^. 

The  spinwave  linewidth  AH^  is  generally  considered  to  be  a  measure  of  the 
spinwave  relaxation  rate,  a  material  parameter,  and  not  expected  to  be  a 
function  of  geometry.  Although  this  expectation  is  fulfilled  for  k  >  0  (Fig.  39b), 
it  does  not  appear  to  be  satisfied  for  k  =  0,  perhaps  because  the  theory  is  not 
strictly  applicable  for  k  =  0.  Even  though  the  results  are  approximately 

corrects  the  AH,  do  not  correspond  to  spinwave  (u/2Y  k  =  0,  9.  )  but  to  mag- 
K  18  K 

netostatic  mode  excitations.  The  coupling  to  these  modes  should  be  strongly 

influenced  by  the  microwave  field  direction.  Such  considerations  are  beyond 

the  scope  of  this  report  and  are  not  directly  related  to  the  present  study.  For 

the  sake  of  comparison,  the  AH,  vs  9,  results  shown  in  Fig.  39a  will  be  used 

*»  __  •  __ 

to  calculate  theoretical  h  ..  values  for  H  >  H  using  the  AH.  determina- 

crit  o  o  ®  k 

tions  for  parallel  pumping  when  parallel- pump  curves  are  to  be  calculated 
and  for  perpendicular  pumping  when  perpendicular- pump  curves  are  to  be 
calculated.  The  meaning  of  theoretical  curves  for  Hq  >  Homin,  however, 
is  to  be  interpreted  in  light  ot  the  above  considerations. 

Using  this  procedure,  AH^  determinations  from  the  sphere  data 
of  Fig.  38  have  been  used  to  calculate  theoretical  threshold  curves  for 
perpendicular  pumping  with  hnearty  polarized  excitation  in  axially  mag¬ 
netized  rods  and  normally  magnetized  slabs.  In  Fig.  40  the  data  for  the 
annealed  materials  are  shown  again  for  comparison  with  theory.  No  theo¬ 
retical  curve  is  shown  for  the  sphere,  because  those  data  were  used  to 
evaluate  AH^.  The  rod  data  are  in  close  agreement  with  the  calculated 
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Comparison 


curve,  which  was  obtained  with  no  adjustable  parameters.  The  difference 
is  <  1  Oe  and  is  probably  related  to  the  effect  of  annealing  on  the  threshold. 
From  Fig.  37,  it  is  clear  that  annealing  serves  to  reduce  hcrjt  for  the  rod 
and  further  annealing  might  reduce  the  threshold  still  further  and  improve 
the  agreement  with  the  theory.  The  comparison  between  the  data  and  theory 
for  the  10-mil  thick,  1-cm  diam  disk  are  also  shown  in  Fig.  40.  The  cal¬ 
culated  curve  does  not  extend  below  1700  0e  since  the  sample  is  not  saturated 
at  lower  field,  and  the  theory  is  not  directly  applicable.  The  agreement  is 
not  as  good  as  it  is  for  the  rod  geometry.  The  theoretical  minimum  is  lower 
than  the  experimental  minimum,  and  the  calculated  curve  does  not  exhibit 
the  double  minimum  indicated  by  the  data. 

The  general  increase  in  the  minimum  hcrit  value  in  progressing 
from  the  rod  to  the  sphere  to  the  disk  data  is  associated  with  the  efficiency 
with  which  the  microwave  field  pumps  the  uniform  precession  amplitude. 

For  samples  with  rotational  symmetry  about  the  static  field  direction,  the 
uniform  precession  amplitude  due  to  the  Larmor  component  of  the  pump 
field  is  giv-.-n  by 


hoeiwt 

ao  "  2(w  -<j) 
o 


where 


u  =  y  (  H  +  (N.-N  )  4irM  J  (20) 

O  O  t  Z  13 

is  the  natural  precession  frequency  for  the  uniform  precession  mode.  In 

the  above  expression,  N  and  N  are  axial  and  transverse  demagnetizing 

factors,  respectively,  h  is  the  microwave  field  amplitude,  and  4*M  is 

o  s 

the  saturation  induction.  In  progressing  from  the  rod  with  Nz  *N  =  1/3, 
to  the  disk  with  Nz  55  1  and  ~  0,  the  uniform  mode  amplitude  decreases 
(Nj-Nz  decreases  from  |  to  -1).  In  effect,  the  efficiency  with  which  the 
microwave  field  pumps  the  uniform  mode  decreases  and  results  in  an  in¬ 
crease  in  the  h  ^  required  for  instability  to  occur. 
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The  differences  between  the  slab  data  and  the  theory  are  probably 
associated  with  magnetostatic  mode  resonances  which  are  excited  in  the 
disk  because  the  pump  field  is  not  entirely  uniform  over  the  sample.  Similar 
to  the  uniform  precession,  these  modes  are  also  coupled  to  the  spinwave 
modes  and  can  lead  to  instability.  For  a  normally  magnetized  slab,  mag¬ 
netostatic  mode  resonances  occur  at  static  field  less  than  that  required  for 
ferromagnetic  resonance.  It  is  suggested  that  the  dip  in  h  near  HQ~2300Oe 
for  the  disk  data  in  Figs.  37  is  due  to  magnetostatic  mode  interactions.  The 
strong  coupling  between  the  external  microwave  field  and  the  various  magneto¬ 
static  resonances  is  presumably  accompanied  by  a  decrease  in  coupling  to  the 
uniform  mode.  This  decrease  may  be  partially  responsible  for  the  large  hcr^t 

values,  in  contrast  to  the  theory,  near  H  =  H  min. 

o  o 

The  increase  in  hcr-t  with  Hq  at  field  strengths  insufficient  to 

saturate  the  material  is  clearly  evident  for  the  sphere  data  for  Hq  <  600 Oe 

and  for  the  slab  data  for  H  <  1700  0e.  The  rod  is  saturated  for  H  »  0,  so 

o  o 

that  no  corresponding  region  is  evident  from  the  rod  data.  This  behavior 
of  hcr^t  below  saturation  can  be  qualitatively  understood  by  considering  the 
parallel  and  perpendicular  pumping  data  shown  in  Fig.  38  for  spherical 
samples.  As  Hq  is  reduced  below  600 Oe  the  sample  demagnetizes,  and 
the  magnetization  direction  in  individual  domains  is  no  longer  constrained 
to  lie  in  the  applied  field  direction.  This  demagnetization  results  in  a  re¬ 
duction  in  the  hcrit  for  perpendicular  pumping  and  an  increase  in  hcrit  for 
parallel  pumping.  The  thresholds  for  the  two  configurations  approach  a 
common  value  in  the  limit  Hq  =  0,  corresponding  to  an  isotropically  de¬ 
magnetized  material.  The  break  in  each  of  the  two  curves  near  HQ  =  600Oe 
is  indicative  of  the  onset  of  demagnetization  processes.  The  demagnetiza¬ 
tion  is  not  completely  isotropic  so  there  are  deviations  from  the  expected 
convergence  at  low  values  of  static  field  as  indicated  by  the  data  for 
Hq  <  200 Oe.  The  transition  between  parallel  and  perpendicular  pumping 
for  both  saturated  and  partially  magnetized  samples  is  considered  further 
in  Appendix  II,  in  conjunction  with  the  oblique  pumping  study. 
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As  shown  in  Fig.  37,  the  relatively  short  1  hr  anneal  at  1000‘C 
had  a  significant  effect  on  h  particularly  for  the  slab  data.  The  rod 
exhibited  a  small  decrease  in  hcr^t  after  annealing.  The  results  for  the 
sphere  are  virtually  identical  before  and  after  the  anneal.  For  the  normally 
magnetized  slab,  h  t  increased  after  annealing.  In  all  cases,  the  general 
shapes  of  the  hcrit  vs  Hq  curves  were  preserved  after  annealing.  The  fact 
that  annealing  has  a  relatively  large  effect  only  on  the  results  for  the  slab 
indicates  that  the  phenomenon  is  associated  with  surface  damage  or  sur¬ 
face  strain.  Consequently,  surface  damage  during  grinding  should  have 
the  largest  effect  on  the  magnetic  properties  for  thin  slabs,  it  should  be 
noticed  that  the  sign  of  the  change  upon  annealing  is  not  the  same  for  the 
two  sample  shapes  where  the  effect  is  noticeable.  Annealing  decreases 
^crit  *°r  axially  niagnenzed  rod  and  increases  hcr^t  for  the  slab.  This 
annealing  characteristic  may  be  associated  with  the  specific  nature  of  the 
surface  strain  interaction  with  the  magnetization.  In  the  case  of  the  rod, 
the  static  field  and  magnetization  are  parallel  to  the  sample  surface,  while 
for  the  slab,  they  are  normal  to  the  sample  surface. 

D.  Parallel  and  Perpendicular  Pumping  in  Thin  Disks 

In  this  study  of  the  influence  of  sample  shape  on  the  instability 
threshold,  one  possibility  for  obtaining  additional  information  is  to  apply 
the  static  field  in  different  directions  with  respect  to  the  sample  principal 
axes.  (Rods  and  disks  may  be  considered  to  be  limiting  cases  of  ellipsoids 
of  revolution.  )  In  spherical  samples,  however  ,  all  directions  are  equivalent. 
For  rods  and  disks,  on  the  other  hand,  parallel  and  perpendicular  pumping 
can  be  accomplished  for  a  number  of  different  field  arrangements  relative 
to  the  sample  axes.  Obtaining  data  for  different  field  arrangements  is 
equivalent  to  interchanging  the  various  sample  shape  demagnetizing  factors 
and,  as  such,  represents  a  useful  technique  to  further  study  the  effect  of 
sample  shape  on  h  The  present  work  was  performed  on  thin  1  cm  diam 

disks  for  different  field  orientations  as  a  function  of  disk  thickness  and 
annealing  history. 
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Five  different  field  arrangements  are  possible  if  the  experiment 
is  restricted  to  parallel  or  perpendicular  pumping.  These  arrangements 
are  indicated  schematically  in  Fig.  41.  The  two  top  diagrams  correspond 
to  parallel  pumping  with  the  static  external  field  Hq  parallel  to  the  linearly 
polarized  microwave  field  hQ.  In  one  case,  both  fields  are  directed  normal 
to  the  disk  In  the  other,  the  fields  are  parallel  to  the  flat  surfaces  of  the 
disk  (in-plane).  The  three  bottom  diagrams  correspond  to  perpendicular 
pumping.  For  a  normally  magnetized  disk,  only  an  in-plane  microwave 
field  can  be  used  for  perpendicular  pumping.  For  an  in-plane  magnetized 
disk,  hQ  may  be  normal  to  the  disk  or  in- plane  and  still  be  perpendicular 
to  the  static  field.  A  condensed  notation  will  be  used  in  the  following  dis¬ 
cussion  to  designate  the  specific  configuration  under  consideration.  The 
symbol  (1,1,  II)  will  indicate  Hq1  disk,  hQX  disk,  and  HqII  hQ,  corresponding 
to  the  columnar  listing  under  the  upper  left  diagram.  Other  configurations 
are  designated  in  a  similar  manner. 

Data  for  these  five  configurations  obtained  on  unannealed  1  cm  diam 
disks  are  shown  in  Figs.  42  to  46.  The  threshold  field  h  ^  is  displayed 
as  a  function  of  static  external  field  with  the  disk  thickness  as  a  parameter. 
Experimental  results  were  actually  obtained  for  10-mil  increments  in  thick¬ 
ness.  Data  are  shown  for  20-mil  increments  in  order  to  simplify  the  pre¬ 
sentation.  First  consider  the  results  for  parallel  pumping  shown  in  Figs.  42 
and  43.  The  data  for  normally  magnetized  slabs  shown  in  Fig.  42  evidence 
a  decrease  in  h  ^  with  increasing  field,  a  sharp  minimum  near  25 Oe 
and  a  rapid  increase  at  higher  field.  The  threshold  for  this  configuration 
also  exhibits  a  strong  thickness  dependence  and  hcrit  is  decreased  by  ap¬ 
proximately  1  Oe  for  each  10-mil  increase  in  thickness.  There  is  no  sig¬ 
nificant  change  in  the  field  dependence  of  h  ^  as  the  sample  demagnetizes 
for  Hq  <  1700  0e  The  results  shown  in  Fig.  43  for  the  other  parallel 
pumping  configuration  (11,  II,  II)  corresponding  to  Hq  parallel  to  the  disk 
plane,  h  parallel  to  the  disk  plane,  and  Hq  parallel  to  hQ  are  quite  different. 
For  this  configuration,  the  minimum  h  ^  occurs  for  Hq  =  ilOOOe  and  the 
increase  in  hcrit  at  higher  field  is  not  as  sharp  as  for  Hq  normal  to  the 
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Fig.  41 
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disk.  The  data  exhibit  no  strong  thickness  dependence.  Furthermore, 
the  hcr^  values  increase  sharply  as  Hq  is  reduced  below  2000e  and  the 
samples  demagnetize. 

Consider  the  results  for  the  three  perpendicular  pumping  con¬ 
figurations.  In  Fig.  44  the  instability  threshold  h  ^  as  a  function  of  static 
field  is  shown  for  the  (1,  ll,  1)  configuration,  with  thickness  as  a  parameter. 
Similar  data  for  the  10- mil  disk  were  shown  in  Sec.  C.  The  threshold  ex¬ 
hibits  a  broad  double  minimum  and  increases  sharply  at  high  field.  In 
addition,  hcr^  appears  to  be  a  decreasing  function  of  thickness  for  Hq 
below  the  high- field  minimum  and  an  increasing  function  of  thickness  above 
the  minimum  As  Hq  is  reduced  below  17000e,  demagnetization  results 
in  an  abrupt  decrease  in  hcrit  which  was  discussed  in  conjunction  with 
Fig.  37.  In  Fig.  45  data  are  shown  for  the  (ll,  11,1)  pumping  configuration. 

The  data  exhibit  only  a  broad  minimum  near  15000e  and  a  small  increase 
at  higher  field.  The  threshold  for  this  configuration  is  not  a  strong  function 
of  sample  thickness.  Upon  demagnetization,  furthermore,  there  is  no 
large  change  in  hcj,^  as  found  for  the  (ll,  H,  H)  configuration  (Fig.  43)  For 
perpendicular  pumping  with  an  in-plane  static  field  and  an  out-of-plane 
microwave  field  configuration  (11,1,1)  as  shown  in  Fig.  46,  h  .  has  a 
minimum  near  14000e,  increases  moderately  at  higher  field,  and  exhibits 
a  strong  thickness  dependence.  As  the  sample  demagnetizes,  moreover, 
hcrit  decreases  sharply,  as  shown  by  the  experimental  points  for  H  <  200 Oe. 
For  this  (11,1,1)  configuration  the  change  in  hcrit  upon  demagnetization  is  a 
sensitive  function  of  disk  thickness.  In  Fig.  47,  data  corresponding  to 
Hq  <  200Oe  are  shown  for  different  thickness  slabs  to  demonstrate  the 
nature  of  this  dependence.  The  sharp  decrease  in  ncr^  begins  near  H  -  50Oe 
for  the  10-mil  disk.  For  the  50-mii  thickness,  the  decrease  begins  at  IZOOe 
Due  to  the  increase  in  the  disk  in-piane  demagnetizing  factor  with  thickness, 
thicker  disks  begin  to  demagnetize  at  a  larger  field  than  do  the  thinner 
disks.  Consequently,  the  field  region  over  which  the  decrease  in  h  ^  occurs 
broadens  with  increasing  thickness 


For  the  particular  slab  configuration  shown  in  Figs.  37  and  44 
(1,  II,  i),  annealing  increased  hcrit  significantly,  particularly  in  the  field 
region  corresponding  to  a  partially  magnetized  material.  Data  before  and 
after  anneal  have  also  been  obtained  for  the  other  four  field  configurations. 
The  only  other  configuration  where  annealing  had  a  large  effect  was  per¬ 
pendicular  pumping  with  the  microwave  field  normal  to  the  plane  of  the 
slab,  the  (11,1,1)  configuration.  Threshold  data  for  two  10-mil  disks,  one 
annealed  and  one  unannealed,  are  shown  in  Fig.  48.  The  instability  thresh¬ 
old  decreases  by  about  1  Oe  after  the  1  hr  anneal  at  1200°C.  Thicker  disks 
exhibited  smaller  decreases,  consistent  with  the  earlier  conclusions  that 
annealing  effects  are  related  to  surface  damage,  and  should  decrease  as 
the  surface-to-volume  ratio  is  decreased.  Annealing  did  influence  h  ^ 
for  the  remaining  configurations  but  the  change  was  quite  small,  on  the 
order  of  a  few  tenths  of  an  oersted. 

The  general  characteristics  of  the  results  shown  in  Figs.  42  to  46 
are  summarized  in  Table  I.  The  following  generalizations  can  be  made: 

1)  The  thickness  dependence  is  significant  only  when  the  microwave  field 
is  directed  normal  to  the  plane  of  the  disk,  Z)  The  minimum  h^.^  position 
is  generally  larger  for  normally  magnetized  disks  than  for  in- plane  mag¬ 
netized  disks,  3)  The  increase  in  hcrit  at  high  field  is  sharp  only  for  the 
normally  magnetized  disks  and  h^^  increases  slowly  for  in- plane  mag¬ 
netized  disks,  4)  Large  changes  m  h  upon  demagnetization  occur  only 
for  disks  initially  saturated  parallel  to  the  disk  plane. 

in  comparing  any  of  the  above  results  with  theoretical  considera¬ 
tions,  dielectric  loading  effects  in  the  microwave  cavity  must  be  taken  into 
account.  The  two  configu'ations  tor  which  the  data  exhibit  a  large  thick¬ 
ness  dependence,  the  ix,J.,n)  and  ill,  1,1)  configurations,  the  disks  are 
positioned  in  the  T Ejq^  microwave  cavity  with  the  microwave  cavity  normal 
to  the  disk  plane.  For  these  two  configurations,  portions  of  the  sample  are 
in  a  region  of  strong  e.ectric  field  directed  parallel  to  the  disk  plane  and  it 
is  expected  that  dielectric  loading  might  significantly  influence  the  microwave 
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TABLE  I 


Summary  of  Experimental  Results  for  Parallel  and  Perpendicular  Pumping  in  Thin, 
1  cm  diam  Discs  of  Poly  crystalline  YIG  at  9.2  GHz 


Change  in 

hcrit  uP°n 
Demagnetization 


Nature  of 

Minimum 

Nature  of 

Thickness 

hcrit 

High  Field 

Configuration 

Dependence 

Position 

Increase 

Parallel  Pump 
(1,1,11) 


Perpendicular 

(1,11,1) 


Perpendicular 

(11,1.  H) 


Large 

Decrease 

2500  Oe 

Sharp 

Increase 

* 

No  Change 

Small 

Decrease 

1100  Oe 

Moderate 

Increase 

Large  Increase 

Moderate 

Decrease 

i 

2500  Oe 
(broad) 

Sharp 

Increase 

No  Change 

Negligible 

1500  Oe 
(broad) 

Moderate 

Increase 

' 

No  Change  j 

i 

1 

Large 

Decrease 

1400  Oe 

Moderate 

Increase 

i 

i 

i 

Large  Decrease 

i 

Unanneoled 

Annealed 


Kig.  48  Threshold  Field  hcrit  as  a  Function  of  the  Static  External  Field  H 
for  a  10-Mil  Disk  and  the  (11,1,1)  Configuration.  Data  before  and' 
after  a  1  hr  anneal  at  1200°C  are  shown. 


field  calibration  in  terms  of  the  incident  power  for  the  cavity.  In  order 
to  check  this  possibility,  data  were  obtained  for  all  the  configurations  dis¬ 
cussed  above  with  a  10-mil  YIG  disk  in  contact  with  dielectric  disks  com¬ 
parable  to  YIG  in  dielectric  properties.  It  was  found  that  thick  dielectric 
disks  had  the  same  effect  in  reducing  the  apparent  hcr^t  for  the  (1, 1,  II) 
and  (11,1,1)  configurations  that  thick  magnetic  disks  did.  For  the  (II,  II,  II) 
and  (11,11,1)  configurations,  where  the  electric  field  is  essentially  zero 
over  the  sample,  the  dielectric  disks  had  no  effect  on  the  results  For 
the  (1,  H,l)  configuration,  where  the  electric  field  is  normal  to  the  disk 
plane,  a  small  decrease  was  noted  for  thick  dielectric  slabs.  On  the  basis 
of  these  results,  it  appears  that  dielectric  loading  of  the  cavity  by  the  large 
disk  samples  is  responsible  for  most  of  apparent  hcr^  thickness  depen¬ 
dence  for  the  (1,1,11)  and  (II,  1,1)  configurations. 

Theoretical  curves  for  comparison  with  the  data  were  calculated 
according  to  the  procedure  described  in  Sec.  C  and  are  shown  as  dashed 
lines  in  Figs.  42  to  46  for  the  configurations  discussed  above.  Consider 
the  comparison  for  the  two  parallel  pumping  configurations  (1,1,  H )  and 
(II,  11,1),  in  Figs.  42  and  43.  In  each  figure  only  one  theoretical  curve  is 
plotted.  For  parallel  pumping,  the  coupling  which  leads  to  instability  is 
directly  between  the  microwave  pump  field  and  the  spinwave  which 
goes  unstable.  This  coupling  is  independent  of  disk  thickness.  There¬ 
fore  h  ^  is  not  expected  to  be  thickness  dependent,  except  for  a  small  shift 
of  the  curve  to  higher  field  with  increasing  thickness  because  of  the  reduced 
demagnetizing  factor  normal  to  the  disk.  For  the  (1,1,  II)  configuration  and 
the  data  in  Fig.  42,  the  observed  thickness  dependence  is  due  to  dielectric 
loading  as  discussed  above,  since  the  microwave  field  is  normal  to  the  disk. 
For  the  thinnest  disk,  where  such  loading  effects  are  less  important,  the 
agreement  is  quite  good.  For  the  other  parallel  pump  configuration  (II,  II,  II), 
the  theoretical  curve  is  the  same  as  in  Fig.  42,  except  it  is  shifted  to  lower 
field  by  4irM  to  account  for  the  decrease  in  the  internal  demagnetizing  field. 
At  fields  below  the  minimum  h  .,  position  H  min  the  theoretical  curve 
agrees  quite  well  with  the  data.  Above  this  minimum,  however,  the  experi¬ 
mental  h  ..  does  not  increase  as  rapidly  as  the  theoretical  curve.  The 
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slow  rise  in  h  ..  at  H  >  H  min  is  a  common  feature  of  the  threshold  for 
cnt  o  o 

all  the  in- plane  magnetized  disk  configurations.  It  may  be  explained  as  one 
or  both  of  the  following:  1)  Magnetostatic  mode  resonances  which  occur  at 
fields  above  the  uniform  resonance  position  for  disks  magnetized  in  their 
plane,  2)  Incorrect  estimates  of  AH^  vs  9^  for  k  =  0.  It  is  important  to 
realize  that  this  portion  of  the  curves  is  not  amenable  to  meaningful  theo¬ 
retical  interpretation  because  the  theory  is  not  strictly  applicable  for  k  -  0, 
as  discussed  in  Sec.  C. 


Now  consider  the  disk  results  for  perpendicular  pumping  summarized 
in  Figs.  44  -  46.  In  these  figures,  theoretical  curves  are  also  indicated  by 
dotted  lines.  For  perpendicular  pumping,  however,  h  ^  depends  on  the 
uniform  precession  amplitude  which  is  influenced  by  the  sample  shape.  There¬ 
fore,  the  threshold  field  should  exhibit  a  dependence  on  disk  thickness  which  is 
related  to  the  efficiency  with  which  the  microwave  field  excites  the  uniform 
precession.  First  consider  the  data  for  the  normally  magnetized  disk  con¬ 
figuration  (1,  II,  1)  in  Fig.  44.  Two  theoretical  curves  are  shown.  The  upper 
curve  was  evaluated  for  10-mil  slab  thickness  and  the  lower  curve  for  50-mil 
slab  thickness.  The  results  for  the  10-mil  slab  were  discussed  in  Sec-  C. 

The  theoretical  thickness  can  be  understood  by  considering  the  uniform  pre¬ 
cession  amplitude  dependence  on  sample  shape.  When  the  sample  does  not 
possess  rotational  symmetry  about  the  static  field  direction  as  was  the  case 
for  the  data  in  Sec.  C,  the  uniform  precession  amplitude  is  given  by 


a 
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hQ(u  +  Y) 
XY  -  u2 
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(21) 


if  the  contribution  of  the  anti-Larmor  component  of  the  linearly  polarized 
microwave  field  is  neglected.  The  uniform  precession  amplitude  is  dis¬ 
cussed  in  detail  in  Appendix  I.  In  Eq.  (21), 

X  =  y  (  H  +  (N  -N  )  4rrM  ]  (22a) 

T  O  X  z  s 

and 

Y  =  y  (  H  +  (N  -N  )  4ttM  I  ,  (22b) 

'  o  y  z  s' 
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whe.-e  N  and  N  are  the  two  transverse  demagnetizing  factors,  N  is 
x  y  z 

the  demagnetizing  factor  in  the  static  field  direction,  and  the  microwave 

field  is  in  the  x-direction.  For  the  (1,11,1)  configuration,  Fig,  44,  with 

the  y- direction  in  the  plane  of  the  disk  and  the  x-direction  normal  to  the 

disk,  N  increases  and  N  decreases  with  increasing  thickness  so  that  Y 

and  also  increase.  With  increasing  thickness,  therefore,  the  efficiency 

with  which  the  microwave  field  pumps  the  uniform  precession  increases 

and  the  threshold  field  decreases.  The  two  theoretical  curves  in  Fig.  44 

express  the  above  considerations  quantitatively.  The  shift  of  the  h  ^ 

minimum  to  lower  field  is  due  to  the  reduction  in  the  static  internal  field 

as  decreases.  The  experimental  h  ^  decreases  with  thickness  for 

H  <  H  min  and  the  minimum  h  .,  position  shifts  to  lower  field  with  in- 
o  o  crit  r 

creasing  thickness,  in  agreement  with  theory.  Other  features  of  the  (1,  11,1) 
data  were  discussed  in  Sec.  C. 


The  data  in  Fig. “45  for  the  ( H ,  II,  1)  configuration  (Hq  and  hQ  both 
parallel  to  the  disk  plane)  indicate  that  for  this  configuration,  hCf^  is  an 
increasing  function  of  thickness.  The  theoretical  hcrit  also  increases  with 
thickness.  For  this  configuration,  the  y-direction  is  normal  to  the  disk 
and,  consequently,  N^,  Y,  and  aQ  decrease  with  increasing  thickness,  so 
that  h  ^  is  expected  to  increase.  As  mentioned  above,  the  increase  in 
hcrit  at  static  fields  above  the  minimum  is  not  as  rapid  as  expected  theo¬ 
retically.  Below  the  minimum,  however,  the  agreement  between  the  theory 
and  the  data  is  quite  good.  Finally  consider  the  results  for  the  (11,1,1) 
configuration  shown  in  Fig.  46  where  Hq  is  parallel  to  the  disk  plane  and 
hQ  is  normal  to  the  disk.  For  this  configuration,  the  data  for  the  10-mil 
slab  are  in  good  agreement  with  the  theory  except  for  the  slow  increase  for 
H  >  H  min.  The  data,  however,  exhibit  a  strong  decrease  in  h  ..  with 
increasing  thickness  which  is  related  to  dielectric  loading  since  hQ  is  nor¬ 
mal  to  the  disk  plane.  No  strong  dependence  is  expected  theoretically 


The  y-direction  is  in  the  plane  of  the  disk  (normal  to  Hq 


also  in  the  disk  plane. 


N  and  N  are  equal  and 
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is  independent  of  thickness. 

In  general,  the  agreement  between  the  theory  and  data  for  all 
these  configurations  is  quite  good,  when  the  theory  is  applicable.  As  pointed 
out,  dielectric  loading,  magnetostatic  mode  resonances,  and  the  possible 
errors  introduced  by  applying  the  theory  for  Hq  >  Homin  modify  the  experi¬ 
mental  results  for  some  configurations. 

The  above  discussion  has  been  concerned  exclusively  with  the 
results,  experimental  and  theoretical,  for  saturated  systems.  It  has  been 
shown  that  threshold  field  data  for  a  large  variety  of  situations  can  be  ex¬ 
plained  on  the  basis  of  spinwave  instability  theory.  In  order  to  apply  any 
of  these  results  to  partially  magnetized  systems,  however,  it  is  necessary 
to  establish  a  connection  between  domain  structure  and  the  threshold  fields. 
Such  a  connection  is  the  subject  of  the  remainder  of  this  section.  In  the 
absence  of  any  quantitative  description  of  the  domain  structure,  the  treat¬ 
ment  will  be  primarily  qualitative.  The  discussion  will  be  based  in  part 
on  the  oblique  pumping  presented  in  Appendix  II  and  primarily  on  the  experi¬ 
mental  results  for  thin  disks.  In  that  appendix,  it  is  shown  that  the  data  for 
partially  magnetized  spherical  samples  can  be  explained  if  for  the  trans¬ 
verse  demagnetizing  factor  fTt  it  is  assumed  that  1/3  <  TTt  <  1/2  rather  than 
taking  =  1/3  as  is  the  case  for  saturated  spherical  samples.  On  the 
basis  of  this  observation,  it  may  be  assumed  that  the  demagnetization  is 
accomplished  by  the  formation  of  elongated  domains  which  satisfy  1/3  <  N 
<  1/2  ,  and  that  the  shape  of  these  domains  influence  h  ^  in  much  the 
same  way  that  sample  shape  does  for  saturated  samples. 

Domain  structure  considerations  may  also  be  invoked  to  account 
for  the  threshold  field  behavior  below  saturation  for  the  thin-disk  results 
summarized  in  Figs.  42  to  46.  First  consider  the  results  in  Fig.  42  for 
parallel  pumping  with  Hq  normal  to  the  plane  of  the  disk,  the  (1,1,  II)  con¬ 
figuration.  The  disk  begins  to  demagnetize  for  Hq  «  47tMg  or  approximately 


101 


1750Oe.  At  higher  field,  the  magnetization  is  oriented  normal  to  the  disk 
and  parallel  to  the  microwave  field  by  the  static  field  Hq.  For  Hq  ^  1750  Oe, 
the  internal  field  is  essentially  zero  and  domains  form  as  demagnetization 
processes  act  to  minimize  the  free  energy  of  the  sample.  As  is  evident 
from  Fig.  12,  h£rit  for  the  (1,1,  II)  configuration  does  not  exhibit  any  sharp 
changes  during  demagnetization,  but  increases  gradually  to  approximately 
7  0e  at  HQ  =  0  (for  the  10-mil  disk),  generally  following  the  curve  for  the 
data  above  saturation  This  smooth  dependence  indicates  that  the  magnetiza¬ 
tion  direction  with  respect  to  the  microwave  field  polarization  direction  does 
not  undergo  any  major  reorientation  upon  demagnetization.  It  may  be  con¬ 
cluded,  therefore,  that  in  the  partially  magnetized  or  demagnetized  material 
the  domains  are  aligned  in  directions  which  are  approximately  parallel  to 
the  disk  normal  Such  a  model  has  been  proposed  by  Kittel.*  For  thick 
slabs,  interior  domains  are  expected  to  be  aligned  along  the  normal  direc¬ 
tion  and  small  closure  domains  are  expected  to  form  at  the  sample  surface 
with  magnetization  directions  parallel  to  the  slab  plane.  These  closure 
surface  domains  represent  only  a  small  part  of  the  sample  volume  and  do 
not  appreciably  influence  h  ^  .  Admittedly,  this  model  is  a  gross  over¬ 
simplification  for  polycrystalline  Y1G  material  which  is  cut  into  thin  disks 
by  mechanical  means,  and  which  clearly  has  considerable  surface  dam¬ 
age  and  interior  defects.  Nevertheless,  such  a  model  provides  an 
approximate  working  basis  for  a  consistent  interpretation  of  the 
results. 


For  parallel  pumping  with  Hq  parallel  to  the  film  plane,  the  data 

in  Fig.  43  follow  the  theory  for  saturated  samples  down  to  Hq  -  200  Oe. 

At  lower  fields  the  sample  presumably  demagnetizes  with  interior  domains 

oriented  along  the  disk  normal.  At  Hq  -  0,  therefore,  the  sample  consists 

of  elongated  domains  with  the  microwave  field  perpendicular  to  the  domain 

axis.  In  terms  of  shape,  this  configuration  is  similar  to  perpendicular 

pumping  in  an  axially  magnetized  rod  with  Hq  ~  0.  The  rod  data  in  Fig.  37 

indicates  that  the  h  ..  for  this  situation  should  be  about  12  Oe,  a  value 
crit 

which  agrees  quite  well  with  the  disk  value  of  hcr^  -  11. 8  0e  at  H  =  0 
in  Fig.  43.  Thus,  the  simple  domain  model  assumed  above  is  consistent 
with  both  sets  of  parallel  pumping  data. 
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This  simple  domain  model  is  also  applicable  to  the  perpendicular 

pumping  configurations  for  the  disk  samples.  The  data  for  H  perpendicular 

to  the  disk  plane  in  Fig.  44  indicate  that  hcr^t  increases  to  about  150e  for 

the  saturated  disk  at  H  =  18000e  where  the  internal  field  is  essentially 

zero,  and  then  decreases  to  12  0e  at  H  =0.  The  value  at  H  =  1800  0e  is 

o  o 

consistent  with  the  theory  for  the  saturated  disk.  The  situation  at  Hq  =  0 
is  much  the  same  as  discussed  above.  The  sample  consists  of  elongated 
domains  oriented  along  the  disk  normal,  the  pump  configuration  corresponds 
to  perpendicular  pumping  in  an  axially  magnetized  rod,  and  hcrit  is  approxi¬ 
mately  equal  to  12  0e  (10- mil  disk).  A  comparable  situation  exists  for  per¬ 
pendicular  pumping  where  both  Hq  and  hQ  are  parallel  to  the  disk  plane,  as 
shown  in  Fig.  45.  The  threshold  is  about  1 1  Oe  at  Hq  -  0  for  the  10-mil 
disk.  The  results  for  H  parallel  and  hQ  normal  to  the  disk  plane,  shown 
in  Fig.  46,  indicate  that  hcrit  conforms  to  theory  for  Hq  >  2000e,  and 
decreases  sharply  to  about  8.  50e  at  Hq  =  0  for  the  10-mil  slab.  According 
to  the  model,  Hq  -•  0  corresponds  to  parallel  pumping  with  domains  oriented 

normal  to  the  disk  and  h  .  should  be  equal  to  the  value  at  H  ~  0  for  the 

crit  o 

(1,1,  II)  configurations.  This  value  is  8.  3  0e,  so  that  the  agreement  here 
is  quite  satisfactory. 

The  above  discussion  clearly  demonstrates  the  influence  of  domain 
structure  on  the  threshold  fields  for  partially  magnetized  magnetic  materials. 
A  qualitative  connection  was  evident  from  the  oblique  pumping  results  for 
spherical  samples  in  Appendix  II.  A  more  quantitative  connection  is  evident 
from  the  results  for  parallel  and  perpendicular  pumping  in  disks.  A  simple 
domain  model  for  the  demagnetized  disks  is  found  to  provide  a  consistent 
explanation  of  the  experimental  results  below  saturation. 

E .  Larmor  and  Anti-Larmor  Pumping  in  Rod  Samples 

The  previous  two  sections  have  been  concerned  with  the  influence 
of  sample  shape  and  orientation  on  the  high- power  thresholds  for  parallel 
or  perpendicular  pumping.  For  all  the  data  reported  above,  linearly  polarized 
microwave  excitation  was  used.  Because  of  the  complex  domain  structure 
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associated  with  most  partially  magnetized  systems,  the  actual  pumping 
environment  is  likely  to  be  more  complicated.  In  order  to  obtain  some 
insight  into  the  effect  on  the  threshold  of  changes  in  the  pump  configura¬ 
tion,  two  experiments  have  been  performed:  1)  The  oblique  pumping  study 
of  Appendix  II,  2)  Perpendicular  pumping  with  circularly  polarized  micro- 
wave  fields  rotating  in  a  Larmor  sense,  the  proper  sense  to  excite  ferro¬ 
magnetic  resonance,  and  in  an  anti-Larmor  sense,  the  opposite  sense. 

This  section  is  concerned  with  the  Larmor  and  anti-Larmor  pumping  re¬ 
sults.  Instability  threshold  data  were  obtained  for  axially  magnetized, 
partially  magnetized,  and  demagnetized  10-mil  diam  rods  of  poly  crystalline 
YIG  material  as  a  function  of  external  field,  average  magnetization,  and 
polarization  sense  for  the  circularly  polarized  microwave  field  normal  to 
the  rod  axis.  The  results  were  obtained  at  5.  5  GHz  and  room  temperature. 

The  transition  between  Larmor  and  anti-Larmor  pumping  was  ac¬ 
complished  by  exciting  the  sample  with  one  sense  of  circularly  polarized 
drive  and  varying  the  static  axiai  field  from  a  large  positive  value  to  a 
large  negative  value.  The  effective  sense  of  polarization  changes  when 
the  field  direction  changes,  because  this  changes  the  natural  precession 
sense.  The  situation  is  complicated  by  the  presence  of  magnetic  hysteresis. 
By  traversing  the  hysteresis  loop  in  both  directions,  i.  e.  ,  from  a  large 
positive  field  to  a  large  negative  field  and  then  back  again,  data  can  be 
obtained  for  both  Larmor  and  anti-Larmor  pumping  and  for  the  two  average 
magnetization  values  which  correspond  to  a  given  field.  In  Fig.  49,  the 
threshold  field,  relative  to  the  minimum  value  near  zero  external  field, 
is  shown  as  a  function  of  external  field  from  -20  to  -tlOOOe.  The  data 
show  that  for  Hq  >  0  the  threshold  is  considerably  lower  than  for  Hq  <  0. 

A  Larmor  rotating  pump  field  is  much  more  efficient  in  exciting  the  uni¬ 
form  precession  and  should  lead  to  a  lower  instability  threshold.  The  data 
for  Hq  >  0  and  a  low  threshold  therefore  correspond  to  Larmor  pumping 
and  the  results  for  Hq  <  0  to  anti-Larmor  pumping.  Hysteresis  is  also 
evident  in  the  data.  Values  of  h  .  measured  with  increasing  Hq  are  dif¬ 
ferent  from  those  obtained  with  decreasing  H  .  Additional  data  are  shown 
in  Fig.  50,  with  an  expanded  horizontal  scale  in  order  to  exhibit  more 
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The  Threshold  Field,  Relative  to  the  Value  at  110  =  0  with  11Q  Increasing, 
as  a  Function  of  .Static  External  Field  Circularly  Polarized  Perpendicular 
Pump  Kxeitation.  The  horizontal  scale  is  expanded  from  that  sliown  in 


clearly  the  nature  of  the  hysteresis  in  hcr^.  For  these  data  Hq  was  cycled 
between  +100  and  -30Oe.  For  negative  Hq  corresponding  to  anti-Larmor 
pumping,  the  threshold  field  is  much  larger  for  increasing  HQ  than  for  de¬ 
creasing  H  .  For  positive  Hq  and  Larmor  pumping,  on  the  other  hand,  the 
threshold  is  larger  for  increasing  Hq  than  for  decreasing  Hq.  This  hysteresis 
in  h  is  due  to  the  different  domain  distributions  which  are  possible  for  the 
same  static  field. 

The  applicability  of  the  theory  to  these  pump  configurations  can 
be  determined  by  comparing  the  ratio  of  the  threshold  for  a  large  negative 
static  field  (anti-Larmor  threshold)  to  the  threshold  for  a  large  positive 
field  (Larmor  threshold)  with  the  theoretical  value,  choosing  the  static 
field  values  so  that  the  material  is  magnetized  to  near  saturation.  Theo¬ 
retical  values  were  obtained  from  the  theory  presented  in  Appendix.  I. 

From  Fig.  49,  the  ratio  between  hcrjt  at  -20  and  +1000e  is  equal  to  4.  7. 

The  theoretical  ratio  is  equal  to  14  0  when  measured  values  of  the  average 
4ttM  at  each  field  are  used  in  the  theory.  At  -20Oe,  4irM  is  -1594G  and 
at  +1000e,  4rrM  is  near  the  saturation  value  of  1750G.  The  agreement  is 
poor.  The  large  discrepancy  which  does  exist  is  probably  due  to  the  non- 
saturated  condition  of  the  rod  at  Hq  -  -20Oe.  The  partial  misalignment 
of  domains  gives  rise  to  stray  fields,  which  lead  to  Larmor  coupling  even 
for  a  microwave  field  rotating  in  an  anti-Larmor  sense.  The  -20Oe  value 
was  the  largest  attainable  static  field  in  the  reverse  direction  for  the  electro¬ 
magnet.  Data  at  higher  static  fields  should  result  in  improved  agreement. 
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IV.  SUMMARY  AND  CONCLUSIONS 


Measurements  of  y\  k' ,  y  y",  k",  and  y"  have  been  made  on 

z  z 

several  materials.  We  have  found  that  Rado's  theory  on  the  permeability 
of  partially  magnetized  ferrites  gives  approximately  the  correct  value  for  x\ 
Some  deviations  were  observed,  however,  the  most  significant  of  these  being 
that  for  some  materials  «'  is  not  a  unique  function  of  4ttM  but  depends  on  the 
hysteresis  history  of  the  sample.  Rado's  theory  does  not  give  the  correct 
results  for  y'  and  y'.  We  have  shown  how  the  rf  demagnetizing  fields  from 
poles  on  grain  and  domain  boundaries  produce  terms  which  Rado  ignored 
in  his  theory  and  which  contribute  primarily  to  y’  and  y^  .  Empirical  for¬ 
mulae  for  y*  and  y^  which  fit  our  data  have  been  presented.  The  results 
of  the  measurement  of  y'\  x",  and  y’  vs  4n-M  for  various  values  of  u,T/u 
have  been  presented  in  the  figures.  At  values  of  close  to  1. 0  their 

dependence  on  4trM  is  substantial.  At  values  of  u^/uj  ^  0.  8  this  dependence 
is  small.  Furthermore,  x"  is  then  much  less  than  y".  In  this  region  we 
feel  a  single  loss  parameter,  y"  at  4nM  =  0,  is  adequate  to  characterize 
the  material.  The  results  of  measurements  of  this  single  parameter  on 
several  very  low- loss  materials  at  various  frequencies  in  this  region 
were  presented. 


The  investigation  of  high- power  phenomena  has  accomplished  the 
following  objectives:  1)  A  general  theory  has  been  developed  for  first- 
order  instability  threshold  fields  in  ellipsoidal  samples  magnetized  to 
saturation  along  a  principal  axis  and  with  an  arbitrary  microwave  field 
configuration,  2)  The  effect  of  sample  shape  and  pump  configuration  on 
threshold  for  saturated  systems  has  been  characterized  and  explained  in 
terms  of  the  theory,  3)  A  connection  has  been  observed  between  the 
domain  structure  and  the  threshold  fields  for  partially  magnetized  systems 
which  indicates  that  domain  shape  plays  the  same  role  in  determining  the 
threshold  for  multidomain  systems  that  sample  shape  plays  for  saturated 
systems. 


The  experimental  data  were  obtained  by  using  a  semiautomated 
instrumentation  system  which  enables  h  ..  determinations  to  be  made 
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accurately  and  quickly,  thereby  facilitating  the  accumulation  of  a  large  amount 
of  data  as  a  function  of  sample  shape,  pump  configuration,  static  field,  and 
magnetization  state.  Perpendicular  pump  data  for  axially  magnetized  rods, 
spheres,  and  normally  magnetized  slabs  indicate  that  gross  changes  in 
sample  shape  have  a  large  influence  on  h  ^  and  that  annealing  also  changes 
hcrit  for  samples  with  large  surface-area-to-volume  ratios.  For  the  nor¬ 
mally  magnetized  disks,  the  presence  of  magnetostatic  mode  resonances 
in  the  resonance  spectra  modifies  the  threshold  field  from  the  theoretical 
value.  Apart  from  modifications  due  to  annealing  and  magnetostatic  mode 
phenomena,  the  data  are  consistent  with  theory  and  indicate  that  the  theory 
is  applicable  for  a  variety  of  sample  shapes.  Parallel  and  perpendicular 
pump  data  for  disk-shaped  samples  have  been  obtained  which  are  also  con¬ 
sistent  with  theory.  For  configurations  with  the  microwave  field  normal 
to  the  disk,  dielectric  loading  in  the  microwave  cavity  must  be  taken  into 
account.  Additional  results  for  perpendicular  pumping  in  rod-shaped  samples 
with  circularly  polarized  excitation  indicate  that  the  threshold  is  much  higher 
for  anti-Larmor  than  for  Larmor  pumping,  and  the  increase  is  slightly  less 
than  expected  theoretically. 

The  results  for  partially  magnetized  samples  in  the  configurations 
summarized  above  indicate  that  the  theory  is  applicable  to  these  multi- 
domain  systems,  if  domain  shape  is  taken  into  account.  The  data  for  spheri¬ 
cal  samples  are  consistent  with  theory,  if  the  samples  are  assumed  to 
demagnetize  by  the  formation  of  elongated  domains  along  the  static-field 
direction.  Similarly,  the  disk  data  can  be  explained  if  a  domain  structure 
consisting  of  rod-shaped  domains  normal  to  the  disk  plane  are  assumed. 
Although  such  domain  models  represent  a  gross  oversimplification  of  the 
demagnetization  processes  involved,  they  are  qualitatively  correct  The 
consistent  correlations  between  the  data  and  the  theory,  which  are  made 
possible  by  these  models,  provide  convincing  evidence  that  domain  shape 
for  partially  magnetized  samples  and  sample  shape  for  saturated  systems 
are  directly  applicable  to  partially  magnetized  materials,  if  reasonable 
assumptions  concerning  the  domain  structure  can  be  made.  On  this  basis, 
predictions  of  peak  power  capability  for  materials  and  geometries  relevant 
to  remanent  devices  should  be  possible. 
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APPENDIX  I 


THEORY  FOR 

THE  FIRST-ORDER  SPIN- WAVE  INSTABILITY  THRESHOLD 
IN  FERROMAGNETIC  INSULATORS  OF  ELLIPSOIDAL  SHAPE 
WITH  AN  ARBITRARY  PUMPING  CONFIGURATION’" 

Carl  E.  Patton 
Raytheon  Research  Division 
Waltham,  Massachusetts  02154 


ABSTRACT 


The  spin-wave  instability  theory  of  nonlinear  effects  in  ferromagnetic 
insulators  developed  by  Suhl  for  perpendicular  pumping  and  by  SchlSmann  for 
parallel  pumping  has  been  extended  to  include  pumping  with  an  arbitrary 
microwave  configuration.  The  theory  is  applicable  to  samples  of  ellipsoidal 
shape,  magnetized  to  saturation  along  one  principal  axis,  but  not  necessarily 
an  axis  of  rotational  symmetry  as  in  the  Suhl  theory.  The  theory,  in  its 
present  form,  is  qui;e  versatile  and  applicable  to  a  large  variety  of  situa¬ 
tions  which  are  not  accessible  on  the  basis  of  the  earlier  calculations. 
Equations  are  developed  specifically  for  first-order  processes.  Several 
examples  are  considered  to  demonstrate  the  versatility  of  the  theory: 

1)  Perpendicular  pumping  with  a  linearly  polarized  microwave  field; 

2)  Oblique  pumping  with  a  linearly  polarized  field  at  an  angle  between  zero 
and  90°  to  the  static  field;  3)  Perpendicular  pumping  with  a  circularly 
polarized  field  rotating  in  the  anti-Larmor  sense. 


INTRODUCTION 


Although  the  spin-wave  instability  theory  of  nonlinear  effects  in 
ferromagnetic  insulators  developed  by  Suhl*  and  SchlSmann  et  al.^  is  quite 
general  in  concept,  the  specific  calculations  were  done  only  for  restricted 
pumping  configurations  and  sample  shapes.  Pumping  configurations  were 
restricted  to  linearly  polarized  microwave  fields  along  the  static  field 
direction  (parallel  pumping)  or  circularly  polarized  fields  perpendicular 
to  the  static  field  direction  with  the  proper  sense  of  rotation  to  excite 
ferromagnetic  resonance  (perpendicular  pumping).  The  results  are  appli¬ 
cable  only  to  samples  of  ellipsoidal  shape  magnetized  to  saturation  along  a 
principal  axis.  For  perpendicular  pumping  the  saturation  direction  must 
also  be  an  axis  of  rotational  symmetry.  The  theory  has  been  extended  to 
include:  1)  Samples  of  ellipsoidal  shape  magnetized  to  saturation  by  an 
applied  static  field  HQ  along  a  principal  axis,  not  necessarily  an  axis  of 
rotational  symmetry;  2)  An  arbitrary  configuration  for  the  microwave 
field.  A  general  expression  for  the  first-order  instability  threshold  field 
as  a  function  of  spin-wave  propagation  direction,  static  field,  sample  shape, 
and  pumping  configuration  is  obtained.  Results  are  presented  for  several 
interesting  special  cases,  primarily  to  illustrate  the  versatility  of  the  theory 
in  evaluating  the  instability  thresholds  for  given  situations.  Three  cases 
are  considered:  1)  Perpendicular  pumping  with  a  linearly  polarized  micro- 
wave  excitation;  2)  Perpendicular  pumping  with  a  circularly  polarized 
microwave  field  rotating  in  the  anti-Larmor  sense  (opposite  to  that  required 
for  ferromagnetic  resonance),  and  3)  Pumping  with  a  linearly  polarized 
microwave  field  at  an  arbitrary  angle  to  the  saturation  direction.  Similar 
results  for  these  special  cases  and  an  infinite  magnetic  medium  have  been 
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reported  by  Kitaev,  et  al.  The  special  case  of  oblique  pumping  has  also 
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been  examined  by  Stern  and  Green  et  al.  The  present  specific  results, 
however,  are  obtained  from  a  general  extension  of  the  original  SchlSmann- 
Suhl  theory  and  are  easily  applicable  to  more  complicated  instability  calcu¬ 
lations. 
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SPIN- WAVE  EQUATION  OF  MOTION 


The  equation  of  motion  for  the  spin-wave  amplitude  a,  ,  first 

1  6  K 
obtained  by  Suhl  and  discussed  by  SchlSmann,  forms  the  basis  of  the 

present  extension  of  the  theory  to  the  general  case.  The  general 

approach  and  much  of  the  mathematics  closely  follow  the  discussion  in 

Ref.  6.  Consider  the  undamped  equation  of  motion  for  the  magnetization 

M(r)  in  terms  of  a  unit  vector  a^r)  *  M(r)/M: 

do(r)/dt  *  -  y  a(r)  XH(r)  ,  (1) 

where  y  is  the  gyromagnetic  ratio,  assumed  positive  for  electrons,  4-rrM 
is  the  saturation  induction,  and  H(r)  is  the  total  internal  magnetic  field. 
This  internal  field  consists  of  four  parts,  the  static  external  field  Hq,  the 
external  microwave  field  hQ,  an  effective  field  due  to  exchange  Hg  (r),  and 
an  effective  field  due  to  dipolar  interactions  (r).  All  fields  except  H(r) 
are  assumed  to  vary  as  elut  where  u  is  the  angular  excitation  frequency. 
By  expanding  the  a(r)  in  a  Fourier  series 

.a(r)  =  .  (2) 

k 

assuming  that  Hq  and  hQ  are  uniform  over  the  sample  volume,  selecting 
the  coordinate  axes  to  be  along  the  principal  axes  of  the  saturated  ellip¬ 
soidal  sample  with  HQ  in  the  z -direction,  and  assuming  that  Mz  »  M  is 
satisfied,  the  equation  of  motion  for  a.^  can  be  obtained.  The  vector  Jt. 
is  the  propagation  constant  or  wave  vector  for  the  spin  wave,  and  the^ 
are  spin-wave  amplitudes  defined  according  to  Eq.  (2).  Following  Schlb- 
mann,^1  the  external  fields  may  be  expressed  in  terms  of  internal  fields. 


(3a) 

(3b) 
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where  4irM  is  the  saturation  induction  and  N  are  the  sample  demag- 

y,  z 

netizing  factors.  The  result  can  be  expressed  conveniently  in  terms  of 
complex  variables 


a,  »  a*  +  i 
k  k  k 


and 


kT  *  kx  +  1  k, 


h  *  h„  +  i  h 
x  y 


(4a) 

(4b) 


The  equation  of  motion  is 


a.  *iY  {  ( H  +  h*  +  4ir M  N  (1  -  a*)) 


+  D  £  (k')2  (ak, 

k* 


o  z  o 

ak-k'  "  V  ak-k»> 


ak  -akh 


with 


+  4»m  £  *•  ■«*>  (k-x  <_k,  -  vk.>/<£.’>2} 


k'*0 


Z  X 

a,  =  o.  T 
k  ko  2 


i  1  V  V 

k' 


k' 


(5) 


(6) 


where  D  is  a  phenomenological  exchange  parameter  defined  by 
Hex(r)  *  DV2a(r)  and  <5ko  is  the  Kronecker  delta  function.  Equation  (5) 
with  h^  =  0  or  h  =  0  provides  the  basis  for  the  calculation  of  perpendicular 
or  parallel  pumped  instabilities,  respectively.  The  inclusion  of  both  h 
and  hQ  in  the  equation  of  motion  for  provides  the  basis  for  the  present 
extension  of  the  theory  to  include  an  arbitrary  pumping  configuration. 

In  a  conventional  low  power  microwave  experiment,  only  the  uniform 
mode  is  excited  to  a  large  amplitude  and  the  a^,  other  than  for  k^=  0,  are 
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negligibly  excited.  In  calculating  the  instability  threshold,  therefore, 
only  terms  to  first-order  in  the  withjc#  0  need  be  considered.  Including 
only  terms  linear  in  the  o^,  the  equation  of  motion  for  with  Jt  *  0  becomes 

2i4>k  *  - 

ak  s  1  [  (Ak  +  Ck}  ak  +  (Bke  +  Dk}  a-kJ  +  i  yho  ak  '  (7) 

where 

Ak  =  yh  +  yDk2  +  jUjjj  sin"  9 k  ,  (8a) 

Bk=I"m3in2’k  (Sb) 

1  '^k  * 

Ck  3  'Iwm  003  9k  ain9k  (e  ao  +  e  ao) 

+  ?wm  I2Nz  +  (3  c°s2dk  -  ni^ol2  +!Vha*  ,  (8c) 

and 

Dk  3  I  ^  “»  -  “m  003  9k  3in  *k  e  ^  *o 

+  ?  “m  (2  cos2  9k  »o  ‘  3i"2  9k  e  k  l“o  I2’  +  ?  vh  %  ’  <8d) 


in  which  um  is  equal  to  y4irM,  9 k  is  the  angle  between  the  spin-wave  wave 
vector  and  the  saturation  direction,  and  4>k  is  the  angle  between  the  pro¬ 
jection  of  jc.  on  the  x-y  plane  perpendicular  to  the  saturation  direction  and 
the  x-axis.  The  and  terms  are  linear  in  the  including  k_=  0  terms, 
the  and  terms  are  linear  in  the  for  Jc,*  0  terms  only  and  include 

a  terms  to  first  or  second  order,  and  the  hz  term  results  from  the  micro- 
o  o 

wave  field  component  along  the  direction  of  the  static  external  field  Hq.  The 

uniform  precession  amplitude  a  is  nonzero  only  for  a  nonzero  microwave 

^  xv 

field  component  perpendicular  to  HQ.  If  hQ  and  h£  are  zero,  aQ  is  zero,  and  the 

C,  and  D,  terms  vanish.  Equation  (7)  then  leads  to  the  instability  threshold 

for  parallel  pumping  calculated  by  SchlSmann.  ,0If  hz  is  zero  Eq.  (7)  is  the 

°  1 

same  as  obtained  in  Suhl's  perpendicular  pumping  theory. 


t 
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In  the  absence  of  any  microwave  excitation,  only  the  and  B  are 

nonzero  and  the  equation  of  motion  can  be  diagonalized  by  a  linear  Holstein - 
Primakoff  transformation 


(9) 

(10a) 

(10b) 

(11a) 

(lib) 
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Gk  *  Dk  Xk  +  Dk  ^k  '  <Ck  +  Ck )  \  ^k 

e24k  -2i\  *  2i<j>,  * 

*  V-  [D*(A*+Ve  +Dk(Ak-yk)e  -<Ck  +  Ck^kJ  * 

(14) 

The  equation  of  motion  given  by  Eq.  (12)  contains  a  contribution  due  to  the 
transverse  components  of  h  in  the  F.  and  G.  terms,  and  a  contribution 
due  to  the  parallel  component  of  in  the  h^  terms. 

GENERAL  FIRST-ORDER  INSTABILITY  THRESHOLD 


On  the  basis  of  Eqs.  (12)  and  (14),  the  threshold  fields  for  first  - 
order  and  second-order  instabilities  can  be  calculated  in  a  straightforward 
manner  as  described  in  Ref.  6.  The  F^  +  yh^  term  in  Eq.  (12) 

leads  only  to  a  slight  modulation  and  shift  of  the  spin -wave  frequency.  The 
G^  -  Yh^  exp  (2i<f>k)/cjk  coupling  term  leads  to  an  unstable  growth  of  the 
b^  with  time  for  sufficiently  high  microwave  power  levels.  For  convenience, 
let 


B, 


(15) 


consider  an  equation  of  motion  of  the  form 

bk  =  1  <ukbk  +  Gk  b*k*  ’  06) 

and  introduce  energy  dissipation  by  replacing  w,  by  fl,  =  u,  +  irj,  where  rj, 

K  K  X  K  K 

is  the  relaxation  rate  for  the  b^.  This  relaxation  rate  is  related  to  the  spin- 
wave  linewidth  AHk  by  AHk  =  2^/y.  The  Gk  can  be  expressed  as 

"f2  Gk(ni  einut 

n  =  -  2 


(17) 


and  the  first -order  instability  can  be  evaluated  by  considering  only  the 
n  *  1  term  of  Eq.  (17),  that  part  of  2^  which  varies  as  e^*.  Similarly, 
the  second  order  instability  is  obtained  by  considering  the  n  =  2  term 
although  thi3  case  will  not  be  examined  here.  Considering  the  n  =  1  term 
only,  the  equation  of  motion  becomes 

bk  =  bk  +^k1}  eUJt  b*-k^  *  (18) 

It  is  convenient  to  rewrite  this  equation  in  terms  of  a  new  variable  bko(t) 
defined  by 

bk  3  bko«>  '“‘/2  •  <>” 

and  examine  the  conditions  under  which  a  stationary  solution  for  the  bkQ  is 
permitted.  By  combining  the  differential  equations  for  b^Q  and  b  which 
result  from  Eq.  (18)  and  (19),  one  obtains  a  second  order  differential  equa¬ 
tion, 


which  has  a  solution  b^Q  oc  exp  (xt)  where 

*  =9k±l  lGk1)|2  '  K  '  W/Z)2]'  •  (21) 

I  1  )  I 

The  stationary  solution  corresponds  to  *  ■  0.  The  smallest  value  of  IG^  | 
for  which  this  solution  is  permitted  occurs  when  the  spin-wave  frequency 
is  equal  to  one-half  the  excitation  frequency  w.  The  quantity  is 

proportional  to  the  microwave  field  amplitude  hQ  and  the  condition  {2^  |  = 

which  corresponds  to  *  *  0  for  =  u/2  determines  the  instability  threshold 
amplitude  for  the  state  (w/2,  k).  Larger  amplitudes  re  suit  in  a  positive  k  and  an 
exponential  increase  in  b,  with  time. 


The  problem  therefore  reduces  to  the  extraction  o(7j^  from 
Eqa.  (8),  (14).  (15).  and  ()7).  In  order  to  obtain  2^  it  is  necessary  to 
consider  the  explicit  time  dependence  of  hQ  and  the  resultant  uniform 
precession  amplitude  aQ.  Let  hQ  be  given  by 


and 


hQ  ■  hQ  ax  cos  (wt  +  «x) 


ho  *  ho  ay  C0S  (wt  +  V 


h^  »  h_  a„  cos  (ut) 
o  o  z 


(22a) 

(22b) 

(22c) 


Equations  (22a-22c)  provide  an  adequate  basis  for  describing  most  exper¬ 
imentally  attainable  drive  configurations.  For  example,  setting  fi  and  6  „ 

222  x  y 

to  zero  and  requiring  that  a  +  a  +  a  equal  unity  corresponds  to  linearly 

x  y  z  , 

polarized  excitation  making  an  angle  *  -  cos  a  with  respect  to  the 

-1  2 

z -direction  and  an  angle  *  *  cos  ax  with  respect  to  the  x-direction.  With 
the  abbreviations 


(23a) 

(23b) 

the  complex  uniform  precession  amplitude  aQ  can  be  calculated  from 
Eqs.  (1).  (2).  (4).  (22).  and  (23)  to  be 


X  *  y(H  +  4wM  Nx) 

and 

Y  =*  y(H  +  4wM  N  ) 


Yh0  iijt  -iut. 

ao  =  ~1~  (qL  *  +  qA  e  } 


where 


(«+ Y)  a  e  x  +  i(u+X)  a  e  . 3 
1L - — ~ 1 - * - 


XY  -  u 


(24) 


qA  =  - 


(w-Y)  ax  e 


-i<5  -i<5 

x  +  i(u-X)  a  e  y 
- 2 - Z - 


XY  -  u 


and 


Damping  of  the  uniform  precession  may  be  taken  into  account  by 
replacing  X  and  Y  by  flQX  »  X  +  ir;o  and  »  Y  +  i  t]q,  respectively. 

The  relaxation  rate  tjq  is  related  to  the  uniform  precession  resonance 
linewidth  AH  by  AH  »  Such  a  procedure  is  reasonable  when  only 

a  Larmor  rotating  pump  field  is  considered,  as  in  the  Suhl  theory.  For 
anti-Larmor  pumping,  the  introduction  of  damping  in  this  way  results  in 
a  physically  unreasonable  negative  energy  absorption  by  the  magnetic 
system.  The  origin  of  this  inconsistency  rests  in  the  Bloch-Bloemberegen 
formulation  of  relaxation  for  the  uniform  precession^  which  is  equivalent 
to  the  above  procedure.  In  this  formulation,  the  magnetization  relaxes 
toward  the  direction  of  the  static  field,  not  the  instantaneous  internal 

field,  and  results  in  a  negative  energy  absorption  for  anti-Larmor  pumping. 
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If  damping  is  taken  into  account  by  means  of  a  Landau-Lifshitz  or  Gilbert 
form  for  the  relaxation,  the  inconsistency  is  eliminated.  For  sufficiently 
narrow  linewidth  materials,  such  fine  points  do  not  significantly  alter  the 
results.  For  first-order  processes  considered  here,  damping  may  be 
neglected  completely  for  narrow  linewidth  materials  without  significantly 
modifying  the  results.  For  second -order  processes,  the  above  procedure 
is  adequate  for  Larmor  pumping  and  damping  can  be  neglected  in  considering 
anti-Larmor  pumping. 

The  contributions  to  Gj  are  determined  only  by  the  explicit  time 
z  K  2 

dependences  of  h  and  G,  .  From  Eqs.  (15)  and  (22),  the  h  contribution 
O  K  ^  O 

is  simply  -  (wm/w)  (Yhy/2)  sin  0^  az  exP  (21$^).  The  G^  contribution  is 
determined  by  Eqs.  (8),  (14),  (22),  and  (24).  The  and  B^  have  no  explicit 
time  dependence.  The  only  terms  in  C^,  D^,  and  their  complex  conjugates 
which  vary  as  elu,t  are  those  which  cortain  aQ  or  aQ  to  first  order  and  do  not 
contain  hQ  explicitly.  After  collecting  all  terms  in  G^  which  vary  as 
the  result  is 


Pf(l)  m 
k  “ 


2i4>k  r  (  , 

e  sin  0^  cos  0^  /  q^  [1(0^)  +  w/2  ]  e 
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In  the  above 
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and  f($k)  is  equal  to  yH  +  yDk2  in  the  dispersion  relation  of  Eq.  (11c).  The 
instability  threshold  field  amplitude  for  spin -wave  state  (uk  *  w/2,  i?^,^) 
is  simply  the  value  of  hQ  for  which  |G^i  3  yAHk/2  is  satisfied.  Let  w  be 
defined  by 


and 


u_  yh 
m  o 

—  ~r 


(27) 


-i*k  *  ,  %  2 


W  *  sin«k  cosdk{qL(f(«k)+u/21  e  K +  qA  (f(«k)  -  u/2J  e  K+sin  akaz} 

(28) 


The  threshold  hc  is  given  by 


h 


c 


wm  j  W  J 


(29) 


First-order  processes  are  restricted  to  external  field  values  for  which  spin 
waves  with  uk  »  u/2  are  allowed.  From  Eq.  (11c),  this  condition  is  equivalent 
to  the  requirement  u/Z  >  yH.  As  yH  approaches  w/2,  0k  approaches  zero 
and  the  sin  9y  factor  in  W  causes  hc  to  diverge. 

The  experimentally  observed  threshold  for  a  first-order  instability 
should  correspond  to  the  minimum  value  of  hc  with  respect  to  k_  for  a  given 
sample  geometry,  microwave  field  configuration,  and  applied  static  field. 

This  theoretical  minimum  can  be  obtained  by  minimizing  jwj  1  with  respect 
to  #k  and  4>k  subject  to  the  condition  that  the  value  of  k  at  the  minimum  be 
real.  This  condition  is  equivalent  to  the  requirement 


t(9k)  -  yH  2  0 


(30) 


i 

I 

i 


li 


Let  and  be  the  values  of  9^  and  ^  which  minimize  Jw|  *. 

If  does  not  satisfy  the  inequality  of  Eq.  (30),  the  which 

actually  corresponds  to  the  instability  threshold  is  given  by  the  0^  solution 
to  f(0k)  -  yH  3  0  and  Jj.  is  zero. 


^crit 


u  _ 

um  }W(0 


*Hk 

(min) 

k 


.  (min/ 
*k 


I 


(31) 


is  the  expected  first-order  instability  threshold  field  for  a  given  sample 
geometry,  microwave  field  configuration,  and  static  field. 


In  the  above  analysis,  a  first-order  spin-wave  instability  threshold 
expression  has  been  obtained  for  1)  samples  of  ellipsoidal  shape  magnetized 
to  saturation  along  a  major  axis  and  2)  an  arbitrary  microwave  pumping  field 
configuration.  No  assumptions  were  made  concerning  sample  symmetry, 
other  than  the  ellipsoidal  shape,  or  microwave  field  geometry,  except  that 
all  fields  were  assumed  uniform  over  the  sample  volume.  Insofar  as  most 
of  the  analysis  is  quite  similar  to  that  previously  reported  by  Suhl1  and 
SchlSmann,^  it  is  important  to  briefly  summarize  the  specific  differences 
between  this  and  the  earlier  calculations  which  make  possible  the  extension 
of  the  theory  to  more  complicated  situations.  The  first  difference  is  the 

Z  X 

inclusion  of  both  a  parallel  component  hQ  and  transverse  components,  hQ 
and  h^,  of  the  external  microwave  field  in  the  initial  equation  of  motion 
given  by  Eq.  (5).  By  considering  only  terms  to  first  order  in  the  due 
to  both  parallel  and  transverse  field  components  in  the  equation  of  motion 
and  applying  a  linear  Holstein-Primakoff  transformation  to  these  terms, 
the  coupling  coefficient  in  the  equation  of  motion  which  leads  to  instability, 
expressed  in  of  Eq.  (16),  is  obtained  for  a  general  microwave  field 
configuration.  The  second  difference  rests  in  the  expression  for  the  uni¬ 
form  precession  amplitude  aQ  which  was  used  to  evaluate  the  coupling 
coefficient  associated  with  the  transverse  components  of  h  .  By  utilizing 
the  general  expression  for  aQ  given  by  Eq.  (24),  which  is  valid  for  an 
arbitrary  transverse  microwave  field  polarization  and  different  transverse 
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demagnetizing  factors,  the  requirements  on  sample  symmetry  and  the 
transverse  pump  polarization  are  eliminated.  These  two  changes  con¬ 
cerning  Ji  and  aQ  constitute  the  major  differences  between  the  present 
calculation  and  the  earlier  theories.  For  simplicity,  the  discussion  has 
been  limited  to  first-order  instability  processes.  Second-order  pro¬ 
cesses  and  the  3teady  state  behavior  above  threshold  have  not  been  con¬ 
sidered  explicitly.  From  the  general  expression  for  of  Eq.  (15),  the 
threshold  for  second-order  processes  with  *  u  and  n  *  2  can  be  easily 
evaluated  in  a  manner  analogous  to  that  presented  here.  From  the  general 
equation  of  motion  of  Eq.  (5),  moreover,  steady  state  behavior  can/be 
calculated  by  including  higher  order  terms  in  the  working  equation  of 
motion.  The  theory  is  quite  general  and  applicable  to  a  large  variety  of 
situations.  Several  simple  applications  of  the  theory  for  first-order 
processes  will  be  considered  in  the  next  section. 


DISCUSSION  FOR  SPECIAL  CASES 

The  above  refinements  in  the  theory  result  in  several  modifications 
cf  the  earlier  predictions  and  extend  the  applicability  of  the  theory  to  more 
general  geometries  and  microwave  field  configurations.  Although  the  result 
expressed  by  Eqs.  (24),  (26),  (28),  and  (31)  are  quite  general,  the  minimi¬ 
zation  procedure  is  most  easily  accomplished  using  numerical  techniques. 

It  is  more  illustrative  to  consider  several  special  cases  for  which  the  equa¬ 
tions  simplify  considerably,  but  which  still  demonstrate  the  new  or  modified 
results  of  the  extended  theory.  Three  specific  situations  will  be  considered. 

1)  Perpendicular  pumping  with  linearly  polarized  microwave  excitation; 

2)  Perpendicular  pumping  with  circularly  polarized  drive  rotating  in  the 
anti-Larmor  sense;  3)  Oblique  pumping  with  linearly  polarized  drive.  In 
case  (1),  the  contribution  of  the  anti-Larmor  rotating  component  of  :he 
linearly  polarized  pump  field  to  h  will  be  demonstrated  and  the  effect 
of  sample  geometry  on  the  threshold  will  also  be  considered.  In  case  (2), 
the  threshold  for  anti-Larmor  pumping  will  be  calculated  and  compared 
with  h  for  Larmor  pumping.  In  casa  (3),  the  dependence  of  b  ..  on  the 
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angle  between  the  polarization  axis  for  a  linear  polarized  pump  field  and 
the  static  field  direction  will  be  demonstrated.  In  all  three  cases,  the 
tjq  contributions  to  and  q^  will  be  neglected  for  simplicity.  As  long 
as  the  difference  between  the  applied  field  HQ  and  the  field  for  resonance 
at  the  excitation  frequency  u  is  large  compared  to  the  uniform  precession 
resonance  linewidth,  the  tjq  terms  do  not  contribute  significantly  to  the 
uniform  precession  amplitude  aQ.  For  first-order  instability  processes, 
the  internal  field  must  be  less  than  u/2y  and  the  above  condition  is  usually 
satisfied.  In  all  cases,  the  dependence  of  AHk  on  k  will  be  neglected  in 
the  minimization  procedure. 

Case  (1).  For  perpendicular  pumping  with  a  linearly  polarized 

microwave  field  in  the  x-direction,  6  ,  a  ,  and  a  in  Eq.  (22)  are  zero, 

x  y  z 

a^  is  equal  to  unity,  and  the  uniform  precession  amplitude  is  given  by 
Eq.  (24)  with 


and 


Y  +  u 
XY  -  J" 


qA 


_  Y  -u 

*  - j- 

XYu 


(32a) 


(32b) 


With  qL  and  q^  from  Eqs.  (32a,  b),  W  in  Eq.  (28)  simplifies  to 


sin«.  cos  0,  -i$,  i<j> 

W» -  - —  {  (Y +w)  [f(d,  ) +w/2]  e  K  +  (Y-w)[f(d,  )-«/£]  e  K} 

u  -  XY  K  * 

-  (33) 


The  $k  spin-wave  angle  determines  the  relative  phases  of  the  two  complex 
terms  in  the  curly  brackets  of  Eq.  (33).  Therefore,  the  magnitude  of  W 
is  maximized  and  the  threshold  hc  is  minimized  for  4>k  *  such  that 

the  phases  of  these  two  terms  are  equal. 


|W(9k,  *^"1)1 


sindk  cos dk 
|uZ  -XYj 


{(u  +  Y)[w/2+f(dk)J  +  |w-Y|  [u/2-f(dk)l} 

(34) 
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The  minimization  of  |W(0k,  *  and  hc  with  respect  to  0^  determines 

the  threshold  as  discussed  above  and  h  . .  is  given  by 


^crit 


u> 


AHk  |u 
(min) 


sin  $ 


2  -  XYj 
cos  9^ 


(35) 


{  (w  +  Y)(u/2  -  f(0jjmin)))  +  |u  -  Y|  I u/2  -  f{0jjm:n,)]>  **  . 

The  value  for  is  a  complicated  function  of  the  sample  shape,  frequency, 

magnetization,  and  applied  field.  The  value  of  $k  \  however,  is  determined 
only  by  the  relative  signs  of  the  real  terms  in  Eq.  (33).  With  f(0k)  as  defined 
by  Eq.  (26),  the  expression  [ f -  u/2]  is  always  nega¬ 
tive.  If  (Y  -  w)  is  also  negative,  Eq.  (34)  corresponds  to  $kmin^  =  0  or 
4>kmin^  a  ir.  If  (Y  -  u)  is  positive,  is  +  ir/2.  For  first-order  processes, 

the  static  internal  field  H  is  constrained  by  the  requirement  H<u/2y.  The 
static  internal  field  value  below  which  <|>km*n)  is  0  or  rr  and  above  which 
<f>kmin)  is  +  it/ 2  is  determined  by 

Y  -  u  =  -yH  +  -y4irMNy,  -  u  -  0  .  (36) 


This  transition  field  is  less  than  u/2y  only  if  the  N  r  demagnetizing  factor 
is  greater  than  For  yttrium  iron  garnet  at  °  ?  GHz,  ♦his  condition 

reduces  to  N  >  0.  94.  The  instability  can  occur  for  =  +  w/2  only  for 

very  thin  slabs  magnetized  in  their  plane  ard  with  the  linearly  polarized 
field  also  in  the  plane  of  the  slab,  perpendicular  to  the  saturation  direction. 
For  most  situations,  therefore,  instability  occurs  for  spin  waves  with  $k 
equal  to  zero  or  tt. 


Equation  (35)  represents  the  instability  threshold  field  expected  for 
transverse  pumping  in  a  saturated  ferromagnetic  ellipsoid  with  a  linearly 
polarized  microwave  field  along  a  principal  axis  in  the  x-direction  and  the 
static  field  along  a  principal  axis  in  the  z -direction.  The  (u  +  Y)  and  (u  -  Y) 


15 


terms  in  the  denominator  of  Eq.  (34)  may  be  viewed  as  due  to  the  Larmor 
and  anti-Larmor  components  of  the  linearly  polarized  field  respectively. 
Previously,  the  usual  procedure  for  obtaining  hcrjt  for  a  linearly  polarized 
transverse  pump  field  has  been  to  apply  the  Suhl  theory  to  the  Larmor 
rotating  component  of  the  field  and  neglect  the  anti-Larmor  component. 

Such  a  procedure  would  result  in  Eq.  (34)  with  the  second  term  in  the 
curly  brackets  omitted.  The  present  result  shows,  however,  that  this 
anti-Larmor  term  has  a  noticeable  effect  on  the  threshold  expression. 
First,  the  anti-Larmor  term  decreases  h  ^  from  the  Suhl  prediction. 
Second,  when  this  term  is  included  the  static  field  dependence  of  the  denom¬ 
inator  cannot  be  separated  from  the  0,  dependence  so  that  for  k  *  0, 
is  a  function  of  static  field.  In  the  Suhl  theory,  is  independent  of 

static  field  for  k  *  0  and  only  a  function  of  u>Aj  .  Finally,  h  is  not 
independent  of  the  ^  spin-wave  angle.  In  order  to  illustrate  these  effects, 
specific  results  for  spherical  samples  are  shown  in  Fig.  1.  Normalized 
curves  of  are  shown  as  a  function  of  applied  field  Hq  normalized 

to  (u/y)  with  um/u  as  a  parameter.  The  solid  curves  correspond  to  hcrit 
when  the  anti-Larmor  contribution  is  included  and  the  dotted  curves  were 
calculated  neglecting  this  contribution.  The  two  effects  pointed  out  above 
are  apparent  from  these  curves.  1)  For  a  linear  pump  field,  hcrit  is 
reduced  somewhat  from  the  value  predicted  by  neglecting  the  anti-Larmor 
component  of  the  drive  and  applying  the  Suhl  theory,  particularly  at  low 
field.  2)  For  static  field  less  than  that  for  the  minimum  threshold, 
hcrit/AHk  is  not  a  linear  function  of  applied  field  Hq  as  expected  from  the 
Suhl  theory. 


In  order  to  demonstrate  the  consequences  of  the  relaxation  of  the 
rotational  symmetry  requirement  m  the  theory,  curves  of  h^^/AH^  as  a 
function  of  Hq/(u/y)  for  perpendicular  pumping  with  linearly  polarized 
microwave  drive  were  also  calculated  for  samples  with  planar  geometry 
magnetized  in  their  plane.  These  results  are  shown  in  Fig.  2.  Two  cases 
were  considered;  1)  hQ  parallel  to  the  plane  of  the  slab;  2)^  perpendicular 
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to  the  plane  of  the  slab.  In-plane  demagnetizing  factors  of  zero  and  an 
out-of-plane  demagnetizing  factor  of  unity  were  assumed.  For  h^ 
perpendicular  to  the  plane  of  the  slab,  the  threshold  amplitude  is  signif¬ 
icantly  higher  than  for  h  in  the  plane.  This  result  can  be  explained  by 
considering  the  effect  of  geometry  on  the  uniform  precession  amplitude 
aQ,  which  is  the  important  parameter  in  determining  the  threshold  for 
perpendicular  pumping.  For  a  parallel  b^,  Y  is  given  by  Y  *  yH  +  am. 

For  a  normal  h  ,  Y  s  yH  is  satisfied.  From  Eq.  (32a),  then,  it  can  be 
seen  that  is  larger  for  parallel  to  the  plane  of  the  slab  than  for Jj0 
normal  to  the  slab.  Consequently,  a  normal  is  less  efficient  in  pumping 
the  uniform  precession  than  an  in-plane  h Q  and  the  hcr for  a  microwave 
field  perpendicular  to  the  slab  is  higher.  Such  differences  are  not  acces¬ 
sible  when  rotational  symmetry  is  assumed  because  all  transverse  direc¬ 
tions  are  then  equivalent. 

Case  (2).  The  instability  threshold  for  perpendicular  pumping  with 
a  circularly  polarized  microwave  field  rotating  in  an  ant.i-Larmor  sense 
can  be  obtained  by  omitting  the  (<j  +  Y)  term  in  Eq.  (33)  and  multiplying 
the  expression  by  two.  The  factor  of  two  compensates  for  the  factor  of 
one-half  between  the  amplitude  of  a  linearly  polarized  field  and  the  ampli¬ 
tude  of  its  Larmor  or  anti-Larmor  rotating  components.  For  comparison, 
the  instability  threshold  for  Larmor  pumping  may  be  evaluated  by  setting 
the  (u  -  Y)  term  in  Eq.  (33)  to  zero,  as  in  the  Suhl  theory,  and  including 
the  factor  of  two.  Theoretical  curves  for  both  cases  are  shown  in  Fig.  3. 
The  upper  set  of  curves  correspond  to  anti-Larmor  pumping  and  the  lower 
set  of  curves  to  Larmor  pumping.  The  curves  were  calculated  for  spherical 
samples.  The  results  in  Fig.  3  illustrate  two  important  points  concerning 
anti-Larmor  pumping.  1)  The  threshold  is  significantly  higher  than  for 
Larmor  pumping.  (Note  that  the  h  axis  is  logarithmic.  )  2'  The 
threshold  does  not,  in  general,  exhibit  the  usual  minimum  which  is  charac¬ 
teristic  of  Larmor  pumping,  although  the  anti-Larmor  curve  does  exhibit. 

a  small,  local  minimum  for  large  u  /w. 

a  m 
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Case  (3).  Finally,  consider  oblique  pumping  with  a  linearly 
polarized  microwave  field  making  an  angle  *  with  the  direction  of  the 
static  external  field  and  polarized  in  the  x-z  plane.  This  field  corres¬ 
ponds  to  a  a  0,  6  a  0,  a  »  sin  *,  and  a  s  cos  *.  Also  assume,  for 
y  x  x  z 

simplicity,  that  the  static  field  is  along  an  axis  of  rotational  symmetry 
so  that  X  and  Y  are  both  equal  to  the  natural  uniform  precession 
frequency  uQ.  With  these  assumptions, 


_  ,  sin  * 

u  -  w 
o 

and 

q  . 

uQ  +  w 


(37a) 


(37b) 


For  this  situation,  W  in  Eq.  (28)  simplifies  to 


W  =  sin  0^  cos 


sin  *  [f(0,  )  +u/2]  -i<j>,  sin  *  [f(0.  )  -u/ZJ  i<f>. 


pin*[f(flk)- 

1 


e  K  + 


u  +  w 
o 


} 


+  sI”2#kaz 


(38) 


As  pointed  out  previously;  f($k)  “  w/Z  is  negative.  For  first-order  processes, 

moreover,  the  sample  is  biased  below  ferromagnetic  resonance  so  that  u 

±i$k  ° 

is  less  than  u.  The  coefficients  of  the  e  terms  in  the  curly  brackets  of 
Eq.  (38)  are  both  negative.  Phases  are  matched  for  all  terms  only  for 
#k  3  it  so  that  |wj  is  maximized  and  hc  minimized  with  respect  to  4>k  at 
4>k^m*n^  =  w.  This  result  is  strictly  valid  only  if  the  damping  of  the  uniform 
precession  is  neglected.  If  uQ  is  assigned  an  imaginary  part  tj  q  corres¬ 
ponding  to  the  uniform  precession  relaxation  rate  as  described  in  the  last 
section,  the  phases  of  the  terms  in  the  curly  brackets  of  Eq.  (38)  are  modi¬ 
fied.  The  change  is  small  for  small  damping. 

If  the  uniform  precession  damping  is  neglected,  is  exactly  ir 

and  |w|  is  given  by 
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\rnek.  ^min)  )| 


sin  9, 


(cos  0.  sin  *  ?  ^ 

— ^ - 2 -  1 2«o  +  u  1  +  COs  *  sin  djc  f  • 

u  -  U).  J 


(39) 


The  threshold  hcrit  is  then  obtained  by  minimizing  |w|  *  with  respect  to  9^ 
as  discussed  above.  In  Fig.  4,  normalized  curves  of  as  a 

function  of  H^fu/y)  are  shown  with  *  as  a  parameter.  Spherical  samples 
and  «m/u  *  0*  54  corresponding  to  yttrium  iron  garnet  at  9.  2  GHz  were 
assumed.  The  curve  for  #  *  0°  corresponds  to  parallel  pumping  and  the 
♦  a  90°  curve  corresponds  to  perpendicular  pumping.  The  curves  illustrate 
the  dependence  of  ilcrit/AHj{  on  pumping  angle.  Actual  h  ^  curves  would  be 
slightly  modified  because  of  the  k  dependence  of  AH.  which  influences 


(min) 

'k 

angles  near  *=  90°.  At  high  field,  h  ^AH^  decreases  with  *  and  is 
most  sensitive  to  angle  near  *  a  0°.  The  threshold  field  minimum  moves 
to  higher  static  field  with  increasing  *. 


0,'*“*“'.  At  low  field,  h  .^/AH^  increases  with  *  and  is  most  sensitive  to 

l  .J 

crifr 


SUMMARY  AND  CONCLUSION 


Spin-wave  instability  theory  has  been  extended  to  include  pumping 
with  an  arbitrary  microwave  field  configuration  and  ellipsoidal  samples 
magnetized  along  a  principal  axis  which  is  not  necessarily  an  axis  of  rota¬ 
tional  symmetry.  The  calculations  have  been  carried  out  explicitly  for  the 
first-order  instability  threshold.  Several  special  cases  have  been  considered 
to  demonstrate  the  versatility  of  the  theory  in  analyzing  situations  previously 
inaccessible  on  the  basis  of  the  original  theoretical  treatment.  Sample 
geometry  and  the  pump  configuration  have  a  strong  influence  on  the  instability 
threshold  amplitude. 
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FIGURE  CAPTIONS 


Fig.  1  First-order  spin-wave  instability  threshold  field  hcrit 

normalized  to  the  spin-wave  linewidth  AH^  as  a  function  of 

the  static  external  field  Hq  normalized  to  w/y,  the  angular 

excitation  frequency  divided  by  the  gyromagnetic  ratio. 

Curves  were  obtained  for  spherical  samples  and  a  linearly 

polarized  microwave  excitation,  with  wm/u  (um  3  y4trM  is 

the  saturation  induction)  as  a  parameter.  The  dotted  curves 

neglect  the  contribution  of  the  anti-Larmor  component  of  the 

microwave  field  to  h  ... 

cnt 

Fig.  2  Normalized  first-order  spin-wave  instability  threshold  field 

h  as  a  function  of  the  normalized  static  external  field  H  . 
crit  o 

Curves  are  shown  for  thin  disk  shaped  samples  with  an  in¬ 
plane  demagnetizing  factor  of  zer  and  linearly  polarized  micro- 
wave  excitation  parallel  to  (hQ!!)  or  perpendicular  to  (hQl)  the 
slab  plane.  The  slab  was  assumed  to  be  magnetized  to  satura¬ 
tion  in  its  plane  and  perpendicular  to  the  microwave  field  direc¬ 
tion.  Sets  of  curves  are  shown  for  w  /u  *  l/2  and  u  /u  =  1. 

m  m 

Fig.  3  Normalized  first-order  spin-wave  instability  threshold  field 

h _ ..  as  a  function  of  the  normalized  static  external  field  H  . 

cnt  o 

Curves  were  calculated  for  spherical  samples  and  a  circularly 
polarized  microwave  excitation  perpendicular  to  the  static  field 
direction,  with  um/w  as  a  parameter.  Results  are  shown  for  a 
Larmor  excitation,  rotating  in  the  proper  sense  to  excite  ferro¬ 
magnetic  resonance,  and  an  anti-Larmor  excitation,  rotating  in 
the  opposite  sense. 


Fig.  4  Normalized  first -order  spin-wave  instability  threshold  field 

h  ..  as  a  function  of  the  normalized  static  external  field  h  . 
cnt  o 

Curves  were  calculated  for  spherical  samples  and  a  linearly 
polarized  microwave  field  at  ingles  of  0°,  30°,  60°,  and  90° 
to  the  static  field  direction.  The  values  indicated  correspond 
to  wm/w  =  0.  54  which  is  appropriate  for  yttrium  iron  garnet  at 
9.  2  GHz. 
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The  threshold  microwave  field  amplitude  hcrjt  required  for  the 
unstable  growth  of  spin  waves  has  been  investigated  for  oblique  pumping 
with  a  linearly  polarized  microwave  field  applied  at  an  arbitrary  angle 
with  respect  to  the  static  magnetic  field.  The  experiment  was  performed 
on  YIG  (1  percent  Dy)  spheres  at  9.  2  GHz  using  a  rectangular  TEJ02 
cavity  and  a  pulsed  magnetron  source.  Data  were  obtained  at  static 
external  fields  of  400,  600,  1100,  and  1500  Oe.  In  general,  hcrit  increases 
with  the  pump  angle.  Theoretical  expressions  for  the  oblique  pumping 
hcrit  have  been  derived  for  ferromagnetic  insulators  of  ellipsoidal  shape 
and  with  axial  symmetry  about  the  direction  of  the  static  magnetic  field 
sufficient  to  saturate  the  sample.  If  the  difference  between  the  angular 
pump  frequency  divided  by  the  gyromagnetic  ratio  and  the  static  internal 
field  is  much  greater  than  the  ferromagnetic  resonance  linewidth  at  the 
pump  frequency,  the  threshold  field  is  determined  by  maximizing  l/h  = 
cos  tj</h||  +  sin  <J»/h^  with  respect  to  the  angle  between  the  internal  field  and 
the  spin  wave  propagation  direction.  Here,  hn  and  h^  are  the  parallel 
and  perpendicular  pump  threshold  fields  of  SchlSmann  and  Suhl.  For 
external  fields  sufficient  to  magnetize  the  spherical  samples  to  saturation, 
the  theoretical  predictions  are  in  good  agreement  with  the  experimental 
data.  For  smaller  fields,  the  agreement  is  poor.  However,  fair  agree¬ 
ment  is  obtained  by  assuming  a  transverse  demagnetizing  factor  Nt  greater 
than  the  value  of  l/3  which  is  appropriate  for  saturated  spherical  samples. 
The  value  of  for  the  best  theoretical  fit  depends  on  the  way  in  which 
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the  sample  is  demagnetized.  This  result  indicates  that  for  samples  not 
magnetized  to  saturation,  the  magnetic  structure  consists  of  elongated 
domains  and  that  the  oblique  pumping  data  can  be  explained  from  the 
theory  for  saturated  systems  assuming  demagnetizing  factors  appropri¬ 
ate  to  the  domain  structure. 
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INTRODUCTION 


The  parallel  and  perpendicular  pumping  of  spin  waves  in  insulating 
ferromagnets  has  received  considerable  attention  in  recent  years.  Suhl* 
has  calculated  the  threshold  microwave  field  amplitude  required  for  the 
unstable  growth  of  spin  waves  with  the  microwave  field  perpendicular  to 
the  static  magnetic  field  for  saturated  samples  of  ellipsoidal  shape  with 
a  rotational  axis  of  symmetry  parallel  to  the  static  field.  SchlSmann  et  al.  ^ 
have  considered  the  situation  for  the  microwave  field  parallel  to  the  static 
field.  There  has  been  little  work  concerning  the  threshold  fields  for  the 
unstable  growth  of  spin  waves  under  so-called  oblique  pumping  with  a 
linearly  polarized  microwave  field  at  an  angle  Osi f>s  90°  with  the  static 

field  direction.  Stern^  has  discussed  oblique  pumping  as  related  to  remanent 

4  - . — 

ferrite  device  operation.  Kitaev  et  al.  have  considered  the  theoretical 

problem  of  pumping  with  an  arbitrary  pump  field,  although  specific  results 
were  shown  only  for  an  infinite  medium.  An  experimental  and  theoretical 
investigation  of  oblique  pumping  is  the  subject  of  the  present  report.  The 
experimental  results  for  saturated  spherical  samples  are  in  good  agree¬ 
ment  with  the  theory.  For  applied  static  fields  insufficient  to  saturate  the 
samples,  fair  agreement  can  be  obtained  by  assuming  transverse  demag¬ 
netizing  factors  which  correspond  more  closely  to  the  domain  shape  than 
the  sample  shape. 
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EXPERIMENT  AND  RESULTS 


The  spin  wave  instability  threshold  for  oblique  pumping  was 

determined  from  microwave  susceptibility  data  obtained  using  cavity 

perturbation  techniques  and  a  rectangular  TEj0£  reflection  cavity  at 

9.  2  GHz.^’  ^  Determinations  of  h  ..  were  made  as  a  function  of 

o  o  crit 

pumping  angle  ip  for  0  £  <p  £  90  for  a  one  percent  dysprosium  doped 

polycrystalline  yttrium  iron  garnet  sphere  with  a  linewidth  of  50  Oe  at 

9.2  GHz  and  a  saturation  magnetization  of  1750G.  This  particular 

sample  was  chosen  for  two  reasons:  (1)  The  rare  earth  contribution 

to  the  spin-wave  linewidth  should  dominate  over  other  contributions, 

thereby  facilitating  the  theoretical  interpretation  of  the  data,  and  (2) 

the  spin-wav;  linewidth  is  still  small  enough  to  allow  measurement  of 

the  threshold  field  h  .4  with  the  maximum  attainable  microwave  field 

crit 

amplitude  of  80  Oe. 


In  Figs.  1  to  4,  the  threshold  field  hcrjt  is  shown  as  a  function  of 
pumping  angle  ip  for  external  static  fields  of  400,  600,  1100,  and  1500 
Oe  respectively.  At  all  but  the  highest  field,  hcrit  increases  with  pump¬ 
ing  angle.  These  data  are  in  qualitative  agreement  with  previous  meas- 
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urements  by  Green  at  al.  which  indicated  that  the  parallel  pump  threshold 
corresponding  to  ip  =  0  was  significantly  lower  than  that  for  ip  =  tt/2  and 
perpendicular  pumping.  The  present  data  on  threshold  field  as  a  function 
of  pump  angle  provide  several  new  conclusions  concerning  spin-wave 
instabilities  in  ferromagnetic  insulators.  First,  these  data  can  be  inter¬ 
preted  in  terms  of  the  previous  calculations  for  parallel  and  perpendicular 
pumping,  if  the  theories  are  extended  to  include  oblique  pumping.  Second, 
the  data  for  samples  in  a  nonsaturated  state  are  also  in  lualitative  agree¬ 
ment  with  the  theory  developed  for  saturated  systems  if  transverse  demag¬ 
netizing  factors  appropriate  to  the  domain  shape  and  not  the  sample  shape 
are  used.  In  the  next  section,  the  theory  of  oblique  pumping  will  be  pre¬ 
sented  and  in  the  following  section,  the  above  data  will  be  interpreted  in 
terms  of  the  theory. 
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THEORY 


The  normal  magnetic  modes  for  ferromagnetic  samples  of  ellip- 

8 

soidal  shape,  have  been  calculated  by  Holstein  and  Primakoff,  who 
obtain  a  dispersion  relation  of  the  form: 

wk  *  Y  fCH.  +  Dk2)  (H.  +  Dk2  +  4trMa  sin2  0R)  ]  *  (1) 

where  is  the  angular  frequency  of  a  spin-wave  with  wave  vector  of  mag¬ 
nitude  k  and  at  an  angle  0.  with  the 

k  7 

direction  of  the  internal  field  H.,  y  gyromagnetic  ratio  (1.76  X  10 

-1  -1  1  -9 

Oe  -sec  ),  D  is  a  phenomenological  exchange  parameter  (*  5  X  10 

2 

Oe-cm  for  YIG),  and  4irM  is  the  saturation  induction  (  a  1750  G).  In 
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Fig,  5,  is  shown  as  a  function  of  k  for  0^  =  0  and  ir/2,  as  calculated 
for  YIG  from  Eq.  (1)  with  an  internal  field  of  180  Oe.  The  region  bounded 
by  these  two  curves  is  called  the  spin-wave  manifold. 

If  a  linearly  polarized  microwave  field  helut  of  frequency  u  is 
applied  to  the  ferromagnetic  sample,  perpendicular  to  the  static  internal 
field  (perpendicular  pumping),  Zeeman,  dipolar,  and  exchange  interactions 
give  rise  to  a  coupling  between  the  uniform  precession  spin  wave  at  u 
driven  by  the  microwave  field  and  other  spin-wave  states.  For  sufficiently 
intense  microwave  drive,  this  coupling  leads  to  the  unstable  growth  of 
certain  spin-wave  amplitudes.  If  the  microwave  field  is  applied  parallel 
to  the  static  internal  field  (parallel  pumping),  the  uniform  precession  ampli¬ 
tude  is  zero.  For  parallel  pumping,  however,  spin-wave  states  are  coupled 
directly  to  the  microwave  drive  and  spin-wave  instabilities  may  also  occur 
for  sufficiently  intense  microwave  field  amplitudes.  For  the  special  case 
of  saturated  ferromagnetic  insulators,  in  the  shape  of  ellipsoids  of  revolu¬ 
tion  about  the  internal  field  direction,  Suhl*  has  calculated  the  threshold 

field  h  .  required  for  the  unstable  growth  of  the  spin -wave  amplitude,  for 

-*•  2 
state  (u^,  k,  0^)  under  perpendicular  pumping  and  Schl5mannetal.  have  calculated 

the  corresponding  threshold  h()  for  parallel  pumping.  The  calculations 

indicate  that  the  threshold  is  particularly  low  if  =  u/2  and  the  pumping 

of  unstable  spin  waves  at  u/Z  has  been  termed  a  first  order  process  (strong 

coupling).  For  a  first  order  process,  the  results  are 


5 


(2a) 


where  <*>m  is  given  by  y4irMg,  uQ  is  the  uniform  precession  resonance 

frequency  u  =  vH.  +  N.u  for  a  given  internal  field  H.  and  transverse 
•'o  itm  1 

demagnetizing  factor  Nt>  AH^  is  the  spin -wave  linewidth,  and  AH  is  the 
uniform  precession  linewidth  at  frequency  u.  The  particular  magnon 
which  has  the  lowest  h^  for  perpendicular  pumping  or  h jj  for  parallel 
pumping  will  be  the  one  which  goes  unstable  first  and  will  control  the 
experimentally  determined  threshold  field  h  For  perpendicular 
pumping  the  threshold  field  corresponds  to  spin  waves  with  0^  »  w/4 
and  for  parallel  pumping  0^  =  ir/2  spin  waves  determine  the  threshold. 
These  two  states  are  indicated  as  black  dots  in  Fig.  5  for  an  excitation 
frequency  of  9.  2  GHz.  In  the  above  and  in  the  discussion  to  follow,  it 
is  assumed  that  u  is  sufficiently  high  or  H.  is  sufficiently  low  so  that 
u/2  is  above  the  0  =  rr/2  branch  of  the  manifold  at  k  =  0.  The  theory  is 
strictly  valid  only  for  k  >  0. 


In  an  oblique  pumping  experiment  both  a  parallel  microwave  field 
component  helu*  cos  tj;  and  a  perpendicular  microwave  field  component 
helu*  sin  i|»  are  present  and  the  possibility  exists  that  some  magnon  state 
with  tt/4  <  0^  <  tr/2  determines  the  threshold  field.  Such  a  state  is  repre¬ 
sented  by  the  open  circle  in  Fig.  5.  The  oblique  pumping  threshold  field 
can  be  calculated  in  a  straightforward  manner  by  considering  the  hamil- 
tonian  for  the  magnetic  system  in  the  presence  of  an  oblique  pumping  field, 
calculating  the  spin-wave  amplitude  equation  of  motion,  and  proceeding 
according  to  the  technique  developed  by  Suhl  and  SchlSmann.  For  the 
present  purposes,  however,  it  is  more  illustrative  of  the  physical  pro¬ 
cesses  involved  to  calculate  the  threshold  from  a  circuit  analogy,  as 
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reported  previously  by  one  of  the  authors.  The  spin-wave  theory  is 
outlined  in  the  appendix.  Assuming  linear  coupling  between  the  pumping 
field  and  the  spin  system,  as  illustrated  by  the  equivalent  circuit  in 
Fig.  6,  and  taking  the  phases 
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of  the  perpendicular  and  parallel  coupling  coefficients,  and  C((,  to 
be  equal,  the  effective  spin-wave  drive  U  may  be  defined  as 

U  ■  Cjj  h  cos  h  sin  i|»  (3) 

where  Cj|  and  depend  on  the  particular  magnon  state  (w/2,  k,  0^) 
being  pumped.  In  Eq.  (3),  h  cos  t|>  and  h  sin  tj;  are  the  parallel  and  per¬ 
pendicular  components  of  the  microwave  field  amplitude,  respectively. 

The  instability  for  a  particular  spin  wave  occurs  when  U  exceeds  an 
effective  drive  threshold  T  which  can  be  determined  by  considering 
the  known  theoretical  threshold  for  that  spin  wave  from  Eq.  (2)  for 
i|»  *  0°  and  if*  «  90°.  A  parallel  pumped  instability  occurs  for  T  *  C|jh|j 
and  a  perpendicular  pumped  instability  for  T  a  C^h^.  With  a  T/h() 
and  a  T/h^,  Eq.  (3)  becomes 

1  _  T  /  cos  4»  L  sin  v 
E  *  0  (“E^  +  ~Ep') 

Instability  occurs  when  U  a  T  is  satisfied  so  that  the  threshold  field  ampli¬ 
tude  hoblique  required  to  drive  state  (u/2,  k.  0^)  unstable  is  given  by 

1  _  cosi|i  .  sin  4  — .  ,,, 

c -  =  +  ~L  T-  •  (4) 

"oblique  nll  °T 

In  the  appendix  Eq.  (4)  is  obtained  from  spin-wave  theory.  It  is  shown 
that  the  equal  phase  assumption  can  be  satisfied  by  correctly  choosing 
the  angle  <6^  between  the  component  of  £  in  the  plane  normal  to  Hq  and 
the  transverse  component  of  the  microwave  field. 


The  oblique  pumping  threshold  field  hcr^t  is  obtained  by  maximizing 
the  right-hand  side  of  Eq.  (4),  thereby  minimizing  ^ob^qUe’  If  the  spin- 
wave  linewidth  AH^  is  independent  of  k  and  0^,  hcr^t  as  a  function  of 
at  constant  external  field  can  be  calculated  from  Eqs.  (2)  and  (4).  If  AH^ 
depends  on  k  and  0^,  the  minimization  procedure  necessary  to  determine 
the  spin-wave  state  (u/2,  k,  0^)  which  gives  the  lowest  threshold  hcrjt  is 
more  complicated.  An  approximate  procedure  was  used  to  take  the  k  and 
0 ^  dependence  of  into  account. 
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It  was  assumed  that  the  k  and  0 k  dependences  of  AHk  could 
be  expressed  in  a  form 


AHk  «  fj(k)  •  f2(0k)  (5) 

with 

f,(k)  ■  A+Bk  (6) 

and 

f2(k)  a  1  +  C{tr/2  -  6k)  +  D(rr/2  -  9^2  .  (7) 


Both  parallel  and  perpendicular  data  indicate  that  AHk  is  approximately 
an  increasing  linear  function  of  k.  The  A  and  B  parameters  were 


obtained  Jfrom  AHk  vs  k  determinations  on  the  basis  hcrjt  data  for 
parallel  pumping  on  the  Dy  doped  material.  These  results  are  shown 


in  Fig.  7.  The  A  parameter  was  deter  uined  to  be  15.  0  and  B  is 
-5 

0.  5  X  10  .  The  parameters  in  (0 were  determined  by  fitting  the 

theoretical  curve  of  hQ^^.C}ue  vs  pump  angle,  calculated  from  Eqs.  (2), 
(4),  and  (5)  -  (7),  to  the  experimental  hcr^t  data  for  a  static  external 
field  of  1100  Oe.  Using  this  procedure,  C  was  determined  to  be  0.  15 


and  D  was  -0.  3.  For  a  heavy  rare  earth  doping,  AHk  is  expected  to 
be  approximately  constant  and  not  exhibit  strong  k  or  0k  dependences. 
This  expectation  is  satisfied  by  the  above  result.  The  theoretical  results 


would  probably  not  be  severely  modified  if  AHk  was  assumed  to  be  constant. 
Consequently,  the  specific  form  for  the  weak  k  and  0k  dependences  should 
not  be  critical. 
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DISCUSSION 


First  consider  the  data  for  external  fields  sufficient  to  magnetize 
the  test  sample  nearly  to  saturation,  as  shown  in  Fig.  3  for  an  external 
field  of  1100  Oe  and  in  Fig.  4  for  a  field  of  1500  Oe.  As  discussed  in 
the  last  section,  a  qualitative  best  fit  between  the  theory  and  data  at 
1100  Oe  was  used  to  determine  the  parameters  A  and  B  in  the  0^  depen¬ 
dence  of  AH^.  As  evident  from  Fig.  3,  this  two  parameter  fit  gives 
good  agreement  with  the  data.  The  theoretical  curve  was  obtained  using 
a  transverse  demagnetizing  factor  of  1/3,  appropriate  for  saturated 
spherical  samples.  If  oblique  pumping  simply  consisted  of  simultaneous 
but  independent  parallel  and  perpendicular  pumping,  one  would  expect  the 
observed  threshold  to  correspond  to  either  h||/cos  ^  or  h^/sin  i|»  whichever 
was  lower.  The  independent  pumping  model,  therefore,  predicts  an 
increase  in  h  ^  with  ij»  to  about  60°,  as  indicated  by  the  h((/cos  ip  curve, 
and  a  subsequent  decrease  at  higher  angles,  as  indicated  by  the  h^/sin  if* 
curve.  The  agreement  between  the  data  in  Fig.  3  and  the  independent 
pumping  model  is  rather  poor.  In  Fig.  4,  the  calculated  hcrit  vs  «|>  curve 
for  a  static  field  of  1500  Oe  is  shown  for  comparison  with  the  data.  The 
theoretical  curve  was  obtained  using  the  AH^.  0 ^  and  k  dependence 
determined  from  the  fit  at  1100  Oe.  There  are  no  adjustable  parameters 
available  to  adjust  the  fit  between  the  data  and  theory  for  Fig.  4.  The 
agreement  is  quite  good.  The  close  agreement  between  the  theoretical 
curves  and  the  data  in  Figs.  3  and  4  indicate  that  for  saturated  ferro¬ 
magnetic  insulators  in  the  shape  of  ellipsoids  of  revolution  about  an  axis 
parallel  to  the  static  field,  the  first  order  instability  threshold  for  oblique 
pumping  can  be  obtained  from  the  extension  of  previous  theories  of  parallel 
and  perpendicular  pumping. 

For  external  fields  less  than  4ffMs/3,  or  about  600  Oe  for  YIG, 
spherical  samples  are  not  magnetized  nearly  to  saturation,  and  the 
samples  consist  of  discrete  domains.  In  the  presence  of  domain  struc¬ 
ture,  the  theory  is  not  strictly  applicable.  Nevertheless,  it  is  informative 
to  compare  the  data  for  external  fields  below  4TrMg/3  with  the  theoretical 
predictions.  Oblique  pumping  h  data  as  a  function  of  are  shown  in 
Figs.  1  and  2  for 
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external  fields  of  400  and  600  Oe,  respectively.  First  consider  the  data 
at  400  Oe.  The  upper  theoretical  curve  was  calculated  assuming  a 
transverse  demagnetizing  factor  of  l/3,  appropriate  for  saturated  spheres. 
The  agreement  between  the  data  and  this  theoretical  curve  is  extremely 
poor.  Better  agreement  can  be  obtained,  however,  if  the  transverse 
demagnetizing  factor  is  assumed  to  be  l/2,  as  shown  by  the  lower  solid 
curve.  With  this  assumption  the  shape  of  the  theoretical  curve  closely 
matches  the  data  although  the  values  are  somewhat  lower.  Some  varia¬ 
tion  in  the  experimental  data  has  been  noted  which  depends  on  the  way  in 
which  the  sample  was  demagnetized.  In  general,  a  good  fit  to  any  particu¬ 
lar  set  of  data  could  be  obtained  with  l/3  <  <  l/2.  Stern  has  previously 

reported  data  on  partially  magnetized  toroids  which  agree  with  theory  for 
*  l/2.  Agreement  for  1/3  <  Nt  <  l/2  indicates  that  the  portions  of  the 
material  which  go  unstable  first  correspond  to  elongated  domains.  Such 
domains  would  have  Nt»  l/2.  Further  support  for  the  proposal  that  the 
domain  structure  is  important  in  determining  the  instability  threshold  is 
provided  by  the  data  at  600  Oe,  shown  in  Fig.  2.  At  600  Oe,  the  internal 
field  is  approximately  zero  for'YICT and  the  material  is  nearly  saturated. 
The  theoretical  curve  shown  in  Fig.  2  which  provided  the  best  fit  to  the 
data  was  obtained  for  Nt  *  0.  38,  slightly  greater  than  the  value  of  l/3 
expected  for  saturated  spheres  due  to  the  residual  domain  structure  at 
600  Oe.  It  should  be  emphasized  that  the  same  AH^  dependence  which  was 
used  to  obtain  the  theoretical  curves  for  external  fields  greater  than 
4irMg/3  was  used  to  calculate  the  curves  shown  in  Figs.  2  and  1.  The  only 
adjustable  parameter  was  Nt.  Demagnetization  effects  due  to  domain 

structure  have  been  used  to  explain  large  low  field  loss  in  ferromagnetic 
a 

insulators.  The  present  data  and  theoretical  agreement  indicate  that 
domain  structure  also  influences  the  instability  threshold.  Demagnetizing 
factors  appropriate  to  the  domain  structure  and  not  the  sample  shape  should 
be  used  in  theoretical  estimates  of  hcrit  for  partially  magnetized  samples. 
Any  quantitative  comparison  for  systems  below  saturation,  however,  will 
require  extension  of  the  theory  to  include  samples  which  are  not  ellipsoids 
of  revolution  and  the  effect  of  the  antilarmor  component  of  the  linear  drive, 
as  well  as  more  precise  information  concerning  domain  structure. 
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SUMMARY 


This  investigation  has  provided  new  information  and  new  under¬ 
standing  concerning  the  microwave  field  thresholds  for  the  excitation 
of  spin-wave  instabilities  in  ferromagnetic  insulators.  (1)  Extensive 
data  on  the  oblique  pumping  threshold  as  a  function  of  pumping  angle 
have  been  obtained  for  external  fields  both  above  and  below  that  re¬ 
quired  to  magnetize  the  sample  to  saturation.  (2)  The  theory  of 
parallel  and  perpendicular  pumping  in  saturated  ferromagnetic  insu¬ 
lators  has  been  extended  to  include  oblique  pumping.  (3)  The  data  for 
saturated  samples  can  be  explained  by  the  theory.  (4)  The  data  for 
unsaturated  samples  can  be  qualitatively  explained  by  the  theory  if 
transverse  demagnetizing  factors  appropriate  to  the  domain  structure 
and  not  the  sample  geometry  are  assumed. 


V 


APPENDIX 


Following  the  treatment  by  SchlSmann,  ^  the  spin-wave  instability 
threshold  field  amplitude  for  oblique  pumping  can  be  obtained  by  con¬ 
sidering  a  hamiltonian  for  the  magnetic  system  of  the  form 

K  *  I"kbkbk 

k 

+  r  1  <fkbobkb-k  +  «■'•> 

k*0 

-  ^  COS  wt  ^  (°kbJ  b!k  +  C.  C. )  (Al) 

k 


where  is  the  spin-wave  frequency,  the  b^  are  spin-wave  amplitudes, 
f^  is  a  coupling  coefficient  for  perpendicular  pumping  (ij<  a  90°),  and  c^  is 
a  coupling  coefficient  for  parallel  pumping  (ip  =  0°)  assuming  a  cos  ut  time 
dependence  for  the  component  of  the  microwave  field  along  the  static  field. 
With  parallel  and  perpendicular  microwave  field  components  h  cos  ut  cos 
and  h  cos  ut  sin  i|j,  respectively,  the  coupling  coefficients  of  Eq.  (Al)  are1 1 


fk*-^  fYH.-lvDk^Je"1** 


sin0k  cos  0^ 


(A2a) 


2*$k  2 

yh  cos  e  sin  0, 


(A  2b) 


The  perpendicular  component  of  the  microwave  field  is  along  the  x-axis 
and  is  measured  with  respect  to  this  axis.  The  bQ  in  Eq.  (Al)  is  related 
to  the  uniform  precession  amplitude  and  is  given  by 


b 


o 


yh_sin  Jj 

2~(w0-(JT 


iut 

t; 


(A3) 
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A 


if  the  antilarmor  term  in  the  transverse  component  of  the  linearly 
polarized  microwave  field  is  neglected.  The  f^  also  neglects  the 
contribution  of  this  term  to  the  coupling.  For  sufficiently  large 
internal  fields,  the  correction  due  to  these  antilarmor  terms  is 
small  and  may  be  neglected. 

The  first-order  instability  threshold  is  obtained  from  the 
equation  of  motion  for  the  b^, 

dbk/dt  a  i  dX/d  b*  (A4) 

*  1  l°kbk  +  *fk  bo  '  ck  COS  ut*  b-k  J 


where  a  +  i  77^.  Damping  is  taken  into  account  by  assigning  an 
imaginary  part  77^  to  the  spin-wave  frequency  The  relaxation  rate 
77^  is  related  to  the  spin-wave  linewidth  by  77^  =  f  y^H^.  ^he  instaljility 
threshold  is  obtained  by  determining  the  microwave  field  amplitude  at 
which  the  equation  of  motion  (A4)  permit  a  stationary  solution 


»  b,  _  e 
ko 


iut/2 


in  the  presence  of  damping.  By  considering  only  terms  which  vary  as 
eiut/Z  db^/dt  and  e  2  in  db*k/dt,  two  homogeneous  linear  equations 
in  bko  and  b  *kQ  result.  The  microwave  field  amplitude  for  which  a  non¬ 
trivial  solution  to  bkQ  and  b_kQ  exists  corresponds  to  the  instability 
threshold.  From  Eqs.  (A 4)  and  (A5),  a  nontrivial  solution  for  bkQ  is 
obtained  only  for 


K  -«/z|  = !  4 


fk  vh  sin 
l(no  -  w) 


where  the  uniform  precession  frequency  has  been  replaced  by  in 

the  expression  for  b  .  The  relaxation  rate  77  is  related  to  the  uniform 

0  0 

precession  linewidth  All  by  77  =  ^  The  right-hand  side  of  Eq.  (.A 6)  is 

proportional  to  the  microwave  field  amplitude  h.  Instability  occurs  for 


the  smallest  value  of  h  if  the  spin-wave  frequency  uk  is  equal  to  one- 
half  the  pump  frequency  u.  If  u>k  *  u/2  is  satisfied,  the  critical  field 
is  determined  by 


\  »  ^  |cos  tj»  sin2  $k  e  *k 

(yH.  +  yDk2  +  u/2) 

+  sin  sin  $k  cos  0k  - ^  -  .  (A7) 

The  lowest  threshold  for  a  given  magnon  instability  is  realized  for  that 
value  of  <£k  for  which  the  phases  of  the  two  terms  inside  the  vertical  bars 
of  Eq.  (A7)  are  equal.  The  threshold  is  given  by 


um  sin2  9. 
m  k 


oblique 


cos 


u  sin  9,  cos  9,  (w/2  +  yDk  +  yH.) 

•  ^  iC  K  1  _•  • 

“AHk  l°o-“|  * 


which  is  equivalent  to  Eq.  (4). 


It  is  important  to  realize  that  ft  -u  is  negative  and  that  the  phase 

-1  °  1 
of  (ft  -u)  ,  which  is  specified  by  the  relative  magnitudes  of  n  =  •*•  yAH, 

o  o  c 

uQ,  and  u,  determines  the  <j>k  which  leads  to  the  lowest  threshold.  For 
narrow  linewidth  materials,  uQ  »  is  satisfied  for  sufficiently  large 
applied  fields  and  hQtjlic}ue  corresponds  to  4>  k  M  *.  In  general,  however, 
the  <£k  for  which  Eq.  (A8)  is  valid  is  given  by 

4»k  =  ir  +  tan  1  f  77 /(u  -  <jq)]  .  (A9) 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Experimental  and  theoretical  microwave  field  threshold 
amplitudes  for  a  YIG  (1  percent  Dy)  sphere  at  9.2  GHz  as 
a  function  of  the  pump  angle  between  the  static  magnetic 
field  direction  and  microwave  field  polarization  direction  for 
an  external  static  magnetic  field  of  400  Oe.  The  theoretical 
curves  (solid  lines)  were  obtained  assuming  transverse  de¬ 
magnetizing  factors  of  l/2  and  l/3,  as  explained  in  the  text. 

Experimental  and  theoretical  threshold  amplitudes  for  a 
YIG  (1  .percent  Dy)  sphere  at  9.  2  GHz  as  a  function  of  pump 
angle  for  an  external  static  field  of  600  Oe.  The  theoretical 
curve  was  obtained  with  a  transverse  demagnetizing  factor 
of  0.  38. 

Experimental  and  theoretical  threshold  amplitudes  for  a 
YIG  (1  percent  Dy)  sphere  at  9.  2  GHz  as  a  function  of  pump 
angle  for  an  external  static  field  of  1100  Oe.  The  lower 
theoretical  curve  was  obtained  by  adjusting  the  dependence 
of  the  spin-wave  damping  on  wave  vector  to  give  the  best  fit 
with  a  transverse  demagnetizing  factor  of  l/3.  The  two 
upper  curves  were  obtained  by  considering  parallel  and  per¬ 
pendicular  pumping  processes  which  operate  independently, 
as  explained  in  the  text. 

Experimental  and  theoretical  threshold  amplitudes  for  a 
YIG  (1  percent  Dy)  sphere  at  9.  2  GHz  as  a  function  of  pump 
angle  for  an  external  static  field  of  1500  Oe.  The  theo¬ 
retical  curve  was  obtained  with  no  adjustable  parameters, 
as  explained  in  the  text. 

Spin-wave  frequency  as  a  function  of  wave  number  mag¬ 
nitude  k  for  ^  parallel  to  the  internal  field  direction  (8^  =  0) 
and  perpendicular  to  the  internal  field  (0^  =  w/2).  The 
parallel  pump  instability  for  a  spin  wave  at  occurs  at 
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Figure  6 


Figure  7 


0^  ~  ant*  P^P^dicular  pump  instability  occurs  at 
9^  2  it/ 4  as  indicated  by  the  black  dots.  The  oblique  pump 
instability  corresponds  to  an  intermediate  state,  ir/4  <0^  < 
ir/2,  as  indicated  by  the  open  circle. 

Equivalent  circuit  for  the  oblique  pumping  of  a  spin-wave 
state  (u^,  k,  0^)  with  a  microwave  field  polarization  di¬ 
rection  at  an  angle  with  respect  to  the  static  internal 
field.  The  parallel  and  perpendicular  pumping  contribu- 
tions  are  characterized  by  sources  h  cos  4<  e  w  or  h  sin  i|j 
e  and  coupling  coefficients  or  C^,  respectively. 

The  experimentally  determined  k- dependence  of  the  spin- 
wave  line  width  for  parallel  pumping  in  a  1  percent  Dy 
doped  YIG  sphere  at  9.  2  GHz. 
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The  microwave  properties  of  partially  mangetized  materials  have  been 
investigated.  The  real  parts  («',  s',  y2 1 >  of  the  components  of  the 
permeability  tensor  have  been  measured  as  a  function  of  the  ratios  of 
saturation  magnetization  to  angular  frequency  'iu)  ana  average  mag¬ 

netization  (-LAM)  to  angular  frequency  and  location'  on  the  magnetic  hy¬ 
steresis  loop  on  a  yttrium  iron  garnet  (TIG)  and  magnesium  manganese 
ferrite  (MgMnP).,  The  off  dlagon  (^..component ,  <  '  ,  gives  good  agreement 
with  Rado's  prediction  that  k‘3H^ — .  The  diagonal  components  u',  y^'  ao 
not  fit  Rado's  prediction  and  the  cause  of  the  failure  of  Rado's  analysii 
for  y',  uz'  is  suggested.  The  dependence  of  y',  y  '  upon  a 

has  been  empirically  fitted  by  certain  mathematical  functions.  The 
imaginary  components  (y1',  c '  ’  ,  y-'1)  have  been  measured  for  .7c  Y. 

1  using  circularly  polarized  cavities.  For  T-*-wlls  «,d,  it  nas  bee'n  found 
[that  *’’<<  u  '  '  and  that  the  loss  of  the  partially  magnetized  state  ca 
I  be  characterized  by  a  single  value  of  u'',  that  for  the  completely  de¬ 
magnetized  state.  For  values  of  —.7  large  rectangular  samples 

and  exact  the  solution  to  the  waveguide  cavity  mode  have  been  utilized 
to  obtain  W ' '  for  the  demagnetized  state. 

The  high  power  threshold  has  been  measured  on  YIG  as  a  function  of 
sample  shape,  the  polarization  of  the  rf  magnetic  field  and  the  orienta¬ 
tion  of  the  rf  field  to  the  static  magnetization.  The  Suhl  and  Scnloe- 
mann  theories  for  hivJi  newer  non-linearity  lave,  been  .eitemlea  to  include 
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arbitrary  ellipsoids,  linearly  polarized 
rf  magnetic  field,  anti-Larmor  circularly 
polarized  rf  magnetic  field,  and  arbitrary 
orientation  between  the  rf  magnetic  field 
and  the  static  magnetization.  .The  experi¬ 
mental  data  and  the  extended  theory  give 
good  agreement  for  saturated  samples.  For 
thin  disc  geometries  it  has  been  found  that 
the  threshold  is  very  critical  to  align¬ 
ment  vhen  the  dc  magnetic  field  is  perpendl 
cular  to  the  plane  of  the  disc.  If  a 
needle-like  anti-parallel  domain  structure 
is  assumed  for  the  demagnetized  state, 
then  the  experimental  results  for  the 
demagnetized  state  can  be  explained  in 
terms  of  the  results  of  the  saturated 
state . 
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